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1. - Introduction 

The transition from coloured partons (quarks and gluons) to colourless hadrons - the so- 
called hadronisation process - is an exemplary process of the fundamental theory of the strong 
interaction, Quantum Chromo-Dynamics (QCD), which still lacks a quantitative understanding 
from first principle calculations. The process by which a highly virtual parton radiates gluons or 
splits into a quark- antiquark pair can be theoretically described by QCD evolution equations such 
as the DGLAP (Dokshitzer-Gribov-Lipatov-AltarelU-Parisi) equations [1-3]. However, the final 
"bleaching" of partons into hadrons takes place at a low virtualities (Q « Aqcd ~ 0.2 GeV) and 
so is dominated by nonperturbative QCD effects which cannot be theoretically addressed with 
the existing perturbative techniques. Modeling and phenomenology - e.g. as implemented in the 
Lund string [4] or cluster fragmentation [5] approaches - are often used to describe hadronisation 
processes. 

One way to study fragmentation and hadronisation is to perturb the environment surrounding 
the hard-scattered parton by introducing a nuclear medium [6]. The nuclear medium provides 
a sensitive probe of parton evolution through the influence of initial-state (IS) and/or final-state 
(FS) interactions. Such IS and FS may result on modifications of the final hadron yield distribu- 
tions compared to "vacuum" production and can help us understand for example the time-scale 
of the hadronisation process [7]. 

Nuclear modifications of hadron production have been indeed observed in Deep Inelastic 
lepton-nucleus Scattering nDIS (^* H- A), in hadron-nucleus (h -H A) and in heavy-ion (A -H A) col- 
lisions, compared to "elementary" DIS on a proton target or proton-proton collisions. In nDIS 
and h-\- A collisions, the medium is the nuclear target itself, also called "cold QCD matter". In 
A + A reactions, the produced parton must in addition traverse the created hot and dense medium 
("hot QCD matter"), be it a hadron gas at low temperature, or a Quark-Gluon Plasma (QGP) 
at high temperatures. In all cases, at high enough pr where hadrons mostly come from parton 



4 



ALBERTO ACCARDI ETC. 



fragmentation, one typically observes two different phenomena: (i) a suppression of hadron mul- 
tiplicities, called hadron or jet quenching, and (ii) a broadening of hadron transverse momentum 
spectra, which induces a local enhancement of the hadron pr spectrum known as "Cronin ef- 
fect" [8, 9]. Such nuclear effects are due to elastic and inelastic interactions of the incoming or 
outgoing partons and/or of the produced hadrons while traversing the surrounding medium. 

In nDIS, the target nucleus allows one to test the hadronisation mechanism and colour con- 
finement dynamics in a clean environment. Knowledge of partonic in-medium propagation 
gained from nDIS can be used in Drell-Yan (DY) lepton pair production in h + A colUsions 
to factor out FS effects (such as medium-induced gluon radiation) from IS effects (such as nu- 
clear modifications of parton distributions). A precise knowledge of parton propagation and 
hadronisation mechanisms obtained from nuclear DIS and DY studies can be very useful for e.g. 
testing and calibrating theoretical tools used to determine the properties of the QGP produced 
in high-energy heavy-ion collisions, as well as to reduce the systematic uncertainties in neutrino 
experiments with nuclear targets. 

This review is structured as follows. For a non technical overview of motivations, lessons 
learned, and an outlook, the remainder of this introduction (where we discuss hadronisation in 
elementary collisions, cold and hot QCD matter) can be followed by a reading of the concluding 
Section 9. The intervening sections detail the experimental and theoretical status, and discuss 
open issues and future experimental possibilities. In Section 2, we define the relevant observ- 
ables and kinematic variables for nDIS and hadronic collisions, comparing the phase-spaces for 
hadron production in these two cases; and we discuss the space-time development of hadroni- 
sation, introducing various estimates of the hadron formation time. In Sections 3, 4, and 5, we 
review the most relevant experimental results. Specifically, in Section 3 hadron production in 
v-nucleus, yu-nucleus, and e-nucleus DIS with emphasis on recent HERMES and CLAS data 
are discussed. Sections 4 and 5 are devoted to high-z^j- hadron production in h + h, h+ A and 
A + A collisions, focusing on recent results from the Relativistic Heavy Ion Collider (RHIC). 
The existing theoretical frameworks for interpreting the experimental data based on partonic or 
hadronic degrees of freedom, are discussed respectively in Sections 6 and 7. Specifically, parton 
propagation and energy loss in both cold and hot QCD matter is addressed in Section 6, where 
data from DIS, DY, h + A and A + A collisions are confronted to different theoretical models. 
In Section 7, the interaction of the prehadronic system and of the formed hadron with the nu- 
clear medium is discussed. Finally, in Section 8, we discuss observables sensitive to the time 
scales and different mechanisms involved in parton propagation and hadronisation, along with 
experimental measurements at future facilities that would help to clarify the whole picture. 

I'l. Parton fragmentation in elementary collisions . 

1 1 perturbative Quantum Chromodynamics (pQCD), collinear factorisation theorems [10] 
allow one to explicitly separate the short and long distance QCD dynamics involved in the mech- 
anism of hadron production from parton fragmentation. In a general Hi -I- H2 inelastic colUsion, 
one writes the inclusive hadronic cross sections for production of a hadron h at large momentum 
transfer or "hard" scale Q^, as 



o-hard(Hi -I- H2 h + X) 



(1) 



hard 



fhi={q,g) Vi=0,N 



where ® denotes a convolution over the kinematical internal variables of the process and A'^ is 
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Fig. 1. - Illustration of universality of PDFs (<^//jv) and FFs (Dj^i,) in leading order processes. Clockwise 
from top left: + annihilation, Deep Inelastic Scattering (DIS), lepton pair (Drell-Yan) emission, and 
hadron production in hadronic collisions. Solid lines indicate leptons, dashed lines quarks. The small black 
disc represents the perturbatively calculable hard interaction coefficient Hhard- 



the number of hadrons in the initial state. In Eq. (1), is the typical hard scale of the pro- 
cess and i/haid is the short-distance and perturbatively calculable hard coefficient function for 
the {fi} j + X partonic process. The long-distance dynamics is factorised into (a) the Parton 
Distribution Functions (PDF) 4>f,\H,{x^i)^ which can be interpreted as the probability of finding 
a parton of flavour fi and momentum fraction x, inside the projectile {H\) and/or target {H2) 
hadron, and (b) the Fragmentation Function (FF) Dj^h, which gives the equivalent "probability" 
that the parton j fragments into the observed hadron h with fractional momentum z- These func- 
tions are non-perturbative and need to be extracted from experimental data. Typically, PDFs are 
extracted from "global QCD fits" of inclusive hadron production in lepton-nucleon DIS [A^ = 1 
and Dj^hiz; Q^) = 1 in Eq. (1)] At large hadron fractional momenta(') z - Phadron/Ppaiton <: 0.1, 
the FFs obey DGLAP evolution equations and are obtained from electron-positron annihilation 
into hadrons [A^ = in Eq. (I)]. The obtained PDFs and FFs are provided by various authors, 
e.g., CTEQ6.6, MRST/MSTW [12-14], and HKNS, DSS, AKK08 [15-19] respectively, to men- 
tion the most recent sets. Once they are known at a given scale their value at any other scale 
can be perturbatively computed by means of the DGLAP evolution equations [1-3]. The factori- 
sation scales and entering PDFs and FFs, as well as the renormalisation scale in the 
perturbative cross section, should be ^ [O^^ i" order to avoid large logarithmic corrections. 
An important consequence of factorisation theorems is that PDFs and FFs are universal i.e. 



(') At small z, successful QCD resummation techniques (e.g. the Modified Leading Logarithmic Ap- 
proximation, MLLA [11]) have been also developed to describe the evolution of a highly-virtual time-like 
partons into final hadrons. 
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process-independent. The measured FFs ine'^ +e~ h + X and the PDFs ine"^ +p h + X can 
then be used to compute observables in any other process, e.g., hadron spectra in proton-proton 
collisions [A^ = 2 in Eq. (1)], or Drell-Yan production of lepton pairs [A' = 2 and /' = i£ in 
Eq. (1)], see Fig. 1. When deahng with hadron production with nuclear systems, universality 
is however experimentally observed to breakdown: the details of the hadron production cross 
sections depend on the colUsion process that yields the final particles, as we discuss next. 

1'2. Parton propagation and fragmentation in cold and hot QCD matter. 

I lie basic assumption behind the factorised form of Eq. (1) is that the characteristic time of 
the parton-parton interaction is much shorter than any long-distance interaction occurring before 
(among partons belonging to the same PDF) or after (during the evolution of the struck partons 
into their hadronic final-state) the hard collision itself. In that case, one can treat each nucleus 
as a collection of free partons, i.e., in the absence of initial-state efi^ects the parton density in a 
nucleus with mass number A is expected to be simply equivalent to that of a superposition of A 
independent nucleons(^); (pajA^x, Q^) - A ■ (f>a/Nix, Q^). In addition, in the absence of final-state 
effects the parton fragments with universal FFs and, therefore, the pQCD factorisation theorem 
for colhsions involving nuclei A predicts that minimum-bias inclusive hard cross sections scale 
respectively as 

do-hardiU h+A^h+X) - A dcThardiU h + p ^ h + X) , 

(2) , 

do-hard{A +A^h + X) = A"- do-hardiP + P ^ h + X) . 

The cleanest environment to test the validity of Eqs. (2) and study possible nuclear modifi- 
cations of hadron production (i.e. "violations" of the expectations given by Eqs. (2)) is nuclear 
Deep Inelastic Scattering (nDIS). In nDIS processes one experimentally controls many kinematic 
variables; the nuclear medium (i.e., the nucleus itself) is well known; and the particle mLiltiplicity 
in the final state is low, leading to precise measurements. The nucleons act as femtometer-scale 
detectors of the scattered hadronising quark, allowing one to study its space-time evolution into 
the observed hadrons (Fig. 2, left). The relevant observable in semi-inclusive nDIS processes 
is the ratio of the single hadron multiplicity on a target of mass number A normalised to the 
multiplicity on a deuteron target. At leading order, this multiplicity ratio corresponds to good 
approximation to the ratio of fragmentation functions (FF) in cold nuclear matter (the nucleus 
A) over that in the "vacuum" (deuteron). Recent HERMES measurements show that this ratio is 
significantly below 1 clearly showing a breakdown of universaUty for fragmentation functions, 
see Section 3 and Refs. [21-24]. 

A complementary means to study parton propagation in cold QCD matter is by measuring 
the Drell-Yan (DY) process in hadron-nucleus collisions: h + A ^ (*f + X, where hadroni- 
sation does not play a role (Fig. 2, centre). If the invariant mass of the lepton pair is large, the 
process can be described perturbatively as a parton-parton scattering producing a virtual photon 
which subsequently decays into the lepton pair. Any modification of this process will come from 
initial-state nuclear interactions of the projectile parton inside the target (as well as from nuclear 
modifications of the PDF, which can be isolated by other means). 

Tests of pQCD factorisation in hot-dense QCD matter can be carried out studying high-/7r 
hadron production in head-on nucleus-nucleus reactions (Fig. 2, right). The suppression of large 



(^) In reality, nuclear PDFs are modified compared to proton PDFs by initial-state "(anti)shadowing" ef- 
fects (see [20] for a recent review). 
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Fig. 2. - Quark propagation inside a target nucleus ("cold QCD matter") in lepton-nucleus {left) and hadron- 
nucleus Drell-Yan (centre) collisions. Right: Hard scattered parton traveling through the "hot QCD 
matter" produced in a nucleus-nucleus colUsion. 



transverse momentum hadron production in A + A compared to proton-proton {p + p) and hadron- 
nucleus h + A collisions at RHIC [25-28] (see Sect. 5), is also indicative of a breakdown of the 
universality of the fragmentation process. The standard explanation is that the observed sup- 
pression is due to parton energy loss in the strongly interacting matter. This assumes of course 
that the quenched light-quarks and gluons are long-lived enough to traverse the medium before 
hadronising, which can be expected at large enough pj because of the Lorentz boost of the hadro- 
nisation time scales. However, dynamical effects may alter this argument (see, e.g., Ref. [29]), 
with hadronisation starting at the nuclear radius scale or before. In this case, in-medium hadron 
interactions should also be accounted for, possibly leading to a different suppression pattern. 
Such mechanisms may be especially important in the case of heavy (charm, bottom) quarks 
which - being slower than light-quarks or gluons - can fragment into D or B mesons still inside 
the plasma [30]. 

In summary, a precise knowledge of parton propagation and hadronisation mechanisms can 
be obtained from nDIS and DY data, allowing one to test the hadronisation mechanism and 
colour confinement dynamics. In addition, such cold QCD matter data are essential for testing 
and calibrating our theoretical tools, and to determine the (thermo)dynamical properties of the 
QGP produced in high-energy nuclear interactions. 

13. Hadronisation and colour confinement. 

Willie not having a direct bearing on the traditional topics of confinement such as the hadron 
spectrum, the hadronisation process nonetheless contains elements that are central to the heart 
of colour confinement, as already emphasised 30 years ago by Bjorken [6]. For instance, in the 
DIS process, a quark is briefly liberated from being associated with any specific hadron while 
traveling as a "free" particle, and it is the mechanisms involved in hadron formation that en- 
forces the colour charge neutrality and confinement into the final state hadron. The dynamic 
mechanism leading to colour neutralisation, which is only implicitly assumed in the traditional 
treatments of confinement based on potential models [31] or lattice QCD [32], can be studied 
quantitatively using the theoretical and experimental techniques discussed in this review. As an 
example, the lifetime of the freely propagating quark may be inferred experimentally from the 
nuclear modification of hadron production on cold nuclei, which act as "detectors" of the hadro- 
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nisation process. Finally, as already discussed, the behaviour of partons propagating through the 
medium created in high energy heavy-ion collisions can give tomographically insight into the 
properties of large-scale deconfined QCD matter (i.e., of the Quark-Gluon Plasma). While still 
at an early stage, the understanding of such elements will ultimately provide deeper insights into 
the confinement-related properties of QCD. 

1'4. Hadronisation and neutrino oscillations. - 

Neutrino oscillation experiments use nuclear targets to enhance the neutrino detection rate. 
Nuclear effects change the topology and total energy of the experimentally measured hadronic 
final-state and are known to be one of the largest sources of systematic errors in current analyses. 
For the lower energy oscillation experiments that use the quasi-elastic channel, several prob- 
lems arise among which the distortion of the knock-out nucleon due to final state interactions, 
the contamination from hadron resonances, e.g. the tt H- — > A process by which a final-state 
pion is absorbed in the nucleus, and the unexplained depletion of low-virtuality events, much 
stronger than Pauli blocking can account for [33]. Experiments such as MINOS measure the 
neutrino energy adding up the muon and hadronic energies, Ey = + Ehad- Experiments such as 
OPERA [34] need to estimate the background to r-neutrino appearance events due to charmed 
mesons production and decay. It is thus crucial to have a good understanding of hadron modifica- 
tions in the nuclear medium and of the space-time evolution of the hadronisation process [35,36]. 
However, at the low hadronic invariant mass involved in these experiments, the theoretical meth- 
ods discussed in this review should be supplemented by those described in Refs. [37-40]. 

2. - Kinematics, observables, and hadronisation time estimates 

In this Section, we cover background material which will be used throughout the review. We 
define the kinematic variables and relevant observables for nDIS and hadronic collisions, and 
compare the phase-space for hadron production in both types of collisions. Finally, the space- 
time development of the hadronisation process is discussed. 

21. Kinematic variables. - 

We discuss the kinematics for hadron production in nDIS, hadron-hadron (h + h), hadron- 
nucleus {h + A) and nucleus-nucleus (A -i- A) collisions, and for the Drell-Yan (DY) process. We 
will explicitly make reference to the leading order (LO) processes in perturbative QCD, but most 
of the definitions are of general nature. 

Throughout this discussion we use light-cone coordinates = (p^ , , pr), where = 
(/j" ± p^)/ V2 and pr = (p^, p^). Our reference frame is such that the z axis is aligned with the 
beam, and a particle moving in the positive z direction has large light-cone plus-momentum. The 
transverse plane is the plane transverse to the beam and p'^p~ = + p\ = nij, where irp- = p^ 
is the invariant mass squared of the particle, and its transverse mass. 

2' 1.1. Deep inelastic collisions. - Deep inelastic scattering at LO in pQCD proceeds by ex- 
change of a virtual photon in the f-channel (Fig. 3 left). The DIS Lorentz invariants are de- 
fined in Table I. Note that the variable xb, and v are not independent but related through 
xb = 2^/(2Mv) in any reference frame. Analysis of inclusive DIS is usually carried out using 
Xb and Q^, because of the xs-scaling of the total cross section in the Bjorken limit: oo, 
Xb fixed. Note that in DIS one can experimentally measure all the listed variables, especially v, 

and Zh, because the initial and final state electron is observable. This is markedly different 
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Variable 


Definition 


Target rest frame form 












Xb 




2Pq 






Q' 










V 




qP 






y 




qP 
kP 




V 

W 






{P + qf 




hf- + 2Mv - 2' 


Zh 




Ph-P 

qP 




V 


Pt 




\Pt\ 







Target mass 

Bjorken scaling variable 

Negative four-momentum squared of 
the exchanged virtual photon 
Energy of the virtual photon in the 
target rest frame 

Fractional energy loss of the incident 
lepton (inelasticity) 
Invariant mass squared of the hadronic 
final state 

Fraction of the virtual photon energy 
carried by the hadron 
Hadron transverse momentum (relative 
to the virtual photon momentum) 



Table I. - Definitions of the kinematic variables for semi-inclusive DIS. The Lorentz invariant definition and 

its form in the target rest frame are provided. Particle 4-momenta are defined in Fig. 3 and 4. All variables 
are experimentally measurable, hence typeset in boldface. Note that Xb - Q^/(2Mv) independently of the 
chosen reference frame. 



from the situation in hadronic collisions, where only final state hadrons can be observed and not 
the partons themselves. The hadron transverse momentum in DIS is defined with respect to the 
virtual photon direction, see Fig. 4. Its analog in hadronic collisions would be the transverse 
momentum of a hadron with respect to the beam axis. 

Nuclear DIS experiments have been performed in fixed-target (ft) conditions in facilities like 
Stanford Linear Accelerator Center - SLAC (E665), Super Proton Synchrotron - SPS (EMC), 
Deutsches Elektronen Synchrotron - DESY (HERMES), Jefferson Lab - JLab (CLAS); and are 
planned in collider mode (cZ) e.g. at the proposed Electron-Ion Collider - EIC or Large Hadron- 
electron Collider - LHeC. The colhding nucleon and lepton momenta are 



(3) 

Prl = 



^ , V2£;v, Or) , kci = {y/2E, 0, Or) 



where E and Et^ are the lepton and nucleon energies measured in the laboratory frame. To 
discuss both modes at the same time, it is convenient to introduce the target rest frame energy of 
the lepton, Etrf 



I E fixed-target 
^"^^ ^"-/-li^ collider mode 

M 

Then the invariant inelasticity y for both modes can be written asy = v/E,rf. 



2' 1 .2. Hadron-hadron collisions. - Parton production in hadronic collisions at leading order 
in the coupUng constant a, proceeds through 2 — > 2 partonic collisions (see Fig. 3 right and 
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k' 




Ph Ph 



Fig. 3. - LO kinematics for parton production in DIS collisions (left) and in hadron-hadron collisions (right). 
Double lines indicate hadrons or nuclei, thin single lines are partons or leptons. The labels define the 
particles 4-momenta. 



Table II). Several LO processes can contribute to a given ij — > /1/2 collision, represented by a 
black dot in the cartoon, see Ref. [41] for details. The momenta of the two nucleons colliding in 
the centre-of-mass-frame (c.m.f.) with energy ^fs/2 each are 



(5) 



2 V2S 



,07 



[ V2S' 



where M is the nucleon mass and 



(6) 



s = s- 



1 + ^/TT^F/s^ 



We will neglect terms of order ff(M^/s) compared to terms of and will use s ^ s. In 

Eq. (5), we explicitly retain the nucleon mass to be able to perform boosts to the rest frame of 
either nucleon. If we assume the partons to be massless and collinear to their parent nucleons, 
their 4-momenta in terms of the parton fractional momenta jc, read / - (^x^ -\/s/2, 0, Oj-) and 

; = (0,X2 V^/2,0r). 

Particle production (partons or hadrons) is described in terms of the particle rapidity and 
transverse momentum. The rapidity of a particle of 4-momentum p and mass squared - 
is defined as 



(7) 



1, Ip' 
xiog — 
2 \p 



= log 



Positive rapidity describes a particle moving in the positive z direction, and likewise for negative 
rapidity. In the non-relativistic limit the rapidity coincides with the particle longitudinal velocity 
/3l in units of the speed of light, y — > Pl. Given the rapidity, one can compute = mr coshy 
and p^ = niT sinhy. Note that under a longitudinal boost of velocity /3, the rapidity transforms 
additively: y' = y - yp, where yp - 0.51og[(l + y6)/(l - jS)] is the rapidity of the particle rest 
frame. As an example, let us boost the target hadron momentum P from the target rest frame, P - 
{Ml V2, Ml y/2, Ot"), to a frame in which P' - ( ^/sJ2, M^l yfls. Or). This boost is accomplished 
by Of = ^fs|M. Likewise, boosting a nucleon from energy V^/2to V//2 requires a - yfVJs. 

Measuring the rapidity of a particle requires measuring two independent variables, say, its 
energy and longitudinal momentum. Not in all experiments this is possible, while just measuring 
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Fig. 4. - Kinematic planes for liadron production in semi-inclusive deep-inelastic scattering and definitions 
of the relevant lepton and hadron variables. The quantities k (k') and E {£') are the 4-momentum and the 
energy of the incident (scattered) lepton; ph is the 4-momentum of the produced hadron, and its transverse 
component relative to the lepton plane is denoted by pr- 



the polar angle 0* between the particle trajectory and the beam axis in the centre-of-mass frame 
is easier. This justifies the definition of the particle pseudorapidity, 



(8) 



T] = -logtan(6l'/2) , 



such that \p\ = Pt cosh 77 and = pr sinh 77. For massless particles it coincides with the rapidity: 
T) = y; for massive particles, they are approximately equal if \p\ » m (and 9 not too small). 
Differential particle distributions in y and tj are related by 



(9) 



dN dN 



^ _ dN E 

dydpj dj]dpj V nil cosh^ y dr]dp\ nij 



In order to compare collider and fixed-target experiments, and difi^erent beam energies, it is useful 
to consider the rapidity in the c.m.f.: 



(10) 



JcmJ. = y- Jem ■ 



The backward rapidity region (target hemisphere) corresponds to y c.m.f. < 0, and the forward 
rapidity region (projectile hemisphere) to yc.m.f. > 0. 

At LO, the 4-momenta of the two produced partons can be expressed in terms of their final 
state rapidities y, and transverse momentum pj (see Table II for definitions) 



(11) 



n 1 Pt v2 Pt -v, , 



and their fractional momenta are 
(12) 
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VaHdble 


Definition 




s 




Nucleon-nucleon centre-of-mass energy squared 


Xi 


= r/r 


Initial-state projectile parton fractional momentum 


X2 


= r/j- 


Initial-state target parton fractional momentum 


PiT 




Final state partons transverse momentum (relative to beam) 


yt -- 


= o.5iog(/:v/r) 


Final state partons rapidity 


ycm 


= 0.5 log (^) 


Rapidity of the centre-of-mass 


Z = 


= pUft 


Hadron fractional momentum relative to parent parton /i 


PhT = 


-- IPhtI 


Hadron transverse momentum (relative to beam) 


yh 


-- 0.51og(p,Vpp 


Hadron rapidity 


u -- 


-logtan(6>*/2) 


Hadron pseudorapidity (6' is the polar angle between the parton and 






the beam in the centre-of-mass reference frame) 



Table II. - Relevant kinematic variables for semi-inclusive parton (top) and hadron (bottom) production. 

Particle 4-momenta are defined in Fig. 3, right. Boldface variables are experimentally measurable; the 
others are theoretically defined in perturbative QCD. 



Finally, the Mandelstam invariants are defined as follows, 

S = {i + jf 

(13) t = {i-hf = {f2-jf 

= (J -/2)' = (/i -;•)'. 

and s + {+ u = Ohy momentum conservation. In terms of rapidities and transverse momentum, 
the Mandelstam invariants read 

S = X\X2S 

(14) i = -p].{l + e^^-^') 

u = -p].{l + e^'-^^) . 

Hadronisation in the coUinear factorisation framework proceeds through independent parton 
fragmentation into a hadron. It is universal, i.e., independent of the process which produced the 
fragmenting hadron, e.g., hadronic or DIS collisions [42J. The hadron fractional momentum z is 
defined by 

(15) pI = Zf^ PhT = zflT ■ 

Therefore the on- shell hadron momentum ph reads 

(16) Ph = \zf^, \^^!'\ zhTy 

The parton and hadron rapidities are related by y\ =yh+ \og{mhT I Pht)- 

The partonic variables pT, y,, x, and the fractional hadron momentum z are not experimentally 
measurable, but are needed in the theoretical computation of the cross section. The experimen- 
tally measurable variables are typed in boldface in Table II. Note that the hadron transverse 
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Variable 


Definition 


LO in OTj 




Xf 
M 

Pt - 


2pl/yfs 
\Pt\ 


= Xi — X2 
= X1X2S 


Feynman x 

Dilepton invariant mass 
Dilepton transverse momentum 



Table III. - Kinematic variables for Drell-Yan dilepton production. The dilepton momentum is p = p{*+ P(-, 
where p(± are the lepton and anti-lepton momenta. A star indicates momenta measured in the centre-of- 
mass frame. The three DY variables are experimentally measurable, hence typeset in boldface. See Table II 
for the definition ofxiy, and s. 



momentum pr in hadron-hadron collisions is defined with respect to the beam axis, so that at 
midrapidity it is the analog of the hadron energy £/, in DIS. 

2'1.3. Drell-Yan processes. - Drell-Yan production of a lepton pair in hadronic collisions 
occurs, at zero-th order in via the quark-antiquark annihilation channel 

as shown in the lower-right panel of Fig. 1. As before, the initial partons carry a momentum 

fraction x\ and X2 of the projectile and target hadron, respectively. At this order, the dilepton 
is produced with zero transverse momentum, with invariant mass M;+/- and with longitudinal 
momentum fraction Xf - 2p*/ ^/s, also called Feynman-x, where p* is the dilepton momentum 
in the centre-of-mass frame. These variables are summarised in Table III. By energy-momentum 
conservation and assuming s, we have 

= s = X1X2S 

(17) 

Xf — Xi — X2 i 

from which it follows immediately that 

(18) ^u = ^ (V^b+4M2/*±x,). 

Usually, Drell-Yan production is measured between the charmonium and bottomonium masses 
(4 < M < 8 GeV/c^) or above the Y. At higher orders in a,, new channels open up, such as 
Compton scattering, qg — » y q. The virtual photon acquires a finite transverse momentum, and 
it is no longer possible to relate xp and M to the momentum fractions of the partons probed in 
the projectile and target hadron. 

2'2. Comparison of hadron-hadron and DIS kinematics . - 

If we consider parton and hadron production at LO in hadronic and DIS collisions, it is easy 
to relate the relevant variables in both processes, that allows one to compare their corresponding 
phase spaces. The discussion closely follows Ref. [43], to which we refer for details. To connect 
the DIS and hadron-hadron kinematics (Fig. 3) we can boost the DIS collision to a frame in which 
the target has energy ^/s/2 per nucleon, and imagine the lepton to be a parton of a coUinear 
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SPS 


FNAL 


RHIC 


RHIC 


LHC 


[GeV] 


17.5 


27.4 


63 


200 


5500 


Ayi 


2.4 


2.0 


1.2 





-3.3 



Table IV. - Rapidity shifts Ayi of the RHlC-equivalent DIS phase space, tabulated for some energies of 
interest. 



phantom nucleon of energy ^/s/2 and with 4-momentum /"* = P"^. Comparing the left and right 
parts of Fig. 3 we can identify 

(19) P = J, P' = 1, k = i, k' = /2. 

The virtual photon momentum q, the fractional momentum Xg of the initial-state lepton and the 
rapidity ye of the final state lepton are identified as follows 

q = k-k' = i-f2, Xe = I P'* = xi, ye = y2. 



In this way, we can relate the DIS to the hadron-hadron kinematics discussed in Sect. 21. As 
an example, it is immediate to see that, in terms of hadron-hadron variables, - -t. The full 
translation "dictionary" from DIS to hadron-hadron variables can be obtained in a straightfor- 
ward way by combining the results of Sect. 2' 1 and the definitions of Tables I-II. 

First, we can express the DIS invariants in terms of parton rapidities and transverse momenta. 
Neglecting target-mass corrections, i.e. up to terms of ^{M^/s), we obtain 

Xb = ^ (e-y^ + e-y^) 
^fs 



(20) 



2M 
1 

y 



1 -I- 

Zh=Z. 

Note that the first three variables are not independent because = 2Mxbv, and that xb - X2 is 
interpreted as the struck parton fractional momentum, as expected in DIS at LO. Note also that v 
increases with increasing pr and increasing y\ . This is because a parton of positive and large yi 
in the cm. frame travels in the opposite direction as the left-moving nucleus, considered as the 
"target nucleus" (see Fig. 3). Hence in that nucleus rest frame it is very fast. Conversely, a parton 
of negative and large ji travels in the same direction as the target nucleus, which means quite 
slow in the target rest frame. It is also interesting to note that up to terms of order ^(M^/s), the 
parton and hadron energy in the target rest frame are £ = v and Eh = ZhV respectively. Finally, 
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Fig. 5. - Left: RHIC-equivalent phase space of nuclear DIS experiments at = 21 .b GeV (HERMES, 
solid line), at E, = 12 GeV (HERMES and JLab, dashed line), and at E, = 280 GeV (EMC, dot-dashed 
line). The dotted line shows the borders of the LO pQCD phase space in p, A + A at top RHIC energy, 
^jSNN = 200 GeV. The two arrows show the location of the midrapidity region at SPS and FNAL (p, A + A) 
fixed-target experiments. The shaded regions show the region of phase-space experimentally explored at 
HERMES [23,44] and EMC [45]. Right: Hadron-hadron-equivalent EMC and COMPASS (f + A) phase 
space at ^fsmi = 27.4 GeV, compared to the SPS and FNAL (p, A + A) phase spaces. 



we can invert Eq. (20) to obtain the hadron-hadron variables in terms of DIS invariants: 



■log 



(21) 



^2 = yi + log 

Z^Zh 



2ME,rf 



y 



with y - v/Etrf. Note that in DIS, the electron energy E,,/, hence the electron Xg, is fixed by 
the experimental conditions; this is different from hadronic collisions where the parton j has an 
unconstrained fractional momentum. Changing the c.m.f. energy to V^^ simply results in a shift 
of the parton rapidity, 



(22) 



yi 



200 GeV 



where Ayi - log( -\fs/ 'sfF). The value of Ay] compared to RHIC top energy 
is listed in Table IV for the experiments of interest in this paper. 

Given a DIS phase space (say a given experiment acceptance region in the (v, Q^) plane), we 
define its hadron-hadron-equivalent phase-space as its image in the (pT,yi) under Eqs. (21). The 
reason for this definition is that for both hadronic and DIS collisions we can identify the parton 
/i of Fig. 3 with the "observed" parton in hadronic and DIS collisions, i.e., the parton which 
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fragments into the observed hadron. Then the variables pr and yi fully characterise the observed 
parton. An analogous definition holds when using xb instead of v as independent variable. As an 
example, the HERMES DIS phase space in the (v, (^) plane is determined by the values of W^,„, 
2mm and j™,;^: 



n2. + jy2 

^min mm . y-i 

<V<ymaxE,rf 



(23) 2M 

Qli„<Q^<M^ + 2My-Wl,„ 



Additionally, one may impose stronger cuts on v, e.g., v > Vmm, as at the EMC experiment, and 
in some HERMES analysis. 

Using Eqs. (21) the hadron-hadwnequivalentDlS phase space in the (yi, pr) plane can be de- 
termined. As an example, in Fig. 5 left, we consider the RHIC-equivalent phase space of the fixed 
target e + A experiments and the planned Electron-Ion Collider (EIC), using ■\/Jm = 200 GeV. 
Note that according to Eq. (22), the hadron-hadron-equivalent phase space at other centre-of- 
mass energies can be obtained by a shift >'i — > >'i + Aji , see Table IV and Fig. 5 right where the 
Fermilab-equivalent phase space is shown. We assume the pQCD formulae used to define the 
hadron-hadron-equivalent phase space to be valid for p > = GeV/c), see Eq. (25) below 
for details. We can see that the HERMES and CLAS experiments, withfifr/ = 27.6 and 12 GeV, 
cover less than one third of the available RHIC pt range at yi ^ -3, with shrinking pT coverage 
at larger rapidity. In the SPS/FNAL midrapidity region it reaches pr = 2.5 GeV/c at most. Since 



(24) >'i<log(— i 



the only way to effectively reach larger values of yi is to increase the electron beam energy 
Indeed, the EMC experiment, withE,^ = 100 - 280 GeV, covers a larger span in rapidity and 
extends to yi > (as would aji + A programme at COMR\SS). Moreover, the increased energy 
allows one in principle to reach much higher pj than at HERMES and CLAS. However, only 
the Pt ^3 GeV/c region has been explored in actual measurements. The proposed Electron-Ion 
Collider (EIC) [46, 47] will be able to efl'ectively study the >'i > region and cover most of the 
RHIC phase space, but only the ji < part of the LHC phase space. 

The direct use of HERMES and CLAS data to understand nuclear effects in high-pr hadron 
production in heavy-ions collisions is therefore not possible. Instead, one needs to use those data 
to understand in detail the hadronisation dynamics and to constrain the various models in the 
context of heavy-ion collisions, by extrapolating them to unmeasured regions of phase space. 

The DIS-equivalent hadron-hadron phase space is defined as the image of Eqs. (25) in the 
(v, 2^, y, Zh) space under Eqs. (20). Then, the hadron-hadron phase space at a given yi is defined 
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Fig. 6. - left: Fixed-yi hadron-hadron trajectories plotted in the DIS-equivalent (y, Q^) phase space for 
RHIC at yjjviv = 200 GeV and various rapidities, for FNAL and SPS at midrapidity. The dashed line 
encloses the HERMES phase space; the dotted line encloses the EMC phase space. The arrow indicates the 
direction of increasing {pj) and (z;,>. Right: Trajectories in the (v, Zh) plane. The arrows indicate increasing 
Pt and Q^. 



by the kinematic bounds on 2 — > 2 parton scatterings [48]: 



|yi|<cosh"' — 



/?() <Pt < 



(25) 



2 cosh(ji) 



log I ^ - e-y^ ] < 3,2 < log - e>'' 
, Pt I \ Pt 



1 + 



PhT 



< z < 1 



where hard scatterings satisfy pj > po, with pa > 1 GeV/c a lower cutoff [48,49]. 

Introduction of next-to-leading order kinematics [50], would relax somewhat these bounds. 
At large rapidity, where the phase space for 2 — > 2 parton processes is becoming more and more 
restricted, 2^1 parton fusion may become the dominant mechanism because it is sensitive 
to much lower fractional momenta .jc, [51]. Hence, at the boundary of the hadron-hadron phase 
space, the presented analysis becomes unreliable. 

This phase-space is 4-dimensional and difficult to directly visualise. A way around this 
problem is to define suitable trajectories in hadron-hadron phase space averaged over 3^2, and 
to project them into the DIS-equivalent {v,Q^) and (v,z;,) phase spaces. We can define a pkj- and 
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SPS 


FNAL 




RHIC 






^/sJ^ = 17.5 GeV 


V^iw = 27.4 GeV 




V^iw = 200 GeV 




yh 











-1 -2 


-3 


PhT 


1-8 


1-12 


1-90 


1-60 1-25 


1-9 



Table V. - Range of hadron transverse momentum (phr, in GeV/c) spanned along trajectories at fixed 
rapidity yi at RHIC top energy s/s^n = 200 GeV and at fixed-target energies sjSf/N = 17-28 GeV. 



jft-dependent average observable as follows 



Jdzdyi dy2 &{pt, yi,y2,z)- 



^^^^ ^^^""^-y = r.-... ... da^s^, - 



J dzdyidyi-^ 



■ dp].dyidy2dz 

where 

ni\ d&AB-.hX _ V ^ d&AB-.f,X 

dpjdy\dy2dz ff ^ dpjdyidy2 



d&AB^fiX is the LO pQCD differential cross-section for production of a f\ parton in a colli- 
sion of hadrons A and B (nucleons or nuclei), and is its fragmentation function into the 

observed hadron. Then we can define fixed-y/, trajectories {{{v)p^^y, {Q^)p.r,y)\ pr ^ Po) and 
{i{y)pT,y^ {zh)pT,y)'^ Pt ^ Po) in the DIS-equivalent phase space. 

As an example, in Fig. 6 we considered hadronic collisions at RHIC top energy = 200 GeV 
and at fixed- target energies s/snn = 17-27 GeV, and plotted the fixed-y;, trajectories in the DIS- 
equivalent phase space. The range of pr spanned along each trajectory is tabulated in Table V. 
The spanned range in is limited by the maximum pr at each rapidity, according to Eq. (25). 
As expected, the smaller the rapidity yi, ~ yi the smaller the spanned v. RHIC trajectories with 
yh < -2 span relatively low values of v < 60 GeV and large values of Zh ^ 0.5, where the EMC 
and HERMES experiments have shown non negligible cold QCD matter suppression of hadron 
production. At higher rapidity, the larger spanned values of v will make cold QCD matter effects 
less prominent. The consequences of these remarks for the interpretation of hadron production 
inh+A and A + A collisions will be further discussed in Section 63.4. 



23. Nuclear modification observables . 



Ill lepton-nucleus DIS, the experimental results for hadron production are usually presented 
in terms of the hadron multiplicity ratio R''^, which represents the ratio of the number of hadrons 
of type h produced per deep-inelastic scattering event on a nuclear target of mass A to that from 
a deuterium target (D). The multiplicity ratio depends on the leptonic variables v and Q^, and 
on the hadronic variables z = Eh/v and p^ defined in Section 21. It is defined as the super-ratio 

(28) R,iz, y, Q,Pr)- [ / [ ^^(,,g2) ' 

where A^;, is the yield of semi-inclusive hadrons in a given (z, v, Q^, p^)-bin, and Ne the yield of 
inclusive deep-inelastic scattering leptons in the same (v,(2^)-bin. Normahsing the hadron yield 
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to the DIS yield allows one to cancel, to a large extent, initial-state nuclear effects such as nuclear 
modifications of PDFs and to isolate final-state nuclear modifications of hadron production as a 
deviation of Rm from unity. A suppression of Rm is experimentally observed to increase with z, 
and to decrease with y, and more mildly with (see discussion in Section 3). When plotting 
Rm as a function of pt, one observes a suppression at small pt and an enhancement above 
/??• ~ 1-5 GeV/c. This behaviour is also known as "Cronin effect" (Sections 3'2 and 4'2). The 
amount of transverse momentum broadening (defined with respect to the direction of the virtual 
photon, see Fig. 4) is quantified via 

(29) A{p\f = {p\)\ - {p\% . 

Here, {p\i\ is the average transverse momentum squared of a hadron of type h produced on a 
nuclear target A 



(30) {p\f 



and {p\)'l) is the same quantity for a deuterium target. 

In high-energy proton-nucleus and nucleus-nucleus collisions what is usually presented is the 
nuclear modification ratio or ratio of the hadron transverse momentum spectrum measured in 
A -I- B at a given rapidity y and impact parameter b (or centrality class) normalised by the nuclear 
overlap function TAsib) - related to the "parton luminosity" atb- over the p+ p spectrum: 



^ ^ AB(/'r,3', ) jj^^Q)^ dpjdy j dpjdy 

TabOj) is computed with a geometrical Glauber eikonal model of the nucleus-nucleus collision 
(see e.g. [52]). In the absence of nuclear effects, one would expect Rab = 1- The behaviour of 
R'^gipT, y, b) as a function of pt will be discussed in Sections 4 and 5. 
For DreU-Yan processes, one defines an analogous nuclear ratio: 

1 do-"?}, I 1 d/crpY 
(jz; kdy - —■ 



A dMdxpdpj I B dMdxpdpj 



The dilepton pr-broadening is defined analogously to Eq. (29), with h = . Results for /?dy 
will be presented in Section 4' 1 . 

Medium modifications of hadron production in nDIS and heavy-ion collisions can also be 
revealed by means of multi-particle azimuthal correlations, which are sensitive to the underlying 
parton-medium interaction and to the properties of the medium. For example, two-hadron cor- 
relations measured at RHIC revealed significant broadening and softening of associated hadrons 
on the away side of a triggered high-pr particle, which is consistent with the observation of 
the hadron suppression in the single inclusive measurements. In nDIS processes, the distribu- 
tion of associated sub-leading hadrons to the leading hadron challenges various models of single 
inclusive hadron suppression. 
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tj, t,„ tl, 

Fig. 7. - Sketch of the time evolution of the hadronisation process with definition of various time scales. 
A quark q created at time in a hard collision turns into a coloured prehadron If,, which subsequently 
neutralises its colour, /j", and collapses on the wave function of the observed hadron h. Gluon radiation 
lasts until colour neutralisation. 



2'4. Hadronisation time estimates . - 

Even though hadronisation is a non perturbative process, a few features can be extracted from 
general grounds. A parton created in a high-energy collision can travel in the vacuum as a free 
particle only for a limited time because of colour confinement: it has to dress-up in a colour- 
field of loosely bound partons, which eventually will evolve into the observed hadron. The same 
dressing process can be expected for partons traveling in QCD matter, yet it will be modified by 
interactions with the surrounding medium. In a deconfined nuclear medium such as the QGP, the 
dressing process might furthermore be delayed until the medium cools down and comes closer 
to the confinement transition. 

The bare parton-medium cross section is dominated by the elastic a + b — > a + b parton- 
parton scattering and gluon bremsstrahlung. The gluon-gluon cross-section is of order cTgg - 
9/27ra^/fi^ ^ 2 mb (for a., = 0.5 at a pr-cutoff of order fJ. - I GeV), but the dressed parton is 
likely to develop an inelastic cross section of the order of the hadronic one, ff(40 mb), becoming 
subject to nuclear absorption similarly to the final state hadron. Hence, it can be viewed as a 
"prehadron" (denoted here by /i,). The prehadron may still be for a short time in a coloured state 
and radiate gluons neutralising its colour before its wave function collapses onto the observed 
hadron wave function. We can therefore identify three relevant time scales, see Fig. 7: (1) the 
"prehadron production time" or "quark lifetime" tpreh, at which the dressed quark develops an 
inelastic cross section, (2) the "colour neutralisation time" f„,, at which gluon bremsstrahlung 
stops, and (3) the "hadron formation time" t/,, at which the final hadron is formed. Typically, 
model applications further simplify the process and merge steps (2) and (3) assuming tpreh - ten- 
Note that the prehadron and the formation times are introduced as a phenomenological tool, 
rather than a well defined quantity, in order to distinguish between the stage in which the parton 
can be described as an asymptotically free particle and treated in pQCD from the stage in which 
colour confinement and non-perturbative interactions kick in and warrant a treatment in terms 
of different degrees of freedom. Such a phenomenology is well suited to the present status of 
the theoretical and experimental investigation, but will need to be substantiated or replaced by 
a more fundamental QCD description. Finally, note that strictly speaking the very question "is 
the prehadron formed within or without the medium?" is ill-posed: in quantum mechanics it can 
happen one way in the amplitude and the other in its complex conjugate, and the interference 
between the two may be non-negligible [53]. 

A simple estimate of the hadron formation time (f/,) can be obtained by defining it as the time 
for the struck partons to build up its colour field and to develop the hadronic wave function [11]. 
In the hadron rest frame this time is related to the hadron radius /?/,, and in the laboratory frame 
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Fig. 8. - Quark hadronisation by emission of one parton {left) and emission of one parton accompanied by 
many soft gluons {right). 



it is boosted by y = Ehlnih. 



(33) 



Eh 

(th) ccRi, 

mh 



In Table VI we show the hadron formation time estimated with Eq. (33) for a typical 7 GeV pion 
{Rh ~ 0.7 fm [54]) at the kinematics conditions found in HERMES (z ~ Zh ~ 0.5, £/, ~ z/,v 
7 GeV) and at RHIC mid-rapidity (z ~ 0.7, E,, * pj. * 7 GeV/c). We find {t^} ^ 35 fm » Ra, 
which points towards a long quark lifetime with hadron formation outside the medium. However, 
for the heavier kaons, rj, and protons (with radii Rfc ~ 0.6 fm [55], R,j ~ R„, Rp ~ 0.9 fm [56] 
resp.) we obtain much shorter formations times (Ik) ~ 8 fm, (f,,) » 9 fm and {tp} « 6 fm, 
which are comparable to the size of the medium. Heavy D and B mesons, with average radii 
Rd = 0.57 fm and Rg = 0.5 fm respectively [57], clearly fragment in-medium {{to.B) 5 2 fm). 

An estimate of the prehadron production time can be obtained by looking at hadronisation 
in light-cone coordinates [30]. Consider a relativistic on-shell quark of mass niq and plus- 
momentum p^, hadronising into a hadron of mass nih and 4 momentum - zp^ ■ Minimally 
prehadron formation, i.e., the formation of a colorless partonic object, proceeds by emission of 
an additional, (typically light) parton to carry away the initial state colour, see Fig. 8 left. The 
process in momentum space is 



(34) 



r 2 ^ 






— > 



zp 



2zp^ ' 



(l-z)p' 



2{\-z)p' 



where we imposed 4-momentum conservation. In time-ordered perturbation theory, the light- 
cone separation Ax^ between the initial and final state can be estimated by the uncertainty prin- 
ciple: 



(35) 



Ax+ 



1/A;,- = 



2z{\-z) p' 



+ {\-z)ml-z{l 



■ z) m\ 



where A/5 = p^, + p^^ - p^. Since Ax+ = (Af/, -i- Az/,)/ V2, and Az/, = {pqlEq)^th with pq (Eg) 
the momentum (energy) of the quark, the prehadron formation time is 



(36) 
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Fig. 9. - Left: Hadronisation in the Lund string model in the target rest frame. Right: Prehadron production 
time {tpreh), and hadron formation times (?/,>, for v = 10 GeV quarks fragmenting into pions with fractional 
momentum zi,, as computed within the Lund string model in Ref. [61]. 



This estimate should be used with care since actually prehadron formation is likely to be accom- 
panied by the emission of many soft gluons, see Fig. 8, and the system of emitted partons has 
an invariant mass Taking this into account, we should add an additional zm\ term at the de- 
nominator of Eq. (35), which would reduce the estimated formation time. In Table VI we show 
the hadron formation time from Eq. (36) for typical hadrons at HERMES and RHIC, where we 
assumed P ^ ^\cd- ^^^^ energy pre -pions are formed outside the medium {{tpre-n) ~ 25 fm), 
while pre-kaons and pre-77 {{tpre-K) ~ 6 - 9 fm) and pre-protons ({tpre-p} ~ 4 fm) have formation 
times comparable to the medium size, and heavy pre-mesons are produced rapidly inside the 
medium. 

At least for light mesons, these estimates can be used to justify the computation of hadron 
quenching in terms of parton-medium interactions alone, as done in radiative energy loss models 
[58-60], see Section 6. 

A successful non-perturbative model of hadronisation is the Lund string model [62], see 
Fig. 9 left. The confined colour field stretching from the struck quark to the rest of the nucleus is 
modeled as a string of tension Ksir ~ 1 GeV/fm, and spans the lightly shaded area in space-time. 
The prehadron formation point is identified with the qq pair production point C, which breaks 
the string in smaller pieces [63]. Hadrons are formed at points P, when a quark and an antiquark 
at the endpoint of a string fragment meet. The subscript / indicates the so-called rank of the 
produced (pre)hadron, counted from the right of the figure. Average prehadron production times 
can be analytically computed [61,63,64] and have the following general structure: 

(37) {tp,,H} = fiZh) (1 - Z;,) — , {th) = (tpreh) + — , 

^str ^str 
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Table VI. - Estimates of typical hadron formation times {th) and prehadron production times {tpreh) for 
pions, kaons, rt, protons, and D and B mesons at HERMES (zh ~ 0.5, v « 14 GeV) and at RHIC at 
mid-rapidity (Py ~ 7 GeV/c, z~ 0.7) obtained with Eqs. (33), (36), and (37). 



where v is the struck quark energy, the function f{zh) is a small correction of (tpreh), which can 
be computed analytically in the standard Lund model [61,63], and K^ir - (I GeV/fm) R^/RI with 
Rh the hadron radius is taken from Ref. [61]. The factor zhv can be understood as a Lorentz 
boost factor; the (1 - Zh) factor is due to energy conservation: a high-z^ hadron carries away 
an energy z/,v; the string remainder has a small energy e = (1 - Zh)y and cannot stretch farther 
than L = e/Kstr- Thus the string breaking occurs on a time-scale proportional to 1 - Zh- The 
resulting pion formation time scaled by a factor v/Kstr is plotted in Fig. 9 left. A typical pion 
produced at HERMES energies (i.e. with fractional energy Zh ~ 0.5 from a parent quark with 
energy v 14 GeV) has (tpreh) ~ 6 fm < Ra and (th) ^ 13 fm > Ra, with similar values at 
RHIC at mid-rapidity. Therefore, the final hadron is typically formed at the periphery or outside 
the nucleus so that its interaction with the medium is neghgible (see Table VI). However, the 
prehadron is formed inside and can start interacting with the nucleus. A detailed space-time 
analysis of hadronisation in the PYTHIA/JETSET Monte Carlo implementation of the Lund 
string model has been performed in [65], with similar conclusions regarding the magnitude of 
the pion prehadron production time. 

In Ref. [66-68] the formation of a leading hadron {zh ^ 0.5) is described in a pQCD model, 
see Fig. 10. The struck quark radiates a gluon. The gluon then splits into a qq pair, and the 
q recombines with the struck q to form the leading prehadron, which later on collapses on the 
hadron wave function. The cross-section can be computed from the modulus squared of the 
sum of the two gauge-invariant amplitudes shown in the figure. At Zh ~ 1, higher twist effects 
spoil this simple mechanism for hadronisation. The dipole model of Ref. [29] approximates the 
described cross-section as a convolution of a Gunion-Bertsch radiation cross section [69] with the 
gluon splitting plus quark recombination process, see Fig. 10 right. The prehadron is identified 
with the qq pair which includes the struck quark, and its production time is identified with the 
time at which the gluon becomes decoherent from the struck quark. Note that, strictly speaking, 
the resolved quark-gluon system may be in an octet state from the production time until gluon 
splitting occurs. However, in this model, gluon radiation is neglected during the octet stage. The 
model can compute the probability distribution in the prehadron production time, see Fig. 11, 
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Fig. 10. - Diagrams for leading hadron formation in ^ + A collisions: Gauge-invariant set [66] (left), and 
dipole-model approximation [29] (right). 



and the average {tpreh) is 

(38) (WJc^d-z/,)^ . 

Q 

The interpretation of the 1 - Zh factor is in terms of energy conservation: the longer the struck 
quark propagates, the larger its energy loss; hence to leave most of its energy to a z/, — » 1 hadron, 
the quark must be short lived. The scale is set by Kjip - Q^. At HERMES, with ~ 10 GeV/fm, 
this model obtains for pion production {tpreh) < 5 fm at Zh > 0.5: pre-pions are formed inside 
the medium. At RHIC, where cc p^, and z/,v « pj- for mid-rapidity hadrons, one obtains the 
counter-intuitive result that {tpreh) ^/pt- prehadrons are formed the quicker the higher their 
transverse momentum [29,70], typically inside the medium. 

In summary, given the model dependences of the theoretical estimates of the hadron forma- 
tion time - and hence of the length of the partonic phase of the hadronisation process - it is 
important to study the kinematic- and flavour-dependences of various hadron production pro- 
cesses through a careful analysis of experimental data and tests of phenomenological models, 
see the discussion Section 8. 

3. - Experimental results in lepton-nucleus deep inelastic scattering 

Deep inelastic lepton-nucleus scattering offers a direct way to study the hadronisation pro- 
cess which follows the hard scattering. In contrast to hadronic collisions, in nuclear DIS no 
deconvolution of the parton distributions of the projectile and target particle is needed, so that 
the experimental observables can be more directly related to the nuclear effects on the quark 
propagation and fragmentation. Moreover, in electron and muon experiments the incoming and 
outgoing leptons are detected and they provide an extra handle to determine the kinematical 
variables of the produced partons/hadrons. 

In the past, semi-inclusive leptoproduction of undifferentiated hadrons from nuclei has been 
studied with neutrinos at FNAL, CERN and Serpukhov; with electrons at SLAC [71], and at 
CERN and FNAL with high-energy muons by EMC [72] and by E665 [73] respectively. Re- 
cently, HERMES has reported more precise data [21-23] on the production of charged hadrons 
as well as, for the first time, various identified hadrons (n^^, n^, if', K^, K , p and p) in 
deep-inelastic positron scattering off nuclei. Finally, high-statistics data have been collected 
at JLab [74] with a 5.0 GeV electron beam on targets of carbon, aluminum, iron, tin, and lead 
at large luminosities (2 10^** cm^^s '), with detection of several identified hadrons, particularly 
TT^, 7T^, , K^, and A. In the next Sections the most significant experimental results for DIS of 
neutrino, muon, electron and positron beams on nuclei are presented. 

31. Hadron production in v-nucleus DIS . - Neutrinos can interact via the exchange of Z° 
(neutral current), or can turn into charged leptons via exchange, while at the energies of 
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Fig. 11. - Probability distribution of the prehadron production time fp^A = t obtained in the colour-dipole 
model of [29]. 



interest here electrons scatter primarily through the exchange of photons. One thus has at hand 
various methods to extract information on the parton production and propagation combining 
results from experiments encompassing all four (neutral and charged) exchanged bosons. In 
contrast to hadron-induced collisions, which preferentially scatter from the front hemisphere of 
a nucleus due to their strong interaction probability, neutrinos and electrons directly scatter with 
partons or nucleons inside the nuclear target because of their weak interaction probability. This 
means that the fragments of the struck parton in neutrino and electron scattering experiments 
should suffer the same final-state interactions. In general, the laboratory energy Ey of the in- 
coming (anti)neutrino cannot be determined directly since there is substantial energy that goes to 
undetected neutral particles. Corrections on the estimation of the neutrino energy reflect on most 
of the kinematic quantities and significantly increase the systematic uncertainties for neutrino 
experiments. 

Neutrino-induced hadron production was first studied in the bubble chamber at FNAL [76] 
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Fig. 12. - Ratio of the differential hadron multiphcity distribution in (anti)neutrino-./Ve over -He collisions 
as function of the transferred quark energy v for fast (zi, > 0.2) charged hadrons [75] compared to model 
predictions for "constituent" (tc) and "yo-yo" (ly) lengths (see text for details). 
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and later in the bubble chamber BEBC at CERN [77]. The main objective of these experiments 
was the study of global properties of hadron production in nuclei. Fast hadron production on 
nuclei was found to be attenuated as compared to that in hydrogen target and to depend both on 
the variables z and v [76]. 

The production of hadrons in charged-current (anti)neutrino interaction was studied with 
higher accuracy by the BEBC WA21/WA59 Collaborations [75]. While in the previous experi- 
ments the analysis was hmited to negative hadrons in order to exclude systematic uncertainties 
due to 'knock-on' protons, charged hadron production was studied in this experiment in terms 
of z/,-distribution normalised to the number of events. A small but significant reduction of fast 
hadron {zh > 0.2) production was found in a neon target as compared to that in a hydrogen target. 
This is shown in Fig. 12 where the ratio of the normalised distribution is presented as function 
of the transferred energy v - Ey - E^, where Zs^ is the laboratory energy of the muon beam. The 
data indicates a significant (10% - 20%) attenuation (/Jjve < 1) of fast charged hadron yields over 
the whole v- range with a stronger attenuation at low v and high Zh- The experimental results 
are compared to theoretical predictions of Ref. [63] in which two hadron formation lengths are 
considered: the constituent (Ic) and the "yo-yo" (ly) lengths. Specifically, the constituent length 
corresponds to the time after which the first constituent of the hadrons is formed, the yo-yo length 
corresponds to the time after which the quark and antiquark meet to form the color singlet [78] 
The /3,-model overshoots the value of ^jve, while a /c-model give a fair description of the data 
thus pointing out that significant interactions of the hadronising system start as early as at the 
constituent point. 

The role of the formation length for the description of hadroproduction in DIS of neutrino 
on nuclei has been investigated by the NOMAD experiment at CERN [79, 80]. The backward- 
going protons and n~ produced in charged current neutrino interaction have been compared with 
intra-nuclear cascade (INC) models. In these models the production of particles in kinematically 
forbidden regions can be seen as the result of multiple scattering and of interactions of secondary 
hadrons with the other nucleons while they propagate through the nucleus. Experimentally it 
has been observed that the cascade is restricted to slow particles only, while the fast ones do not 
re-interact inside the nucleus. A proposed explanation for this effect is that, since the formation 
time is proportional to the hadron energy (via the Lorentz time-dilation factor y = Elm), the 
INC process is restricted to slow hadrons which have formation lengths smaller than the nuclear 
radius [79]. 

The inclusive spectra of hadrons have been also measured with the aid of the SKAT propane- 
freon bubble chamber irradiated with a beam of 3 to 30 GeV neutrinos from the Serpukhov 
accelerator [81]. The ratio of the yields of charged hadrons in the subsamples of nuclear interac- 
tions Ba and of deuteron interactions Bd is shown in Fig. 13 as a function of Zh (left) and y (right) 
respectively. The left part of Fig. 13 displays the ratio measured at {A) = 28 and energies in the 
range 2 < v < 1 5 GeV, compared with the multiplicity ratio obtained from deep inelastic scatter- 
ing of positron on nitrogen nuclei (A = 14) at higher transfer energies of 7 < v < 24 GeV [21]. 

The suppression of the hadron yield appears more pronounced for the most energetic hadrons 
{zh > 0.6). The V dependence of the ratio for leading charged hadrons with Zh > 0.5 is shown in 
the right part of Fig. 13, which shows the SKAT data at (v) = 3.3 and 7.7 GeV, along with data 
from e*-^^N interactions in the region v > 8 GeV. The reduction of the yields occurs for hadrons 
produced in processes where the parent quark has the lowest energies. The neutrino data, within 
their larger uncertainties, show the same trend as observed with electron beams. 

3'2. Hadron production in e-nucleus ( SLAC) and fi-nucleus ( CERN) DIS. - In contrast to the 
neutrino scattering experiments, in electron- and muon-induced DIS the incoming and scattering 
leptons are detected and provide a well defined reference system for the measurements of the 
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Fig. 13. - Ratio of the yields of charged hadrons in neutrino DIS on nuclei A over deuteron D targets 
(closed circles) [81] and in e* DIS on nitrogen and deuteron nuclei (open circles) [21]. Left: z;, -distribution 
for hadrons produced in processes with transferred energy 2 < v < 15 GeV. Right: v-distribution for hadrons 
with Zh > 0.5. 



kinematics of the outgoing parton/hadrons. 

Electroproduction of hadrons from nuclei was studied for the first time at the SLAC using 
a 20.5 GeV/c electron beam incident on different targets (H, D, Be, C, Cu, Sri) in the late sev- 
enties [71]. A single arm spectrometer was used to measure the scattered electrons and the 
produced hadrons. The ratio between the number of single inclusive hadrons detected per gram 
per square centimeter per incident electron for nucleus to the analogous number for deuterium 
was measured. The results were presented in the original paper [71] as a function of the trans- 




Fig. 14. - Ratio of differential charged hadron multiplicities as function of Zh measured ine+A interactions 
at SLAC [71] and derived in Ref. [82]. The curves represent a fit to the data in the functional form {A/iy^^K 
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actions by EMC [45]. 
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in /i-A collisions at EMC [45]. 



verse momentum pr of the hadron in the photon-nucleon system and Zc.m. = PilPmax, where pi is 
the longitudinal hadron momentum and pmax - y/s/2. The ratio of differential charged hadron 
multiplicities as a function of Zh, derived in Ref. [82] from the values published in the original 
paper, is shown in Fig. 14. This experiment shows an attenuation of the electroproduced hadrons 
which clearly increases with the size of the target nucleus. In addition the nuclear attenuation 
increases for increasing hadronic momenta. This measurement performed prior to the discov- 
ery of the EMC effect, does not account for modifications of the nuclear PDFs and, thus, needs 
correction because a semi-inclusive cross section ratio was measured instead of the multiplicity 
ratio as defined in Eq. (28). An estimate of the correction for this effect based on Ref. [21] results 
in a 4% increase of the original ratio. 

After the pioneering measurements with muon beams performed at FNAL [83, 84], hadron 
production in DIS of muons was studied by the EMC-coUaboration first in the v range 50- 
140 GeV [72] and later on in the v range 20-220 GeV [45]. Their large single-arm spectrometer 



PARTON PROPAGATION AND FRAGMENTATION IN QCD MATTER 



29 



allowed this experiment to carry out a precise determination of the muon momentum with a good 
acceptance for forward produced hadrons. In addition, due to the simultaneous measurements of 
nuclear targets and deuterium, most of the systematics uncertainties cancel in this measurement. 

The ratio of the hadron multiphcity distribution as function of Zh is shown in Fig. 15 for 
different targets. The corresponding mean value of the multiphcity ratios and of the v variable 
are noted on the figure. For large nuclei (Cm, Sn)a small but distinct reduction of the fast hadron 
production compared to that of deuterium is observed, whereas for carbon the ratio is compatible 
with unity over the whole range in Zh- The left part of Fig. 16 shows the hadron multiplicity 
ratio as a function of v in two bins. The same variation with v is seen in all intervals, thus 
no trend in is observed. The ratios show a gradual decrease with decreasing v below 60 GeV, 
whereas they slowly approach unity for higher v. The depletion of the fast hadron multiphcity 
in muon interaction with heavy targets is only ~10% even in the low-v bin. The pr dependence 
was measured as shown in the right part of Fig. 16 for two v intervals. At high pj the ratio rises 
above unity in both v intervals. The observed trend is consistent with the Cronin enhancement 
reported in hadron-nucleus collisions [8, 91 (see Section 4). Since in lepton-nucleus, at variance 
with hadron-nucleus, collisions neither multiple scattering of the incident particle nor interaction 
of its constituents can contribute to the Cronin effect, the observed enhancement can be ascribed 
to rescattering effects in the final-state only. 

In conclusion, the results from the experiments performed at SLAC with electrons and at 
CERN with muons have shown that the multiplicity ratios mainly depend on the energy of the 
virtual photon v, and on the fraction Zh of this energy carried out by the final hadron. These 
experiments demonstrate that the nuclear effects on the ratios decrease with increasing v. This 
has been confirmed by the Fermilab experiment E665 [73] in DIS of 490 GeV muons off xenon 
and deuterium targets. No nuclear dependence in the z/j-distributions of the forward-produced 
hadrons was found in the y-range from 50 GeV to 500 GeV. Thus we can conclude from the 
performed measurements that the transfer energies v where nuclear effects are the largest are in 
the range from a few GeV to few tens of GeV. Such an energy dependence is easily understood 
within parton energy loss models: since the energy lost by quarks propagating through QCD 
media, e, is independent of their energy (i.e. v in this context) in the high-energy limit, the 
relative energy-loss e/v - which controls the amount of hadron suppression in nuclear DIS - 
vanish at high energy making the ratio tend to 1. In the hadron-oriented picture, the larger v 
(parton energy) the longer the formation time, which therefore decreases the amount of nuclear 
absorption. 

3'3. Hadron production in e-nucleus DIS at the HERMES experiment. - The influence of 
the nuclear medium on lepto-production of hadrons has been recently extensively studied by the 
HERMES experiment at DESY in semi-inclusive deep-inelastic scattering of 27.6 GeV positrons 
off deuterium, nitrogen, neon, krypton and xenon targets. The data were collected in the v-range 
3-23 GeV using high density gas targets internal to the positron storage ring. During these high- 
density runs HERA operated in a dedicated mode for the HERMES experiment. 

For the fist time the v-dependence of the multiplicity ratio has been measured for various 
identified hadrons (n^, n~, tt", A"^, K~, p and p) and neutral pions [22], as shown for a kryp- 
ton target in Fig. 17. The identification of charged pions, kaons, protons and antiprotons was 
accomplished using the information from the RICH detector [85], which replaced a threshold 
Cerenkov counter used in the previously reported measurements for the charged hadrons multi- 
plicity on '^A^ [21]. The corresponding z/,-dependences of R'Ij with v > 7 GeV are shown in the 
right part of Fig. 17. In the bottom panels the average values for and Zh or v are displayed for 
all the presented data. 

The results presented in Fig. 17 and the average values reported in Table VII for zh > 0.2, 
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Fig. 17. - Multiplicity ratios for identified pions, kaons, protons and antiprotons from a Kr target as a 
function of v for Zh > 0.2 (left), and as a function of z/, for v > 7 GeV (right). In the upper right panel the 
multiplicity ratio for identified pions from a ''^A' target are also shown. The closed (open) symbols represent 
the positive (negative) charge states, and the crosses represent n° mesons. In the bottom panels the average 
Zh and V values are displayed: pions and kaons (protons and antiprotons) are shown as closed (open) circles; 
the average values are indicated by the open stars referring to the right-hand scales. The inner (outer) 
error bars represent the statistical (total) uncertainties. Multiplicity ratios for negative kaons and antiprotons 
at the highest z^-bins are not displayed due to their poor statistical significance. 



show that the multiplicity ratios for positive, negative and neutral pions as well as for negative 
kaons are similar. However, for positive kaons is significantly larger. An even larger differ- 
ence is observed between protons and their antiparticles compared to the meson case. These dif- 
ferences in of positive and negative kaons, as well as those between protons and antiprotons, 
are still present at Zh > 0.5. This is shown in the last column of Table VII, where the average 7?^ 
values are reported for Zh > 0.5, i.e. when emphasising leading hadrons. In addition the Zh > 0.5 
range is most suitable to compare of mesons and baryons as this comparison is performed 
at the same average v as shown in the bottom right panel of Fig. 17. The difference observed in 
mesons and baryons multiplicity ratios resembles the anomalous baryon enhancement reported 
at intermediate pr's {pr ~ 1.5-8 GeV/c) in proton-nucleus and heavy-ion collisions [86, 87] 
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hadron 


</?^> (ZH > 0.2) 


«> (zi, > 0.5) 




0.775 ± 0.019 


0.712 + 0.023 


n~ 


0.770 ± 0.021 


0.731 ±0.031 




0.807 ± 0.022 


0.728 + 0.024 


K* 


0.880 ± 0.019 


0.766 ± 0.024 


K- 


0.783 ± 0.021 


0.668 ± 0.036 


P 


0.977 ± 0.027 


0.816 ±0.029 


P 


0.717 ± 0.038 


0.705 ± 0.067 



Table VII. - HERMES multiplicity ratios for various hadron species produced in e + Kr over e- 
deuterium collisions integrated over fractional hadron energies zh > 0.2 and 0.5 for struck quark energies 
V > 7 GeV [22]. Total experimental uncertainties are quoted. 



(see Sections 4 and 5). 

Recently, new data on neon and xenon and more krypton data have been collected at HER- 
MES [23]. This allows for a multidimensional analysis of the multipUcity ratio thus reducing the 
correlation between i.e. the v and Zh variables shown in the bottom panels of Fig. 17. In addition, 
the dependence of the of the other variables like Q^, p\ as well the mass number dependence 
can be studied. 

The pt dependence of the multiplicity ratio for identified charged hadrons in shown in Fig. 1 8 
for different nuclei. A nuclear enhancement is observed at p\ > 0.7 (GeV/c)^ similar to the 
Cronin effect observed in hadron-nucleus collisions (see Section 4'2). The conventional expla- 
nation of the Cronin effect \n p + A collisions ascribes this effect to the multiple scattering of 
projectile partons within the target nucleus (see Section 6'1). HERMES data highlight the role 
of partonic final-state multiple scattering (whereas Drell-Yan data are sensitive to initial-state 
parton multiple scattering, see Section 4'1) although explanations of the HERMES p^ broaden- 
ing in terms of prehadronic final-state interactions also exist, see Section 7'4. The larger Cronin 
enhancement for protons compared to mesons and antiprotons (note also the already discussed 
difference in the multiplicity ratio of p and p) is also seen mh + A and A + A collisions, where 
both baryons and antibaryons are enhanced compared to mesons, see Section 4. 

The two-dimensional analysis of the charged pion multiplicity ratio is presented in Fig. 19, 
where the multiplicity ratio is shown in three z/, -ranges as a function of v, and p^. The 
leftmost panels indicate that the dependence on v hardly depends on Zh- The 2^-dependence is 
similar for the various Zft-bins therefore, the dependence on Zh is not affected when integrating 
over Q^. The data in the rightmost plots indicate that the increase of for Kr and Xe at large 
p\ is smaller for larger Zh- Such a dependence was predicted in Ref. [29] and is consistent with 
the idea that the rise of R\ at large p], is of partonic origin. 

An analysis of the hadron p^-broadening has been recently performed at HERMES [88], 
which measured the observable A{p^}- {p\))^- iPj)'^, introduced in Section 23, for different 
hadron species and nuclear targets. The -broadening measurement is expected to provide 
new insights on the space-time evolution of the propagating quark [29], and on the multiparton 
correlation function inside the nucleus [89] (see Section 8' 1.3). 

The HERMES results for pj-broadening are shown in Fig. 20 for ;r* and produced on 
He, Ne, Kr, and Xe targets. The panels presented in Fig. 20 show (pj) for deuterium (top row) 
and the pr-broadening (remaining rows) as a function of v, Q^, x and Zh- The data do not reveal 
a significant dependence on v in the kinematic range covered. An increase of the broadening 
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Fig. 18. - HERMES multiplicity ratio versus for identified charged hadrons on krypton, neon and helium. 



with is observed, the behavior as function of x is very similar to the behaviour , due to the 
strong correlation between x and in the HERMES kinematics; hence it can not be excluded 
that the dependence is actually an underlying x dependence or both a and x dependence. 
The /7j-broadening is seen to vanish as z approaches unity while the {pj} for deuterium is 0.2 
or higher in the highest energy bin. The results on the mass number (A) dependence of the 
pr-broadening are shown in the left part of Fig. 21. The broadening is similar for tt'^ while is 
systematically higher for K^; it increases with A, however the uncertainties of the data do not 
allow to firmly establish the mass number dependence. In addition, the preliminary results for 
the HERMES ratio (p^)^ / {Pt)^ displayed in the left part of Fig. 22 as a function of Bjorken 
xb for charged hadrons, thus showing a non-negligible x-dependence of the broadening in the 
x-range explored. In particular the /:> 7- -broadening appears to increase from low to high xg, with 
a tendency to flatten out at xg > 0.2. 

As discussed in detail in Section 7, the new HERMES data allow one to definitely rule out 
nuclear absorption models based on one-step hadronisation process in which the struck quark 
propagates in the nucleus, interacts with the surrounding nucleons with perturbative cross section 
o-g, and fragments into a leading hadron in the vacuum. Instead, they support a more complex 
picture where the space-time evolution of the fragmentation process - as encoded e.g. in the Lund 
string fragmentation model - is largely modified by the surrounding matter. However, in spite of 
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the theoretical effort in developing new calculations for describing the wide range of data shown 
in this section, the observed (fj., z, p\) kinematical dependences can be described both in term of 
interaction of the intermediate prehadronic stage in absorption models (see Section 7), as well as 
in terms of parton energy loss calculations as discussed in Section 6. It is therefore important to 
extend the nDIS kinematical region with new measurements and to analyse other observables in 
order to disentangle the relative role of perturbative parton energy loss and prehadron absorption. 

3'4. Hadron production in e-nucleus DIS at CLAS/JLab. - Hadron production has been mea- 
sured with a 5.0 GeV electron beam with the CLAS detector at Jefferson Lab. Approximately 
25 fb"' of integrated luminosity was taken among the three primary solid targets, carbon, iron, 
and lead. In each case a 2-cm long liquid deuterium cryo-target was located in the beam simulta- 
neously, providing a normalisation for the nuclear ratios. A small amount of data was also taken 
on aluminum and isotopically enriched tin. The small spacing of 4-cm between the deuterium 
and solid targets limited the acceptance differences between the two targets for the large CLAS 
detector. 
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Particle identification for hadrons consisted of the standard CLAS instrumentation, which 
uses tracking and time-of-flight (TOF) systems, an electromagnetic shower calorimeter (EC) and 
a gas Cerenkov counter (CC). Identification of positive pions and protons was possible through 
the full momentum range using a combination of TOF and CC with the EC to reject positrons, 
while a somewhat more limited range of momentum was available for negative pions. Neutral 
pions and rj mesons were measured in their two-photon decay mode in the EC, while neutral 
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Fig. 21. ~ Mass number dependence of hadron pr-broadening in nuclear DIS. Left: HERMES data for n* 
and production with Ne, Kr and Xe targets. The inner error bars represent the statistical uncertainties, 
the total bars represent the statistical and systematic uncertainties. Right: CLAS preliminary data for n'^ for 
C, Fe and Pb targets threefold differential in Q^, v, and zi, as shown in the legend. The errors shown are 
statistical only. 



kaons were also measured in the tt'^tt^ channel, all over the full momentum range, but requiring 
a background subtraction. This very large data set is currently under analysis, and preliminary 
results are only available for positive pions and for neutral kaons at the present. The kinematical 
range of the data was for v = 2 - 4 GeV, - I - 4- GeV^; the full range of zi, and pj was 
available for analysis. Analysis cuts include y - v/v„,„^ < 0.85, > 4.0 GeV^, and target 
vertex cuts. While the CLAS/JLab data have a much more limited range in v, they offer two 
orders of magnitude more integrated luminosity than the HERMES data. This feature provides 
access to three-fold dimensional binning for at least positive pions, and should provide a first look 
at hadron formation in some previously unmeasured hadrons, such as rj, A, and K^. Because of 
the limited range in ^4-10 GeV^, the initial analysis is being focused on the region 
Zh = 0.4 - 0.7. All the (preliminary) data shown hereafter have not been corrected for acceptance 
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Fig. 22. - Preliminary x^-dependence of the A/D average (p^> ratio at HERMES (left) and CLAS (right). 
For the HERMES data the inner error bars represent the statistical uncertainties, the total bars represent the 
statistical and systematic uncertainties. For the CLAS data the errors shown are statistical only and the data 
have not been corrected for acceptance or radiative effects. 
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Fig. 23. - Left: CLAS preliminary data for the z = Zh dependence (left) and p\ = p\i^ of the hadronic 
multiplicity ratio for for carbon, iron, and lead targets. The data with solid (hollow) symbols correspond 
to the range v = 2.2 - 3.2 (3.2 - 3.7) GeV and 2^ = 1.0 - 1.3 GeV^. The errors shown are statistical only. 



or radiative effects; however, these coiTections are known to be rather small. 

Figure 23 shows preliminary data for the CLAS hadronic multiplicity ratio [90]. The data 
are shown in just two plots for compactness, but the statistical accuracy is adequate to divide 
the data further into multidimensional bins in Q^, v, and p^. In qualitative terms, these data are 
remarkably consistent with the HERMES results. The drop in multiplicity ratio with increasing Zh 
(and slow rise with increasing v) as well as the rise with increasing are observed in the JLab 
data as they are with the HERMES data (Figs. 17, 18, 19). The dependence on Q^, although 
visible, is also small in these data, as inferred from the HERMES studies. Thus, it is hoped that 
the two datasets can be inter-compared quantitatively in detailed model studies. 

In Fig. 2 1 right, are shown preliminary CLAS data for transverse momentum broadening (de- 
fined in Eq. (29)). Three sets of data are shown; each data point is binned in a multidimensional 
bin in Q^, v, Zn+ as well as A'^^. The available statistical sample is adequate to make more than 
two dozen sets of points with good statistical precision. The naive expectation that Ap^ is linear 
with the nuclear radius i.e. with A'^-*, is seen to approximately hold, although the data for the 
heaviest nucleus (lead) appear to undershoot the linear behavior. This flattening behaviour sug- 
gests the possibility that the partonic-level multiple scattering presumed to cause the broadening 
does not continue uniformly through the largest nucleus. If this picture is correct, then the pro- 
duction length can be estimated from the data in a rather direct fashion, using the well-known 
nuclear densities and sizes. 

In the right part of the Fig. 22 the preliminary jr^-dependence of the pj broadening at CLAS 
is shown for positive pions. Note that the xb range is different from HERMES. As can be seen, 
the Xb dependence is reasonably consistent for these two datasets, which are integrated over all 
other variables. The eiTors shown are statistical only; the CLAS data shown are not coiTected for 
small effects due to acceptance or radiative processes. 

3'5. Di-hadron correlations in nuclear DIS. - 

In nuclear DIS, double-hadron leptoproduction on A^, Kr and Xe relative to deuterium have 
been studied [24] via the ratio 

^A?-^.>0.5(^2)/*2\ 



(39) 



^zi>0.5 
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Fig. 24. - Double-ratio Rii, as a function of zi for ''*A' (squares), Kr (circles) and Xe (triangles) with 
Zi >0.5, compared to various theoretical predictions (solid for '"'A', dashed for Kr and dotted for Xe): 
BUU prehadronic transport model [91] (left), hadron absorption model [91] {centre), parton energy loss 
calculations [92] (right). 



where z\ and zi coiTespond to the leading (largest z) and sub-leading (second largest z) hadrons, 
respectively. The quantity c/A^''^"^ is the number of events with at least two detected hadrons 
in a bin of width dzi at zi with z\ > 0.5. The quantity N~'^^^'^ is the number of events with at 
least one detected hadron with zi > 0.5. The label A(D) indicates that the term is calculated for 
a nuclear (deuterium) target. 

If partonic energy loss of the struck quark were the only mechanism involved, it would be 
naively expected that the attenuation effect does not depend strongly on the number of hadrons 
involved, and the double-hadron to single-hadron ratio for a nuclear target should be only slightly 
dependent on the mass number A. On the contrary, if final hadron absorption is the dominant pro- 
cess, the requirement of an additional slower sub-leading hadron that is more strongly absorbed 
would suppress the two-hadron yield from heavier nuclei [29], so that this ratio should decrease 
with A. 

Results from HERMES are presented in Fig. 24, and show that the double-hadron ratio R2h 
is generally below unity with no significant difference between the '"^A^, Kr and Xe. The nuclear 
effect is much smaller than for the single-hadron attenuation measured under the same kinematic 
conditions. The displayed model computations, based on prehadron absorption [91] or parton 
energy loss [92,93], can reproduce the general trend in Z2 of the data, but predict a significant 
A-dependence which is not seen in the data (^). 

The small nuclear effect and its relative independence on A may simply point toward a strong 
surface bias of the photon-hadron interaction point. Such a bias is very natural in NLO dihadron 
production, in which the virtual photon's hard scattering produces two partons with energy vi_2 < 
V, which independently hadronise with fractional momentum zi_2 = 21,2^1,2/^- 

4. - Experimental results in hadron-nucleus collisions 
4'1. Drell-Yan production. - 

The production of dileptons with high invariant mass through the qq — > /^/" Drell-Yan (DY) 
process has been measured extensively in hadronic collisions, typically focusing on dilepton 
invariant masses between the charmonia and bottomonia masses (4 < M/+/- < 9 GeV/c^) or 



(^) Alternative string model approaches to two-particle correlations are discussed in [94, 95]. 
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measured by E866/NuSea [100]. SoUd curves are computations using the EKS98 shadowing parametrisa- 
tion [101]. Right: pr distribution in p+^^'^W normalised by p +' Be measured by E772/E866 [100, 102] 
compared to the calculations of [103]. 



above (Mz+,- > 1 1 GeV^) - see Refs. [96-98] for a review. NLO pQCD calculations describe 
well the experimental mass and transverse momentum distributions in hadronic collisions [99]. 
In hadron-nucleus collisions several nuclear effects are expected to modify the pQCD expectation 
of a linear dependence of the cross section with the mass number A (Eq. (2)). First of all, the 
modifications of parton densities in nuclei (nPDF) - such as shadowing at small parton fractional 
momentum X < 10"^, or the EMC effect at large x > 0.1 (see [20] for a recent review)- affect DY 
production in h+A with respect to that inh + p collisions. On top of nPDF effects, the projectile 
(anti)quark propagating through the nucleus may experience multiple scattering and lose some 
energy before the hard process takes place. In that sense, DY production data inh+A colUsions 
is particularly appropriate to study parton propagation in cold nuclear matter (see theoretical 
discussion in Section 63.3). 

Drell-Yan lepton pairs have been measured in p+A and n~ +A colUsions in fixed target experi- 
ments attheCERN SPS (NA3 [104], NAIO [105],NA38 [106]) and atFermilab (E772 [102,107], 
E866/NuSea [100] and older experiments [108]). The first measurements of DY production 
in nuclei were performed by the NA3 and NAIO collaborations using 150 GeV, 200 GeV and 
280 GeV pion beams on fixed-target nuclei ( yfsmi = 16.5-23 GeV). NA3 reported on the nu- 
clear dependence of DY production as a function of the longitudinal momentum-fraction of the 
incoming parton, x\ [104], while NAIO measured the pr-broadening of lepton pairs [105]. These 
data cover a phase-space range (0.2 < xi < 0.9) ideal to probe energy loss effects, even though 
the experimental uncertainties are quite large at the edge of phase space. Later, the Fermilab 
E772 experiment measured similarly DY production on various nuclear targets (Z), C, Ca, Fe, 
W) as a function of Xf using the 800 GeV proton beam (i.e. yfsNN - 38.7 GeV) [102]. More 
recently, the E866/NuSea collaboration reported on very high-statistics data in the same kine- 
matical region of E772 and using Be, Fe, and W nuclear targets [100]. 

Although the total cross section was found to approximately scale with A and therefore shows 
no nuclear effects within the experimental uncertainties, more differential measurements such 
as mean transverse momentum or longitudinal momentum (or xp) distributions actually exhibit 
a small, yet significant, nuclear dependence. Moreover, the data show /^^-broadening, also re- 
ported in hadron production in semi-inclusive lepton-nucleus DIS and hadron-nucleus collisions. 
Figure 25 shows the DY distribution ratios as a function of M, X2, xp and pr from the E866 [105] 
and E772 [102] experiments. At Xf > 0.2, more than 90% of the DY cross-section is due to the 
scattering of a quark from the hadron with an antiquark from the nucleus. Hence the DY process, 
unUke DIS, yields direct information on the nuclear modifications of the antiquark distribution. 
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Fig. 26. - Transverse momentum broadening {p])^ - (p^^)^ of Drell-Yan dileptons as a function of the 
invariant mass (NAIO dimuon data in ;r + A collisions [105], left) and as a function of A (E772 p + A 
data [96], right). 



The experimental results show no nuclear enhancement of the q distribution at moderate X2, con- 
trary to the predictions of "pion cloud" models which explain the EMC effect in terms of nuclear 
enhancement of exchanged mesons [96, 109-1 11]. Unless a large quark energy loss compensates 
for the predicted antiquark enhancement, DY data put strong constraints on nuclear models de- 
scribing the EMC effect [109-1 13]. At variance with the /?7-integrated ratios just discussed, the 
dilepton pr-broadening is mostly sensitive to the parton rescattering dynamics in the target nu- 
cleus. Measurements from NAIO as a function of M and from E772 as a function of A are shown 
in Fig. 26. While no significant dependence on the DY mass is observed by NAIO, the clear 
increase of the pj-broadening with the atomic mass number in the E772 data is qualitatively 
consistent with initial-state parton rescatterings. 

4'2. High-pT hadron production: "Cronin effect" . - 

High-energy proton-nucleus colhsions in fixed-target experiments at FNAL [8,9, 1 14], SPS [115] 
and at HERA-B [116] have observed an enhancement of the single inclusive hadron production 
yield relative to proton-proton collisions for transverse momenta above pr ~ 1 .5 GeV/c (Fig. 27). 
Such a high-pr enhancement, called "Cronin effect" [8,9], is evidenced by a nuclear modification 
factor that exceeds unity, and has also been observed in nDIS, see Figs. 16, 18 and 23 (right). 
The fact that RpA < 1 below pr ~ I GeV/c is simply because the incoherent binary scaling 
assumption behind Eq. (31), is not valid for soft particle production in hadronic collisions. 

The HERA-B collaboration has recently studied the production of mesons and A", A" 
baryons in /7+C,Ti,W interactions at y/s^N - 41.6 GeV [116]. The Cronin effect is clearly 
observed for all three species for transverse momenta above pj ~ 1.5 GeV/c (Fig. 28). The mass 
number dependence is parameterised as CTpA = CTpj^ ■ A" where cr^^v is the proton-nucleon cross 
section. The values of a are above one, in particular for the baryon species (A", A"). 

At RHIC, inclusive hadron production in minimum-bias deuteron-nucleus {d + Au) collisions 
at y^iviv = 200 GeV also features a small Cronin effect above pt = 1.5 GeV/c at mid-rapidity 
(Fig. 29). The enhancement is smaller, or practically absent, in the case of mesons (RdAH ~ 1) 
but is visible, of order R^au ~ 1 -4, for baryons. The enhancement peaks in the range pr ~ 2.5- 
4 GeV/c, and then starts to decrease and disappears beyond 8 GeV/c, an observation also apparent 
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Cronin enhancement in liigh-pr hadron production in proton-nucleus collisions at FNAL [8,9, 
compared to the predictions of [117], and at FNAL and SPS [115] (right). 



in the fixed-target data (Fig. 27). This peak structure is also progressively suppressed as the 
rapidity of the hadrons is increased [25, 118]. 

Such a transverse momentum broadening of high-/??- hadrons in hadron-nucleus collisions 
is usually interpreted (see Section 6'1) as due to either (i) multiple elastic scatterings of the 
incoming or outgoing parton inside the cold nuclear medium [117, 124], or (ii) recombination 
at the hadronisation stage of the scattered parton with other final state partons created in the 
collision [125, 126]. The larger Cronin effect for baryons than for mesons is naturally explained 
by the latter mechanism: the combination of three quark momenta boosts up the final baryon 
spectra more than the two-quark coalescence into mesons. The lower Cronin enhancement of 
all hadrons, in general, at RHIC centre-of-mass energies compared to fixed-target results can be 
explained by the steeper parton spectra at lower energies, which makes it easier to get a relatively 
larger boost for the same amount of kj "kick". The fast disappearance of the Cronin enhancement 
with increasing rapidities is well accounted for by models based on non-linear QCD evolution of 
the gluon densities in the nuclei [127-130]. The position of the maximum is thus connected to 




Fig. 28. - Cronin enhancement in high-pj- K° (left) and A (right) production in p+A collisions measured 
by HERA-B at ^/s^ = 41.6 GeV [1 16]. Note that o- « 1.05 corresponds to R^a = A" ' a 1.2 for Ti or W 
(with A = 22, 74). 
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Fig. 29. - Nuclear modification factor RdAu(PT) measured in d + An collisions at -yTsjviv = 200 GeV at 
RHIC. Mesons: 7t°, rj and oj (PHENIX [119, 120], top-left), and charged pions and kaons (STAR [121], 
top-right). Hadrons: tt^, and p + p (PHENIX [122, 123], bottom-left), and n^, h^, and S± (STAR [121], 
bottom-right). 



the rapidity-dependent value of the "saturation momentum" Qs(y) (see Section 61). 



5. - Experimental results in nucleus-nucleus collisions 
5'1. High-pr hadron production . - 

The production of hadrons at large transverse momentum has been since long proposed as 
a valuable "tomographic" probe of the hot and dense QCD matter produced in heavy-ion col- 
lisions [131]. If the final-state system is dense enough, the hard scattered partons will be at- 
tenuated while traversing it, resulting in a variety of "jet quenching" phenomena, e.g., suppres- 
sion of leading hadron spectra [132], distortion of azimuthal correlations between back-to-back 
jets [133, 134], modifications of the energy-particle flow within the final jets [135, 136]. The 
study of these modifications provides valuable information on the (thermo)dynamical properties 
of the produced system such as the initial gluon density dN^/dy, or the q transport coefficient 
characterising the "scattering power" of the medium. A detailed recent review of jet quenching 
results and phenomenology can be found in [137]. Here we highlight the main findings. 

The dominant contribution to the energy loss of partons is believed to be medium-induced 
gluon radiation as described in the Gyulassy-Levai-Vitev (GLV) [138, 139] and Baier, Dok- 
shitzer, Mueller, Peigne and Schiff (BDMPS) [140-142] (or LPCI [143]) formalisms (see Sec- 
tion 6'2). In the GLV approach, the initial gluon density dN^/dy of the expanding plasma (with 
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original transverse area Ax. = ttR^ x 150 fm^ and thickness L) can be estimated from the mea- 
sured energy loss AE: 

(40) AE oc ale R L, 

Aj. dy 

where Cr is the Casimir colour factor of the parton (4/3 for quarks, 3 for gluons). In the BDMPS 
framework, the transport coefficient q - characterising the squared average momentum transfer 
of the hard parton per unit path-length: q = m^/A - rn^^pcr, where mo is the medium Debye 
mass, p its density, and cr the parton-matter interaction cross section - can be derived from the 
average energy loss via: 

(41) {AE)^asCR{q}L\ 

For example, for an equilibrated gluon plasma at T = 0.4 GeV with coupling a,, = 0.5 - i.e. with 
density pg = 16/7T- ^(3) ■ ^ 15 fm'^, Debye mass mp = {47Ta,y^^T ^ 1 GeV/c^ and LO 
perturbative cross section o-gg * 1.5 mb, one finds q ^ 2.2 GeV^/fm [144]. 

Experimentally, the standard method to quantify the medium effects on the yield of a large- 
Pt particle produced at rapidity y inaA + A reaction is given by the nuclear modification factor 
RAA(PT,y', b), see Eq. (31), which measures the deviation of hadron spectra in A -H A collisions 
at impact parameter b from an incoherent superposition of spectra in nucleon-nucleon collisions 
(Raa = !)■ If the A + A and p + p invariant spectra are both a power-law with exponent «, i.e. 
1/Pt dN/dpT Pj", the fraction of energy lost eioss - ^Pt/Pt, can be (grossly) estimated from 
Raa eioss ^ I - R]^''^^ [145]. 

Among the RHIC highlights is the observation (Fig. 30, left) of a strong suppression of 
the products of parton fragmentation such as high-pj- neutral mesons (n", rf) [146, 152, 153] 
and charged hadrons [149, 154-156]. The Raa ~ 1 perturbative expectation which holds above 
Pt ~ 4 GeV/c for other hard probes such as "colour blind" prompt photons [157, 158] (Fig. 30, 
left) or for mesons in d + Au reactions (Fig. 29, top), is badly broken in central Au + Au col- 
lisions where one measures Raa ~ 0.2. The measured tt" power-law spectral exponents of 
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n » 8 and the Raa ~ 0.2 imply an average fractional energy loss of the high-pj- hadrons, of 
eioss ~ ^Pt/Pt ~ 0.2 [145]. Such a significant suppression was not observed at SPS where 
- even after reevaluating the p + p baseline spectrum [159, 160] - the central Pb + Pb me- 
son spectra show an Raa around unity (Fig. 32). Yet, this does not exclude the possibihty of 
energy loss at SPS, since the factor of the ~50% Cronin enhancement observed in p + Pb at 
yfsNN ~ 20 GeV [161] likely compensates for the same amount of final-state suppression in the 
hot medium [159, 162] or in the cold nuclei [43]. 

As discussed next, most of the empirical properties of the suppression factor are in quantita- 
tive agreement with the results of the non-AbeUan parton energy loss models presented in more 
detail in Section 6 2. 

5' 1 . 1 . Magnitude of the suppression and medium properties. - The Au + Au high-pr suppres- 
sion can be well reproduced by parton energy loss calculations in a very dense system with initial 
gluon rapidity densities dN^/dy a; 1400 (GLV curve in Fig. 30, left) [139], plasma temperatures 
T « 0.4 GeV (AMY model) [163], and time-averaged transport coefficients (q) ^ 13 GeV^/fm 
(PQM model in Fig. 3 1 , left) [ 1 20, 1 64] , or initial-time transport coefficients §o ~ 1 - 1 8 GeV ^fm 
(ASW curve in Fig. 31, right). 

The consistency between the extracted q, dN^/dy and T values in the various models has 
been studied e.g. in [137, 165, 166]. Whereas the agreement between the fitted thermodynamical 
variables dN^ jdy and T is good, the values of the transport parameter q extracted from the data 
within the various models differ by factors of 2 - 3, and are much larger than the LO BDMPS 
estimate § « 2 GeV^/fm at T = 0.4 GeV given before(^). The discrepancy with the perturbative 
estimate is a factor K - 3.6 in the case of ASW [166] and an order of magnitude regarding 
the value extracted within PQM. At least part of the uncertainty is due to the relative insensitiv- 
ity of the q parameter to the irreducible presence of hadron from (unquenched) partons emitted 
from the surface of the plasma [167]. Although relating transport properties to thermodynamical 
quantities is model- (and medium-) dependent, an accord between the fitted q and dN^/dy val- 
ues - via q = mjj per with p oc dN^/dy - can only be seemingly achieved for parton-medium 
cross-sections much larger than the o-gg - €?(1 mb) LO perturbative expectation. Such an ob- 
servation lends support to the strongly-coupled nature of the QGP produced at RHIC [168]. It is 
however not satisfactory to obtain parameter values which are typically non-perturbative based 
on a purely perturbative framework, see e.g. the discussion in [144]. Additional constraints on 
q can be placed by requiring also the model reproduction of the suppressed dihadron azimuthal 
correlations (see Section 5'3). 

51.2. Universality of (light) hadron suppression. - Above pr ~ 5 GeV/c, ti^ [153], t] [147], 
and inclusive charged hadrons [149, 156] (dominated by n- [156]) show all a common factor of 
~5 suppression relative to the Raa = 1 perturbative expectation which holds for hard probes, such 
as prompt photons, insensitive to final-state interactions [157] (Fig. 30, left). Such a "universal" 
hadron deficit is consistent with in-medium partonic energy loss of the parent quark or gluon prior 
to its fragmentation in the vacuum. The high-pj- n^/jj ratio is indeed found to be independent 
of the collision system (or centrality in Au + Au colUsions) within uncertainties (Fig. 30, right). 
Since both hadrons have similar valence quark content but very different masses (the is four 
times heavier than the n^), and different meson-nucleus cross sections, a natural explanation 
for the observation is that quenching happens at the parton level: the medium induced energy 



It is also interesting to note that the BDMPS transport coefficient for cold nuclear matter « 
0.6 GeV/fm^, extracted via Eq. (58) from hadron quenching in nDIS is also a factor 10 larger than the 
perturbative estimate ~ 0.05 GeV/fm? discussed in Eq. (54). 
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Fig. 31. - Left: Raa(Pt) for pions in central An + Au collisions (triangles) [120] compared to PQM predic- 
tions [164] for varying values of the q coefficient (red curve, best fit for <§> = 13 GeV^fm). Right: Central 
and semi-central pion suppression at PHENIX [120] compared to AMY, HT and ASW models within a 
common space-time evolution medium profile [166]. 



loss experienced by the parent parton is the same independently of the meson which is to be 
produced, and hadronisation occurs in vacuum, according to the same fragmentation functions 
extracted from e^e" and p + p collisions. 

5' 1.3. Centre-of-mass energy dependence. - As one increases the collision energy in A -i- A 
collisions, the produced medium reaches higher energy and particle densities, the system stays 
longer in the QGP phase, and correspondingly the traversing partons are more quenched. Fig- 
ure 32 compiles the measured Rj\j^(^pj) for high-/?^ tt" in central A + A collisions at sfst^ ^17.3 
and 200 GeV compared to parton energy loss calculations that assume the formation of a QGP 
with initial gluon densities dN^ /dy ~ 400, 1400 [139, 169] or, equivalently, averaged transport 
coefficients (q) * 3.5, 13 GeW^/fm [164] respectively. Table VIII collects these results as well as 
those for central Am +Au reactions at y/s^N ~ 62 and 130 GeV [152, 170]. For each collision en- 
ergy the derived values for dN^/dy are consistent with the final charged hadron density dNch/drj 
measured in the reactions(^). This is expected in an isentropic expansion process, where all the 
hadrons produced at midrapidity come directly from the original gluons released in the collision: 



dm N,o, 

(42) K — 

ay Nch 



dNch , „ dNch 



drj drj 



This relationship is relatively well fulfilled by the data as can be seen by comparing the fourth 
and fifth columns of Table VIII. 

5' 1 .4. Transverse momentum dependence. - At RHIC top energies, the hadron quenching fac- 
tor remains basically constant from 5 GeV/c up to the highest transverse momenta measured so 
far, pj « 20 GeV/c (Fig. 30). As can be seen in Eq. (45), the suppression is roughly propor- 
tional to the z-slope of the fragmentation functions, dD(z, Q^)ldz. Consequently, even though 



(^) The charged particle multiplicity itself follows a logarithmic dependence on the cm. energy [171]: 
dN^iJdr] » 0.75-(A'parr/2)-ln( -s/snn [GeV]/1.5) (Npun is the number of participant nucleons in the collision). 
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Fig. 32. - Nuclear modification factor, Raa{pi), for liigh-px :^"'s at CERN-SPS [159, 160] and RHIC [120] 
compared to the predictions of the GLV parton energy loss model (dN^/dy = 400, 1400) [139]. The bottom 
bands show the predicted R^a for high-px charged hadrons in central Pb + Pb collisions at ^snn =5.5 TeV 
as given by the GLV (dNg/dy = 2000-4000) and the PQM (q « 30-80 GeV^fm) models [172]. 



the relative parton energy loss, elkj, becomes smaller at higher pr (leading naively to a smaller 
suppression and an increase of Raa with pr), the larger steepness of the partonic spectrum due to 
the restricted phase space to produce hight-p^ partons at RHIC leads to a significant suppression 
even at large transverse momenta. Indeed, full calculations [139, 164, 167, 174] including the 
combined effect of (i) energy loss kinematics constraints, (ii) steeply falling pr spectrum of the 
scattered partons, and (iii) ff(2Q%) p^-dependent (anti)shadowing differences between the pro- 
ton and nuclear parton distribution functions (PDFs), result in an effectively flat Raa{Pt) as found 
in the data. The much larger kinematical range opened at the LHC TeV-energies [172] will allow 
one to test the pj-dependence of parton energy loss, and the associated radiation spectrum, over 
a much wider domain than at RHIC. As can be seen in Fig. 32 (yellow bands) the PQM model 
seemingly predicts a slower (and smaller) rise of Raa{pt) than the GLV model. 

5'1.5. Centmlity (system-size) dependence. - The volume of the produced plasma in a heavy- 
ion collision can be "dialed" modifying the overlap area between the colliding nuclei either by 
selecting a given impact-parameter b - i.e. by choosing more central or peripheral reactions - or 
by colliding larger or smaller nuclei, e.g. Au {A - 197) versus Cu {A - 63). The relative energy 



Table VIII. - Suppression factors measured in central A + A collisions in the range ^[snn ~ 20 - 200 GeV, 
and initial gluon densities dN^ /dy [139, 169], and transport coefficients (q) [164, 173] obtained from parton 
energy loss calculations reproducing the observed high-pj tt" suppression at each ^ s^n- The measured 
charged particle densities at midrapidity, dN"^'' / dr]\,po [171], are also quoted. 



W(GeV) RAA(n",pT-^ 4 GsY/c) <§)(GeVVfm) d^/dy dN'J /dri\„=o 



SPS 17.3 -1.0 [159, 160] 3.5 400 312 ±21 

RHIC 62.4 -0.4 [170] 7. 800 475 ± 33 

RHIC 130. -0.3 [152] 11 1000 602 ± 28 

RHIC 200. -0.2 [153] 13 1400 687 + 37 
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Fig. 33. - Left: Centrality (N,,„,i) dependence of the high-pr n° suppression in Cu + Cu and Au + Au at 
200 GeV [175]. Right: Raa(Pt) for in central Cu + Cu collisions at ^fs^ =22.4, 62.4 and 200 GeV 
compared to GLV calculations with initial gluon densities dN^ /dy a 100 - 370 [176]. 



loss depends on the effective mass number Agfj or, equivalently, on the number of participant 

2/3 2/3 

nucleons in the collision A^part, as: A/i/Zi oc A^^. oc A^p^^^ [164, 177]. The measured T^AACpr) 
in central Cu + Cu at 22.4, 62.4, and 200 GeV [176] is a factor of (AaJAq,)^^^ ~ 2 lower than 
in central Au + Au at the same energies (Fig. 33, right) Yet, for a comparable Npan value, the 
suppression in Au + Au and Cu + Cu is very similar (Fig. 33, left). Fitting the Npan dependence 
to R/^/^ - {\ - K Npg„)"^^ yields a - 0.56 ±0.10 [120], consistent also with parton energy loss 
calculations [164, 177]. 

5' 1.6. Path-length dependence. - The analytical quadratic dependence of the energy loss on 
the thickness of a static medium L, Eq. (41), becomes effectively a linear dependence on the 
initial value of L when one takes into account the expansion of the plasma, see Eq. (40). Ex- 
perimentally, one can test the dependence of parton suppression on the plasma thickness (L) by 
exploiting the spatial asymmetry of the system produced in non-central nuclear collisions. Par- 
tons produced "in plane" ("out-of-plane") i.e. along the short (long) direction of the ellipsoid 
matter with eccentricity e will comparatively traverse a shorter (longer) thickness. PHENIX has 
measured the high-pj neutral pion suppression as a function of the angle with respect to the 
reaction plane, Raa(Pt, <f>) [145, 178]. Each azimuthal angle (p can be associated with an average 
medium path-length via a Glauber model. The energy loss is found to satisfy the expected 
AE oc L dependence, Eq. (40), above a "threshold" length of L ^ 2 fm, interpreted in [179] as 
due to a geometric "corona" effect. 

5' 1.7. Non-Abelian (colour factor) dependence. - The amount of energy lost by a parton 
in a medium is proportional to its colour Casimir factor: Ca = 3 for gluons, Cf - 4-/3 for 
quarks. Asymptotically, the probability for a gluon to radiate another gluon is Ca/Cf = 9/4 
times larger than for a quark and thus g-jets are expected to be more quenched than q-jets in 
a QGP One can test such a genuine non-Abelian property of QCD energy loss by measuring 
hadron suppression at a fixed pr for increasing cm. energy [173, 180]. At large (small) x, the 
PDFs are dominated by valence-quarks (low-.x: gluons) and consequently hadroproduction will 
be dominated by quark (gluon) fragmentation. Figure 34 (left) shows the Raa for 4-GeV/c pions 
measured at SPS and RHIC compared to two parton energy loss curves [180]. The lower (upper) 
curve shows the expected Raa assuming a normal (arbitrary) behaviour with AEg/AEq - 9/4 
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Fig. 34. - Left: Excitation function of the nuclear modification factor, Raa('^Snn), for n° in central AA 
reactions at a fixed pr = 4 GeV/c [173], compared to non-Abelian (solid) and "non-QCD" (dotted) energy 
loss curves [180]. Right: Raa(Pt) for pions and (anti)protons in central Au + Au at ^fs^N = 200 GeV [181]. 



{^Eg — AEq). The experimental high-p^ tt data supports the expected colour-factor dependence 
of RAAiyf^ [173]. 

A second test of the colour charge dependence of hadron suppression has been proposed 
based on the fact that gluons fragment comparatively more into (anti)protons than quarks do. 
One would thus naively expect R'^ < R'^. The STAR results (Fig. 34, right) are however 
seemingly at variance with this expectation: pions appear more suppressed than protons at high- 
er [181]. Yet, the use of (anti)protons as a reference for perturbative particle production is ques- 
tionable. First, it is worth reminding that the assumption of production from parton fragmentation 
in vacuum may well not hold for protons which have estimated formation time-scales a factor 
~5 shorter than for pions (Table VI). Second, p,p are already found to be Cronin-enhanced in 
d + Au compared to p + p collisions by a factor ~ 50% - 100% (see Fig. 29 bottom-left) for pj's 
as large as 7 GeV/c [182]. It is likely that there is an extra mechanism of baryon production, 
based e.g. on in-medium quark coalescence [183-185] (see Section 6'1), which compensates for 
the energy loss suffered by the parent partons. Finally, another explanation has recently been 
proposed in Ref. [186]: protons could be produced in a compact colour-singlet configuration 
(through higher-twist processes) and therefore, because of colour transparency, escape more eas- 
ily the dense medium than pions do. 

5'2. Heavy flavour production . - 

As seen in the previous Section, most of the empirical properties of the quenching factor 
for light-flavour hadrons - magnitude, pr-, centrality-, yfsNN- dependences of the suppression - 
are in quantitative agreement with the predictions of non-Abelian parton energy loss models. A 
robust prediction of radiative energy loss models is the hierarchy A/Sheavy-Q < A/siight-q < AEgi 
gluons, which fragment predominantly into light hadrons, are expected to lose more energy than 
quarks because of the larger color coupling to the radiated gluon; whereas massive c, b quarks 
are expected to lose less energy than light quarks due to a suppression of small-angle gluon 
radiation already in the vacuum ("dead-cone" effect) [187, 188]. Grossly, the in-medium radiative 
energy loss is suppressed by a factor ffimalMo) (where mo ~ 1 GeV/c^ is the medium Debye 
mass) [189], i.e. it is a factor about 25% (75%) less for a charm (bottom) quark than for a 
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Fig. 35. - Nuclear modification factor R'^ of heavy-quark decay electrons in central Au + Au collisions at 
yjjvjv = 200 GeV [190, 191] compared to various radiative+elastic energy loss models for c and b quarks 
(left) and to a model of heavy-meson (D, B) dissociation in a plasma (right) [30]. 



light-quark. 

Yet, RHIC measurements [190-192] of high-pj- electrons from the semi-leptonic decays of 
D- and B-mesons (Fig. 35) indicate the same suppression factor for light and heavy mesons: 
l^AAiQ) ~ RAAiq, g) ~ 0.2. In order to reproduce the high-pr open charm/bottom suppression, jet 
quenching models require either initial gluon rapidity densities dN^ jdy ~ 3000 [193] which are 
inconsistent with the total hadron multiplicities, see Eq. (42), as well as with the dN^/dy ~ 1400 
needed to describe the quenched light hadron spectra. Various explanations have been proposed 
to solve this 'heavy flavor puzzle' (see e.g. [194, 195]). 

• First, precise comparisons between theory and data require a better determination of the 
relative contribution of the c and b quarks to the measured non-photonic electron yields [193, 
196]. If only c quarks (roughly three times more suppressed than the heavier b quarks) 
actually contributed to the measured high-pj- decay electron spectrum, then one would 
indeed expect /?^y\(c) ^ 0.2 [197]. However, indirect measurements from PHENIX [198] 
and STAR [199] have confirmed the similar production yields of from D and B decays 
above pr ~ 5 GeV/c predicted by next-to-leading-log pQCD [196,200]; 

• The second mechanism points to an additional contribution from elastic (i.e. non-radiative) 
energy loss for heavy-quarks [201,202] which was considered negligible so far [132]. 
Recent works [203-206] have shown that a proper evaluation of the QCD running coupling 
substantially increases the amount of collisional energy loss suffered by the heavy-quarks. 
AEcoii can indeed be a significant contribution for heavy quarks (see 'rad.-i-el.' curves in 
Fig. 35, left); 

• Two works [207, 208] have argued that in a plasma the large charm-quark coalescence into 
Ac baryons (with a small semileptonic decay branching ratio) would deplete the number of 
open-charm mesons, and coiTespondingly reduce the number of decay electrons, compared 
to p + p collisions; 

• The assumption of vacuum hadronisation (after in-medium radiation) implicit in all parton 
energy loss formalisms may well not hold in the case of a heavy quark. The discussed 
quark-hadronisation time estimates (see Section 2'4) are inversely proportional to the mass 
mi, of the final produced hadron: the heavier the hadron, the fastest it is formed. From 
Eq. (35), one can see that in the rest frame of the hot QCD medium of the fragmenting 
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Fig. 36. - Comparison of the azimuthal di-hadron correlation dNp„ii /dAipdrj for p + p (open symbols) and 
central Au + An (histograms) at yfs^ = 200 GeV for p'^''^ = 5-10 GeV/c and increasingly smaller (right to 
left) values of [211]. 



heavy-Q, the formation time of D- and B-mesons [30] is of order(^) (f/,) x 0.4 - 1 fm, 
respectively. Thus, one may need to account for both the energy loss of the heavy-quark as 
well as the possible dissociation of the heavy-quark meson inside the QGP. The expected 
amount of suppression in that case is larger and consistent with the data (Fig. 35, right). 



5'3. High-pT di-hadron correlations . - 

At leading order, the parton-parton 2 — > 2 scatterings are balanced in pr i e. they are back-to- 
back in azimuthal angle, Aip = n, modulo some smearing due to the partons intrinsic transverse 
momentum. Such azimuthal correlation is smeared out if one or both partons suffer rescatter- 
ings inside the plasma. The dijet-acoplanarity arising from the interactions of a parton in an 
expanding QGP is {k^)med - {in^/^)L ln(L/To) oc qL [209] and, thus, the final azimuthal cor- 
relations between the hadrons issuing from the fragmentation of quenched partons will show a 
dependence on the transport coefficient and thickness of the medium: d-Npairldh.(p - f{q,L). 
Jet-like correlations in heavy-ion collisions have been measured on a statistical basis by select- 
ing a high-pj- trigger particle and measuring the azimuthal (A^ = <f> - <Prng) and pseudorapidity 
(At; = T]—T]frig) distributions of its associated hadrons (p^'™' < p'"^): C(Acf>,Ar]) - jf—'j^j^- In 
p+p collisions, a dijet signal appears as two distinct back-to-back Gaussian-like peaks at A0 x 0, 
Arj ^ (near-side) and at A(p ^ n (away-side). At variance with such a topology, early STAR 
results [210] showed a monojet-like topology with a complete disappearance of the opposite-side 
peak for p"/"" ^2-4 GeV/c. 

Fig. 36 shows the increasingly distorted back-to-back azimuthal correlations in high-/?j- trig- 
gered central Au + Au events as one decreases the pr of the associated hadrons (right to left). 
Whereas, compared to p + p, the near-side peak remains unchanged for all pr's, the away-side 
peak is only present for the highest partner pj's but progressively disappears for less energetic 
partners [211, 212]. The correlation strength over an azimuthal range A(p between a trigger hadron 
h, and a partner hadron in the opposite azimuthal direction can be constructed as a function of 



(*) Note that in the laboratory system there is an extra Lorentz boost factor: t^^,^ = 7e ' ''"j 
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Fig. 37. - Le/f.' D''^(zt) distributions for d + Au and Au + Au and Iaa(zt) ratio for central An + Au at 
200 GeV [214], compared to HT calculations for varying eg energy losses [213]. Right: Data vs. theory 
values for the fitted eo parameter [213]. 



the momentum fraction zt - p^^""^ / p'^^'^ via a "pseudo-fragmentation function": 

r ,n<. r r d^o-''/!'/dp':'>^dp"'""dA(f> 

(43) - (zt) ^ \ dpj. dpj dAcf> . 

J J JA(A>i3o° dcr^^/dpj'^ 

shown in Fig. 37 (top-left) compared to various values of the eo parameter characterising the 
amount of energy loss [213] (see the Higher- Twist model discussion in Section 6'2.1). Simi- 
larly to Raa{Pt), the magnitude of the suppression of back-to-back jet-like two-particle corre- 
lations can be quantified with the ratio Iaa{zt) — Daa(zt)/Dpp{zt)- I"^' is found to decrease 
with increasing centrality, down to about 0.2 - 0.3 for the most central events (Fig. 37, bottom- 
left) [210,214]. 

The right plot of Fig. 37 shows the best eo a; 1.9 GeV/fm energy loss value that fits the 
measured Raa and Iaa factors. Due to the irreducible presence of (unquenched) partons emitted 
from the surface of the plasma, the single-hadron quenching factor Raa(Pt) is in general less 
sensitive to the value of eo than the dihadron modification ratio Iaa^Zt)- The combination of 
Raa{Pt) and Iaa{zt) provides more robust quantitative information on the medium properties. 

Since energy and momentum are conserved, the "missing" fragments of the away-side (quenched) 
parton at intermediate p^'s should be either shifted to lower energy {pj < 2 GeV/c) and/or scat- 
tered into a broadened angular distribution. Both, softening and broadening, are seen in the data 
when the pr of the away-side associated hadrons is lowered (see two leftmost panels of Fig. 36). 
As a matter of fact, the away-side hemisphere shows a very unconventional double-hump angular 
distribution with a "dip" at A(f) » n and two neighbouring local maxima at » tt ± 1.1 - 1.3. 
Such a "volcano"-like profile has been interpreted as due to the preferential emission of energy 
from the quenched parton at a finite angle with respect to the jet axis. This could happen in 
a purely radiative energy loss scenario due to large-angle radiation [215], but more intriguing 
explanations have been put forward based on the dissipation of the lost energy into a collective 
mode of the medium in the form of a wake of lower energy gluons [216] with Mach- [217, 218] 
or Cerenkov-like [219,220] angular emissions assuming that the hard parton velocity exceeds 
the speed of sound or speed of light, respectively, of the QCD medium. Theoretically, it is 
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Fig. 38. - Prompt photon spectrum compared to preliminary PHENIX data. Left: Spectrum computed in 
Refs. [222-224], Plot taken from Ref. [225]. Right: Spectrum computed in the GLV formalism including a 
variety of effects. Plots taken from Ref. [226]. 



unclear if such partonic collective wake(s) and cone survive both hadronisation and the final 
hadronic freeze-out [221]. More involved studies, e.g., accounting for the plasma expansion and 
the hadronic phase evolution, are needed before a final conclusion can be reached. 

5'4. High-px photon production . - 

Photon production in cold and hot QCD media is a particularly interesting probe because, 
unlike the gluon, photons can escape without final-state interactions, therefore caiTying infor- 
mation about the medium at the location of its production. Moreover, it avoids the ambiguities 
related to the hadronisation process, so that it is in principle a cleaner probe of the underlying 
dynamics. However, several mechanisms proposed over the last few years showed that prompt 
photon production could actually be somehow affected by the quark-gluon plasma formation in 
heavy-ion collisions, and in principle by cold nuclear matter as well. 

First, prompt photons are also produced from the collinear fragmentation of quarks and glu- 
ons produced in the hard process. Such fragmentation photons are thus sensitive to energy 
loss of their parent fragmenting parton. The formation time needed to produce such a photon- 
parton system with a small invariant mass exceeds by far the typical lifetime of the hot medium: 
the multiple scattering of the hard parton in the medium is followed on a much larger time- 
scale by the parton-to-photon fragmentation process in the vacuum. Consequently, one expects 
that those "fragmentation" photons should be as quenched, at least qualitatively, as the hadron 
yield [222,227,228]. 

Second, the multiple scattering incuiTed by the hard partons traversing the produced medium 
induces the emission of soft gluons - leading to the usual "jet quenching" - as well as soft 
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photons [163, 224, 229] (although with a probabihty in principle governed by the much smaller 
ttem)- Within a path-integral picture for the parton energy loss mechanism [143, 230], Zakharov 
computed this medium-induced photon bremsstrahlung contribution at RHIC [231]. The en- 
hancement of photon production is particularly noticeable in the moderate p^^ range, say below 
20 GeV, raising hopes that it could be measured. It has also been proposed that large pt par- 
tons may couple to the thermal quarks and gluons in the medium through Compton scattering 
or qq annihilation. This so-called jet-photon conversion mechanism, [232], would enhance the 
production of photons, leading to a nuclear production ratio, /?aa, larger than one. Each one of 
these individual mechanisms, with opposing (quenching and enhancement) effects, is unfortu- 
nately under poor theoretical control and the medium effects on (prompt) photon production still 
remain rather model-dependent. 

Experimentally, the photon quenching factor has been measured at RHIC by the PHENIX 
collaboration, up to large pr ~ 20 GeV/c [225]. Despite the stiU rather large error bars of 
those preliminary data, it seems that the prompt photon /?aa is consistent with unity in the range 
p-j- » 4-15 GeV/c, while a suppression is reported, /?aa - 0.6, in the highest pt bin. Therefore it 
appears that the photon enhancement predicted either due to parton multiple scattering [224] or 
jet-photon conversion [232] is not seen at RHIC. The BDMPS energy loss calculation, Ref. [228], 
supplemented by the proper treatment of protons and neutrons in the Au nuclei (leading to a 
suppression at large pt due to smaller electric charge of d compared to u valence quarks), is on 
the contrary able to reproduce the shape and magnitude of the PHENIX data. 



6. - Parton propagation and energy loss 

6'1. In-medium parton propagation and the Cronin effect . 

'['lie Cronin effect is the enhancement of single inclusive hadron production at intermediate 
Pr'sipr ~ 1 - 8 GeV/c) observedinhadroproductionin e-H A (Fig. 18) and/n-A(Fig. 27) colli- 
sions, as well as in Drell-Yan events (Fig. 25, right). Theoretically, such an effect has been usually 
explained in terms of initial- or final-state multiple scattering of the parton prior to its fragmen- 
tation leading to a broadening of the transverse momentum of the produced hadrons [124]. More 
recently, a modification of the hadronisation mechanism due to the recombination of the scattered 
parton with other partons produced in the collision has been proposed [125, 126] and accounts 
for the experimental data as well [233]. We review here both theoretical interpretations of the 
Cronin enhancement. 

6' 1.1. Parton multiple scatterings. - Parton multiple scatterings have been discussed in the 
pQCD factorisation formalism, in the colour dipole model and in the Colour Glass Condensate 
approach. We briefly review their main features, similarities and differences. 



• Glauber-Eikonal models. 

The Glauber-Eikonal (GE) approach [49,234-237] to the Cronin effect treats multiple 2 — > 
2 partonic collisions in coUinearly factorised pQCD. The cross-section for the production 
of a hadron with transverse momentum px and rapidity yinh+A colUsions at fixed impact 
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parameter b is written as [49] 
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where the crossed-circle symbols denote appropriate integrations and summations over 
parton flavours / and j, (pnh are the PDFs in the projectile hadron, dcriN is the cross section 
for a parton ; scattering on a nucleon and TA(b) the A nucleus thickness function at impact 
parameter b (see Ref. [49] for details). At moderate pr, the accumulation of transverse 
momentum due to multiple scatterings leads to an enhancement of transverse momentum 
spectra, and to a suppression of the low-pr region due to energy-momentum conservation. 
At high pt the binary scaled p + p spectrum is recovered: no high-/?r suppression is pre- 
dicted in this approach, except as a consequence of the nuclear modifications of PDFs. In 
early appUcations, the GE series (44) has been directly evaluated only up to the n - 3 par- 
ton scatterings, under severe approximations, and only for ^/s < 40 GeV [114, 234, 235]. 
Instead of evaluating the full GE series, other approximated GE models modify the pQCD 
rates through the inclusion of a phenomenological nuclear broadening of the intrinsic par- 
ton momentum kj [139, 238-240], ignoring however the unitarity constraints built into the 
GE multi-scatterings. The GE series is directly computed via a numerical convolution of 
elementary parton-nucleon processes, assuming a decoupling of the transverse and longi- 
tudinal kinematics in Refs. [49, 162, 241], and with exact energy-momentum conservation 
up to 3 scatterings in Ref. [242]. 

The formulation preserves unitarity and is directly constrained to reproduce the abso- 
lute normalised spectra in p + p collisions, and quantitatively incorporates kinematic 
phase-space limitations at large pr- Thus, the /57-broadening is computed without ad- 
justable parameters, rather than assumed as an input. Midrapidity pion production data 
at '\/snn = 20 - 200 GeV can be well described (Fig. 39 left), but the large hadron sup- 
pression reported at forward rapidity at RHIC [243, 244] (Fig. 39 right) can be reproduced 
only within gluon saturation approaches [245, 246] or using extreme nuclear shadowing 
parametrisations [247] . 

• The colour dipole model 

In Ref. [117,249], the GE series is formulated in terms of the multiple scattering of a 
colour dipole on the nuclear target. The colour dipole cross section is determined phe- 
nomenologically by fits to lepton-proton and proton-proton scattering data [250,251], and 
hadron production is determined from the overlap of the nuclear broadened dipole wave 
function with the hadron light-cone wave function. The computation of nuclear efl'ects 
is carried out with no tunable parameters, and leads to a good description of midrapidity 
pion data (Fig. 27 left). This model allows for the inclusion of coherent multiple scatter- 
ing [117], relevant at RHIC to LHC energies, and is equivalent to the GE model in the 
incoherent scattering limit. It can also describe in a unified formalism the /^j- -broadening 
of DY lepton pairs [103, 252] and that of hadrons produced nDIS [29] (see Section 7'3). 
This approach also allows for a description of the forward rapidity hadron suppression at 
RHIC due to energy conservation and Sudakov suppression at the edge of the phase space, 
alternative to gluon saturation [253, 254]. Hadron suppression is predicted at forward ra- 
pidity also for lower energy collisions, -\/%7v = 63 or 130 GeV, where gluon saturation 
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Fig. 39. - Left: GE multi-scattering computation [124] of neutral pion RdAuiPr) at midrapidity in different 
centrality bins compared to PHENIX data at ^/s^ = 200 [122]. Plot taken from Ref. [122]. Right: 
BRAHMS data on negatively charged hadron RdAuipr) at 77 = 3.2 compared to several models. From top 
to bottom: pQCD computation supplemented by nuclear shadowing effects by Guzey etal. [50]; GE model 
computation by Accardi et al. [124]; CGC model computations by Tuchin et al. [245] and by Jalilian- 
Marian [246]. Plot from Ref. [248]. 



effects are not expected to play a major role because of the rather large X2 probed at these 
energies. 

• Colour Glass Condensate 

The Colour Glass Condensate (CGC) is an effective theory for the gluon field in hadron 
and nuclei at small x, see Ref. [255] for a recent review. At low-jc the gluon occupation 
number is so large that the gluon field can be treated semi-classically and computed as a 
solution of the Yang-Mills equation of motion in the presence of random colour sources. 
Gluons with momenta lower than a scale Qs are in the "saturation" regime with a den- 
sity high enough that 2 — > 1 gluon-gluon fusion processes limit a further growth; on 
the contrary at larger momenta the gluon field is in the DGLAP "dilute" regime. An in- 
termediate "geometric scaling window" extends at momenta Qs < pr < Q^/Qo, where 
quantum effects from the saturation region further modify the evolution from the pertur- 
bative to the saturation region. Observables are computed as an average over the colour 
sources density p with a weight Wv[p], depending on the gluon rapidity y = log(l /x). The 
quantum evolution of the gluon field with y is captured by the non-linear JIMWLK evo- 
lution equation [256-260]. Using a Gaussian approximation for the weight W, known as 
McLerran-Venugopalan model, gluon production in p + A collision can be interpreted as 
multiple 2 — > 1 partonic scatterings [261]. A comparison to the 2 — > 2 multiple scatterings 
included in the GE models is discussed in Ref. [51]. 

When the correlations in the nuclear gluon field are local (i.e. in the incoherent scattering 
approximation), the gluon /77-spectrum shows a peak structure which increases in mag- 
nitude and moves to higher pr with increasing rapidity [262]. This agrees with RHIC 
data at midrapidity, but is at variance with forward rapidity measurements. In the deep 
saturation regime, the solution of the JIMWLK equation describes high density partons 
with non-local correlations, and suppresses gluon production over the whole pr range, 
irrespective of the rapidity (see Fig. 40). The resulting picture [263, 264] is that at the 
values of x probed at RHIC at midrapidity, the nucleus wavefunction has not yet reached 
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Fig. 40. - Cronin effect in ttie CGC formalism with naive quantum evolution (left) and with full evolution 
in the "deep saturation" limit (right). Figure taken from Ref. [262]. 



the saturation regimes - is small, at most 2 GeV^ - hence a description in terms of 
pQCD multi-scatterings is vaUd (and perhaps numerically more accurate [51]). At larger 
rapidities, the saturation scale increases and the nucleus wave function undergoes a longer 
quantum evolution, resulting in hadron suppression on a large pj interval. Hence, gluon 
saturation may have indeed been revealed in the observed forward rapidity hadron sup- 
pression at RHIC. 

To make quantitative contact with phenomenology, one approximates the proton as a dilute 
colour source, so that gluon production in p + A collisions can be explicitly written in a 
^r-factorised form [262]. The nucleus wavefunction is computed in terms of a colour 
dipole forward scattering amplitude, which is modeled in order to capture the analytic 
key features of the CGC, and to incorporate the transition to the semi-classical regime at 
lower rapidity [127, 245, 246, 265]. Examples from the computations of Refs. [245, 246] 
are shown in Fig. 39 (right). 

An important limitation of parton multi-scattering models is their inability to explain the 
strong flavour dependence of the Cronin effect and the corresponding "baryon anomaly" (tt"^ < 
<s: p,p at intermediate pr) observed at RHIC. Indeed, if the pr-broadening has a partonic 
origin, the only flavour dependence of the Cronin effect can be due to the different contribution 
of gluons and quarks to the final-state hadron, which cannot explain the large difi"erence between 
protons and pions [49]. 

6' 1.2. Final state parton recombination. - A description of hadronisation as parton recom- 
bination, as opposed to the parton fragmenting into the observed hadron, was proposed long 
ago to explain hadron production at large Feynman xp m h + A collisions [266]. The idea 
that two or three partons can recombine, or coalesce, into a meson or baryon has been re- 
vived [125, 126, 183-185, 267-269] because of recent experimental findings at RHIC, chiefly, the 
strong baryon enhancement (already observed inh+A collisions at Fermilab, see Section 4), and 
the scaling of the elliptic flow of intermediate pr hadrons with the number of their constituent 
quarks [270,271]. The main idea is that whenever there is a large reservoir of partons in the 
final state, there is a lesser need to produce additional ones through parton splitting as assumed 
in parton fragmentation. Recent reviews of recombination models can be found in [233, 272]. 

In heavy-ion collisions at RHIC, a hot and dense medium made of deconfined soft quarks 
and gluons is created in the first instants after the collision. The parton reservoir is made of 
these thermal partons, which are distributed in momentum according to an exponential spectrum. 
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Fig. 41. - Left: Ratio of baryon yields to meson yields for central Au + Au collisions at RHIC 
i/iAW = 200 GeV. The recombination model computations by Greco, Ko, Levai [185, 267] and by Fries, 
MuUer, Nonaka, Bass [184,268] for pjn- (left) and A/2A:^ (right) are compared to PHENIX and STAR 
data. Plot taken from Ref. [233]. Right: Hwa-Yang recombination model computation [125] of the Cronin 
effect for midrapidity pion and proton md+Ait collisions at RHIC ^[snn = 200 GeV, compared to PHENIX 
data [1 19, 122]. Plot taken from Ref. [273]. 



~ e ''/^ with T X 0.2 GeV. At large momenta partons come dominantly from hard interactions 
and are distributed according to a power law, ~ l/Pj, with « ^ 7 - 8. Parton fragmentation 
requires the production of a hard parton with large momentum pj- > pj to fragment into the 
hadron. On the contrary, parton recombination is kinematically favoured since 2 or 3 partons of 
smaller momenta pr < Ph can coalesce into the observed hadron. At large pj, the fragmentation 
dynamics is supposed to dominate over recombination since thermal partons are suppressed, 
while recombination is expected to be important at low pT- In an intermediate pt range, the 
recombination of one hard parton together with one or two thermal partons becomes the dominant 
channel to produce a meson or a baryon. Recombination readily explains the Cronin effect of 
mesons and baryons at intermediate pj'- the larger the number of soft thermal partons involved 
in the recombination process, the larger the momentum broadening, leading to a stronger Cronin 
effect for baryons than for mesons. As can be seen in Fig. 41 left, these models turn out to 
reproduce well the data. 

In /:> + A collisions, the parton reservoir is made of the soft partons created in sequential 
nucleon-nucleon scatterings [125]; in particular, they can recombine among themselves into the 
observed low-p^ hadrons. Fitting the exponential soft parton spectrum to the observed yield of 
I0W-/77 hadrons in p+p collisions, this recombination model can explain well the large difference 
in the baryon and meson Cronin effect ind+Au collisions at RHIC (Fig. 29). Additionally, it can 
describe the observed suppression of intermediate pj hadrons in forward rapidity bins, and also 
the negative/positive rapidity asymmetry of the Cronin effect at moderate 77 because its magnitude 
tracks the (soft) hadron rapidity distribution, which decreases in rapidity at > -2 [274]. The 
model is quite successful at RHIC energy in the range -0.75 < 77 < 3.2. An interesting check 
would be to measure hadron pj spectra at < -2, where the soft hadron yield starts decreasing 
with decreasing 77. Therefore, the magnitude of the Cronin effect should peak aX tj ^ -2, and 
slowly decrease as 77 is reduced. 

61. 3. Origin of the Cronin effect. - In the parton recombination model, the Cronin effect is a 
final-state (FS) effect at the hadronisation stage rather than due to initial state (IS) or final-state 
parton rescatterings in the target nucleus. However, contributions to the Cronin effect from IS 
parton rescatterings cannot be entirely discarded. In fact, in e + A collisions the multiplicity of 



PARTON PROPAGATION AND FRAGMENTATION IN QCD MATTER 



57 



soft hadrons in the final state is much lower than in p + A collisions, making the recombination 
mechanism much less effective, even though a large Cronin effect is observed in the HERMES 
and EMC data, see Section 3'2. (Another potential difficulty of the recombination models in 
nDIS would be to account for the anti-proton Cronin effect which is similar to that of the mesons, 
unUke in h+A colUsions where p « p.) The same processes which cause the Cronin effect in e+A 
collisions are also likely to be at work in h +A collisions in addition to possible PS recombination 
effects. Indeed the hard partons are traversing two media, the cold nuclear target and the soft 
partons cloud created in the collisions, and both can contribute to the Cronin effect. Since these 
two media are well separated in space-time because of the formation time oc 1 /Aqcd for the soft 
partons, multiple parton scatterings may add to the -broadening from FS recombination. 

The effect of a short quark life time, as indicated by HERMES data in nDIS, has not been 
included so far in recombination model computations. At RHIC energy the large DIS -equivalent 

V = 1^(1000 GeV) (see Section 2'2) ensures that parton fragmentation starts outside the cold 
nuclei. However, at small enough pj, fragmentation may still start inside the soft parton cloud, 
thus reducing the quark path-length in the soft cloud, therefore reducing the probability of soft 
parton pick-up and the effectiveness of the recombination mechanism. At fixed-target energies 

V = 10-100 GeV hadronisation may in fact start inside the cold nucleus target, before the parton 
has any time to travel through the soft parton cloud, thus preventing altogether the recombination 
mechanism. 

We finally remark that a precise understanding of the mechanisms underlying the Cronin ef- 
fect in h + A collisions is crucial for measurements of hadron quenching in A H- A collisions at 
the low SPS energy, a/^aw = 17.3 GeV [161]. In A -H A collisions, two competing effects modify 
hadron pj spectra: (i) the pr broadening caused by both targets, hence larger than inh+A col- 
lisions, and (ii) the hadron suppression due to parton energy loss in the produced medium. The 
Cronin effect is observed inh + A collisions to grow with decreasing ^Jsm = 200 - 27 GeV: the 
steeper transverse momentum spectra at lower energies translate into a larger relative enhance- 
ment of the yields for the same kx "kick". Hence, at SPS energy one can expect an even larger 
enhancement than observed at Fermilab. Therefore the Pb + Pb/p + Pb ^ 0.5 nuclear modifica- 
tion ratio observed by the WA98 collaboration [161] might be caused by a large parton energy 
loss mPb + Pb that outweighs a large Cronin effect [159, 161]. A numerical evaluation of jet 
quenching at fixed-target energies requires a precise theoretical and experimental control of the 
underlying Cronin effect (including basehne p + p pr spectra at SPS energies). 

6'2. Energy loss in hot QCD matter . - 

6'2.l. Formalisms. - The idea of parton energy loss in hot QCD matter was first discussed by 
Bjorken in the early eighties [131]. This process was then revived a decade later when for the first 
time Thoma and Gyulassy [132] and Gyulassy, Pliimer, and Wang [275,276] computed pertur- 
batively the radiative energy loss of high-energy quarks in a QGP. Since then many approaches 
have been developed to determine the gluon radiation spectrum, dl/dw, of a hard parton under- 
going multiple scattering. For a short discussion and comparison of the different energy loss 
formaUsms, see Ref. [165]. Let us briefly recall the main assumptions made in each of the four 
main frameworks detailed there: 

• Opacity expansions (BDMPS and GLV) 

BDMPS developed the first perturbative framework to describe the medium-induced gluon 
emission process from the soft multiple scattering of hard partons in both cold [277] and 
hot [278] QCD matter. The calculation assumes that the number of collisions, or opacity, 
is large: n = L/A » 1, where L is the medium length and A the parton mean free path. The 
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typical momentum exchange in each scattering is given by the Deb ye mass of the medium, 
H, and the hard scale of the calculation is n x //^ = §L » A^^^, where the transport 
coefficient q = fi^/A represents the scattering power of the QCD medium. The Landau- 
Pomeranchuk-Migdal (LPM) effect - basically the destructive interference of the gluon 
radiated on several scattering centres - takes place whenever the gluon lifetime t = tj/kj 
exceeds its mean free path, that is for gluon energies larger than ~ jj^A. Soft gluon emission 
is assumed, a> ~ qL^ « E, although corrections (co/E) were proposed in [108]. On the 
contrary, the formaUsm by Gyulassy, Levai and Vitev (GLV) first took into account one 
hard scattering in the medium [279] (i.e. first order in n) and from this a recursive approach 
has been used to determine the gluon spectrum at any opacity [60, 138]. Such an expansion 
has also been rederived by Wiedemann in Ref. [280] from the light-cone path-integral 
(LCPI) approach to energy loss developed by Zakharov [143, 230] (and equivalent to the 
BDMPS framework described above, see [141]). All these approaches model the medium- 
modified fragmentation functions as a convolution of (vacuum) fragmentation functions, 
D'^iz, Q^), and an energy-loss probability distribution, ^{e) [281]: 



is the rescaled momentum fraction carried away by the hadron h in presence of the QCD 

medium, and Ef is the energy of the parent parton in the medium rest frame. In heavy- 
ion collisions, for instance, Ef - kj for partons produced at mid-rapidity, while in semi- 
inclusive DIS to LO in the strong coupling constant the quark energy is given by the photon 
energy, Ef - v. The energy loss probability distributions <^(e), also called quenching 
weights, will be discussed in Section 6'2.2. 

Higher-twist formalism (HT) 

In the HT approach [92,282,283], power corrections proportional to IjQ^ (where Q is 
the typical radiated gluon virtuahty) and enhanced by the medium length L are included 
to the leading-twist total cross section in DIS, assuming that Q is much larger than the 
energy-scale exchanged jj. between the hard parton (with energy E) and the medium: E » 
Q'S> fx. This approach was first proposed in the context of DIS (see Section 6'3.2) where 
the strength of the higher-twist terms is monitored by a single parameter, C, which can 
be adjusted to fit the data just like q in BDMPS or the initial gluon density dN^ jdy in 
GLV. For hot QCD media, a phenomenology based on a z-shift prescription for the FF, 
z z/il - e/Ef) similarly to Eq. (45), was developed afterwards [282]. The energy loss 
€ can be related to C and is adjusted to fit the data. 

Thermal field theory (AMY) 

The formalism of Arnold, Moore, and Yafie (AMY) addresses the production of thermal 
photons in a finite-temperature QCD medium [284,285] and was then extended to describe 
gluon emission from the scattering of "hard particles", with energy (T), on softer modes, 
€>" igT). The photon/gluon emissions that are collinearly divergent were resummed to all 
orders in , leading to a suppression as compared to the leading-order result because 



(45) 




where 



z* = z/(I - e/Ef) 
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of the LPM effect already discussed. This calculation should be accurate for asymptoti- 
cally large temperatures, for which all the relevant (hard, soft, ultra-soft) scales are well 
separated: T » gT » g^T . 

An alternative Q^-shift prescription for medium modified fragmentation functions - which 
is quite different from the z-rescahng prescription (45) used in most apphcations - has been 
discussed in [29, 70]. Assuming that hadronisation occurs outside the medium, as in the above- 
mentioned energy loss models, an upper bound on hadron quenching can be obtained by consid- 
ering 

(46) D^(z,e^)=Dj(z,e2 + Ap2), 

with Ap^ the parton in-medium pr-broadening related to q [140]. Numerical estimates show 
that such an effect is much smaller than observed in experimental data, which is interpreted 
as showing that mechanisms other than energy loss are at play in quenching hadron spectra. A 
detailed vacuum hadronisation model which tries to avoid this ambiguity [29] is briefly discussed 
at the end of Section 2 4, and its consequences for hadron quenching in cold nuclear matter will 
be reviewed in Section 7'3. 

6'2.2. Quenching weights. - The quenching weights appearing in Eq. (45) have been com- 
puted numerically within the BDMPS approach by Salgado and Wiedemann (SW) for light par- 
tons [286,287], and later on extended by Armesto, Salgado and Wiedemann (ASW) to include 
heavy quarks [288,289]. They have also been studied in the GLV formalism in Ref. [290]. 

The ASW quenching weights are determined using the Poisson approximation of independent 
gluon emission [291], 

It depends on the characteristic radiated gluon energy ojc = qL/2 and on the medium size param- 
eter R = (o^L, where L is the parton in-medium path-length. The approximation of an asymptoti- 
cally large medium, considered by BDMPS (supplemented in [59] with finite energy corrections 
of order ^ ^ao/Efj) would correspond toR oo. At finite R, the parton has a probabihty po of 
not interacting with the medium and therefore not to suffer any energy loss. Correspondingly, 
the quenching weight can be spUt in a discrete and continuum parts, 

(48) l3^ie)=po6ie) + Pie). 

The quenching weight is computed for a static and uniform medium. In heavy-ion collision, the 
longitudinal expansion of the medium is taken into account by rescaUng the transport coefficient 
according to an approximate scaUng law discussed in [292]: 



(49) q(L)^^ driT- to) m 

where a characterises the power-law time-dependence of the medium number-density, p(t) oc 
T~", and §(t) = §(to) (to/t)". The purely longitudinal (or Bjorken) expansion corresponds to 
a - I, and is often assumed in phenomenological applications. When fo <c L, Eq. (49) reduces to 
(q) ^ 2q{to) to/L [293]. In nDIS, the medium is static but non-uniform, and an analogous scaling 
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law is proposed in Ref. [43]. Recently, a simple prescription has been given by Arnold in order to 
determine dl/dw in a finite and expanding medium [294]. Therefore, applying this recipe using 
e.g. hydrodynamical space-time evolution will allow eventually for the computation of more 
realistic quenching weights. Also, Peshier has proposed a useful way to compute iteratively the 
quenching weights in Ref. [295] avoiding the need to compute inverse Laplace transforms (as 
done e.g. in [59,286]). 

6'2.3. Phenomenology. - Many phenomenological applications of the above formalisms have 
been carried out in order to describe or predict the production of large-/?7 hadrons or jets in 
heavy-ion collisions (see Section 5 and Ref. [137] for a recent review). All models coincide 
in characterising the system with very large initial gluon densities {dN^ jdy ^ 1400), transport 
coefiicients {q) ^ 13 GdV^/fm and/or high temperatures T x 400 MeV, in order to reproduce 
the existing high-pr suppression. Yet, when trying to compare all model predictions through 
a common q coefficient, differences of a factor of 3 - 4 appear [164] . Recently, progress has 
been made towards a more reaUstic implementation of energy-loss scenarios including a fuU hy- 
drodynamical expansion of the produced medium - constrained by soft/global observables; first 
in [296-298] and more recently in [299-301]. These approaches aim at a consistent description 
of soft and hard probes and at a realistic extraction of the medium-parameter from fitting the 
large-/7r suppression data (see e.g. Fig. 31). A recent quantitative comparison of energy loss 
schemes under identical conditions (i.e. same medium evolution, same choice of parton distribu- 
tion functions and scales, etc.) is presented in [302]. 

6'2.4. Modified DGLAP evolution. - As an alternative to the discussed energy loss formalisms 
based on an energy-rescaling of vacuum fragmentation-functions, attempts to reformulate par- 
ton energy loss in pQCD as a medium modification of the Dokshitzer-Gribov-AltareUi-Parisi 
(DGLAP) evolution of fragmentation functions have been recently suggested [303-305]. 

In the perturbative description of fragmentation processes, the produced parton of time-like 
virtuality Q radiates gluons in order to reduce its virtuality down to a soft scale Qo = GeV) 
where hadronisation takes place. The Q^-evolution is governed by the DGLAP equations [1-3], 
which control the probability that a quark, say, branches into a quark and a gluon in going from 

+dQ^ to Q^. Gluons also can spht into a or a qq pair. The offspring partons can in turn spUt, 
iteratively. After each branching the scattered parton gains some pj, which can be controlled by 
using pr-dependent evolution equations [306]. Physically, the probabilistic picture of DGLAP 
evolution leads to a parton shower which can be simulated in a Monte Carlo generator. 

In the presence of a QCD medium, the parton shower may be modified in basically two 
ways: (i) the splitting probability is enhanced, mocking radiative energy loss, and (ii) the partons 
can rescatter, leading to a stronger showering process and a pr-broadening of both the leading 
particle and the shower. Various models for the medium-modified spHtting enhancement have 
been proposed. In [304], under the assumption that splittings are independent, the in-medium 
spUtting functions are 

(50) P{z) = P'^\z) + ^P{z) , 

where z is the fractional momentum carried away by one of the split partons, and f ^^'^ is the vac- 
uum splitting probability. The additional term Af can be calculated from the medium induced 
gluon radiation spectrum, dlldio, previously discussed. Medium-modified fragmentation func- 
tions at a scale Q are computed by evolving input fragmentation functions at a scale Qo according 
to the DGLAP equations with medium-modified splitting functions, see Fig. 42 left. A virtue of 
this approach is that at large virtualities, Q » A , the usual medium-modified fragmentation 
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Fig. 42. - Left: Medium over vacuum ratio of the gluon fragmentation functions plotted for a medium witli 
q= \Q GeV^/fm (green) and ^ = 50 GeV^/fm (red), and for two different medium lengtlis: L = 2 fm (solid) 
and 6 fm (dashed). Figure taken from Ref. [304]. Right: Single inclusive hadron distribution as a function 
of ^ = log(£jet/p). e*e^ data from TASSO [307] and OPAL [308] compared to vacuum-FFs (solid curves) 
and to medium-FFs (dashed/dotted curves, obtained with /i^ej = 0.8 in the LO splitting functions). Figure 
taken from Ref . [303]. 



functions defined by the quenching weights, Eq. (45) is formally recovered. In the approach of 
Borghini and Wiedemann [303], the medium effect is argued to modify the splitting functions by 
enhancing its singular part; for example, the quark splitting function is 

(51) f,,(z) = C^(^^^-^-(l+z)), 

where f,„ed - gives back the vacuum An interesting application is the medium mod- 

ification of the single-inclusive energy distribution of hadrons inside a jet, dN/d^, where ^ - 
log(£'jet/£/i). This leads to a distortion of the usual hump-back plateau, predicted within the 
Modified Leading Logarithmic Approximation (MLLA) of QCD and observed experimentally 
in e^e , DIS and hadronic collisions (for a review see e.g. [309]). As can be seen in Fig. 42 
(right), the number of highly-energetic particles (small ^) is suppressed while the soft gluon 
yield is enhanced at large ^ due to energy-momentum conservation. The inclusion of 2 — » 2 elas- 
tic rescatterings, which accounts for elastic energy loss, can be accomplished by supplementing 
the DGLAP evolution equations with a gain and loss term, which describes partons scattered into 
and away from a given kinematic variable bin, and is also suitable for a Monte Carlo interpreta- 
tion [305]. 

Several parton showers in the medium have been recently developed [310-314]. As com- 
pared to analytic calculations, parton showers have many advantages such as conserving energy- 
momentum throughout the evolution. They allow one to directly compare their multiple-differential 
hadronic distributions to experimental data, and thus give a better access to the microscopic dy- 
namics. They also allow to study the particle and energy flow inside a jet. As an example, 
first results from the JEWEL parton shower [310], indicate that the distribution of 1, 2 and 3 
jets events (reconstructed using a given granularity parameter, jcut) is sensitive to the elastic or 
inelastic nature of parton rescatterings. 

Applications of modified DGLAP equations have been studied only in the context of A -H 
A collisions (see, e.g., [215,315]) testable at RHIC and LHC. It would be very interesting to 
also study jet modifications in e -H A collisions, which would be accessible at the Electron-Ion 
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Collider (Section 87.4), and the hadron pj -broadening already under study at HERMES and 
CLAS (Section 8' 1.3), although the lower virtualities of the latter limit the amount of QCD 
evolution (parton radiation) accessible. 

63. Energy loss in cold QCD matter . - 

Although the density of scattering centres in cold QCD matter is much smaller than in a de- 
confined state of matter such as a quark-gluon plasma, it is not a priori excluded to observe effects 
of parton energy loss in large nuclei [277]. As a matter of fact, several of the above-discussed 
parton energy loss formalisms should prove more appropriate to describe the rescattering of a 
hard parton off static nucleons rather than off quarks and gluons carrying thermal momenta, 
p = (T). Moreover, the nuclear density is well known and the medium does not expand while 
the hard parton traverses it, unlike for the QGP whose energy density drops rapidly with proper 
time. For these reasons, cold QCD matter is an ideal testing ground to compare different energy 
loss formalisms and test the approximations made in their phenomenological applications. 

One pragmatical approach, consists in comparing the nuclear DIS and h+A hadroproduction 
data to NLO pQCD calculations (with nuclear PDFs) and simply encode any observed modifica- 
tions of the final yields for various species in properly fitted medium-modified parton-to-hadron 
FFs [316]. More commonly, however, one tries to get a more physical insight on the mechanisms 
affecting hadron production in cold nuclear matter by describing parton multiple scatterings in 
nuclei within the BDMPS framework [277], as done for hot and dense QCD matter. In such an 
approach [278] one derives a simple relationship between the energy lost by the hard parton (per 
unit-length) and its transverse momentum broadening, [141] 

(52) -- . — ip% 

independent of the parton species. The transport coefficient of cold nuclear matter, q - IJ^/A, 
governing the amount of energy loss or momentum broadening can be estimated perturbatively. 
It is related to the nuclear matter density, p ^ 0.17 fm"^, and the gluon density, G, in a nu- 
cleon [277]: 

(53) r'" = ^^^pxG(x,e^) 

where Cr is the colour charge of the parton, (AT^ - l)/2Nc and Nc for quarks and gluons respec- 
tively. The value x at which G should be evaluated is parametrically (niNty^, where the scale 
{ is larger than the mean free path of the hard parton. A, and smaller than the medium length, 
L. The virtual scale which enters G is of ^ (§L). Using - 1/2 at such low scales, and 
xG{x, Q^) ^ 1, the gluon {Cr = Ca = Nc) transport coefficient is roughly given by 

(54) ^ 0.05 GeV^/fm . 

This value is much smaller than the leading-order BDMPS estimate for a hot plasma at RHIC 
temperatures, ^ 2.2 GeV^/fm, given in [144]. 

When discussing energy loss in cold nuclear matter, however, one should make an impor- 
tant distinction between initial-state and final-state energy loss [141, 317], because the medium- 
induced gluon radiation can interfere with the radiation originating at the hard scattering (see 
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Fig. 43. - Sketch of parton propagation in cold nuclear matter in the nucleus rest frame, (a) Initial- and 
final- state interactions in /i + A collisions in the nucleus rest frame, (b) Absence of initial state interactions 
in / + A collisions, (c) Absence of final-state interactions in DY events, (d) The theoretical case of an 
asymptotic parton penetrating the nucleus with no hard interactions. The nucleus is drawn as an oblong oval 
for illustrative purposes. The possibility that hadronisation starts inside the target nucleus is considered in 
this cartoon, see Section 7. 



Fig. 43). Inh+A and A+A collisions, both initial- and final-state interactions occur. If the rapid- 
ity of the parton which fragments in the observed hadron is far enough from projectile rapidity, 
initial-state and final-state radiation should not interfere, and may be treated independently. In 
semi-inclusive DIS, only final-state interaction may take place because the projectile lepton does 
not interact strongly with the nucleus. Similarly, in Drell-Yan lepton pair production only initial- 
state interactions occur; an advantage of the Drell-Yan processes is that conversion rate to a 
lepton pair is perturbatively calculable, unlike the non-perturbative fragmentation function con- 
verting a parton into an observed hadron. The theoretical important case of an asymptotic parton 
penetrating the nucleus with no hard interactions is also considered in the figure 43 (d). Because 
of the interference with the hard radiation, each of these energy-loss problems should be con- 
sidered and solved separately, in order to probe the properties of cold nuclear matter such as its 
transport coefficient q via measurements of the parton energy loss -dE/dz- 

Baier et al. give in Ref. [141] the expression for the energy loss of an asymptotic quark: 

,,,, dE a, Nc , 

which is a factor 3 smaller than for a quark produced inside the medium [277]. Indeed, a quark 
coming from -oo has had the time to construct its gluon field and therefore will only start to 
radiate (medium-induced) gluons only after it has experienced a first single scattering in the 
medium, after a time t - {A). On the contrary, a quark produced in the medium immediately 
radiates soft gluons to get rid of its virtuality; this gluon emission interferes constructively with 
gluon emission stimulated by the medium. 

In Ref. [317], initial- and final-state energy loss in cold QCD matter has been evaluated in 
the GLV formalism. The resulting fractional energy loss is shown in Fig. 44 as a function of 
the parton energy fijet, for two sets of average momentum transfer squared fi^ and parton mean 
free path A coiTesponding to a fixed value of q - i-P'/A - 0.12 GeV^/fm, and compared to 
the asymptotic parton case. There are large differences in initial- and final-state energy losses, 
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Fig. 44. - Fractional energy loss Ap* /p* for massless quark partons versus the parton energy £jet, in cold 
nuclear matter of length L = 5 fm. Two different sets of typical momentum transfer per scattering and gluon 
mean free path have been used for comparison. Plot taken from Ref. [317]. 



which are process-dependent and need to be correctly accounted for to experimentally access 
the properties of the nuclear medium. Furthermore, a single parameter such as the transport 
coefficient q - fP'/A may not adequately describe the stopping power of cold nuclear matter, as 
shown by the different energy loss of the solid and dashed lines. 

Keeping this in mind, we review in the following the existing phenomenological applications 
of energy loss computations in cold nuclear matter for semi-inclusive nuclear DIS and DY mea- 
surements. Table IX and Fig. 49 summarise various estimates of -dE/dz from the existing data 
(see also the short reviews in Refs. [318, 319]). 

63. 1. Nuclear DIS in the BDMPS formalism . - In the BDMPS approach of Ref. [59, 320], 
extended in [43, 321] to include a realistic treatment of the nuclear geometry, the reduced quark 
energy at the time of hadronisation is translated into a shift of Zh in the vacuum fragmentation 
function D ^ D [281] via quenching weights. At leading order (LO) in perturbative QCD, the 
hadron multiplicity is then computed as follows: 



where D is the medium-modified fragmentation function, Eq. (45), computed with Ef = v. 

The LO computation is known to underestimate the experimentally measured average (v)^ 
and in each Zh bin [64]. The problem, is likely to be solved at NLO. An effective way of 
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circumventing it at LO is to approximate Eq. (56) by 

The value of the average variables {x)z and in each z/,-bin is taken from the measured 

values. This procedure is used in all computations based on LO cross-sections [29, 58, 61, 320, 
321]. The only parameter of the computation, namely the transport coefficient at the centre of 
the nucleus, is found to be 

(58) q = 0.6 GeV-/fm 

in order to reproduce the latest pion quenching data on the Kr target from HERMES [23], see 
Fig. 45. When comparing the theoretical v-distribution with data, mind that (z/,)v ■x, 0.3, where 
the model tends to slightly overestimate the data. Note that the extracted cold matter q value is a 
factor 10 larger than the perturbative estimate q x 0.05 GeV^/fm from Baier et al. [277]. 

An appropriate treatment of the medium geometry is important to extract the quenching pa- 
rameter from the data [43, 321]. In Fig. 46, we compare experimental data for pion production 
on Kr to realistic and approximate geometries, using SW quenching weights. The crudest ap- 
proximation is to use an average quark path-length L x {3/4)Ra, asymptotic quenching weights 
(R — > oo) and a constant nuclear density, corresponding to a constant transport coefficient q. 
This approximation is commonly considered, but cannot reproduce the z-dependence of the data 
(dotted line) even if § » 0.15 - 0.20 GeV^/fm is adjusted for the curve to touch the data. 

We can improve these approximations in several steps: using asymptotic quenching weights 
with a variable medium length (dashed line), non-asymptotic quenching weights using a finite 
R - oJcL, either with fixed average L (dot-dashed line), or with variable L depending on the 
position of the 7*-quark interaction point (solid line). The largest effect is given by the use 
of non-asymptotic quenching weights, mainly because of a non-zero probability of no energy 
loss, see Section 6'2.2. Modelling the full geometry instead of assuming a mean length gives 
a smaller but still important effect. In particular, when geometry is properly taken into account 
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Fig. 46. - Realistic and approximate treatments of the nuclear geometry with quenching weights in the 
model of Ref. [43, 321]. The upper two lines are computed with R = ai^L, the lower two lines with R ^ oo. 
Solid and dashed lines: variable quark path-length. Dotted and dot-dashed lines: fixed L = (4/3)Ra- 
Experimental data are for quenching of /r^ (full circles) and tt" (open circles) on Kr from HERMES [23]. 



there seems to be no need to invoke a finite quark lifetime to explain the large-z data, as proposed 
in Ref. [320]. 

Energy conservation, namely the constraint e < Eg, is not always fulfilled in the quenching 
weights because of the approximations involved in their determination. It can be imposed from 
the outside by cutting the quenching weight ai e - Eg and reweighting it to conserve probabil- 
ity [164]. Alternatively, one can cut the single gluon radiation spectrum aio) - E, and consider 
ff{(ji)clEq) corrections [59]. These corrections tend to reduce the quenching. Their effect in- 
creases with q and A, and decreases with z/,. Their magnitude is about 5-10% for Kr targets, 
and typically comparable to the experimental error bars in the figure. Thus they are subleading 
compared to the effect of correctly implementing the nuclear geometry. 

The flavour dependence of hadron quenching in the BDMPS formalism has been discussed 
in [320,322,323]. We can approximate the vacuum FF at large zi, by D^izh) (1 - ziif% and 
read the value of the exponents from the global fits of Refs. [324] (Table 2 there). Taking into 
account that n"^ production is dominated by m, d quarks, and K"^ production is dominated by m, s 
quarks, at 2^ = 2 GeV^ we have 

(59) Kd^Zh)'^{l-Zh)' 

(60) DX,{zk)^{l-Zhf\ 

So at HERMES (Q^ » 2.5 GeV^) we may expect a slightly stronger suppression for large-z/, 
pions than for kaons (this is opposite to what is asserted in Ref. [29], which assumes fi" « 0.5 
and a; 0.8 on the basis of Regge phenomenology). Using the flavour separated global fit of 
Kretzer [325] one has fi^ < fi'^ , which translates into a stronger quenching for K than for , 
compatible with data. The described flavor dependence of hadron attenuation in the BDMPS 
formalism is in agreement with HERMES data. 

6'3.2. Nuclear DIS in the higher-twist formalism . - In Refs. [58,282,283] the medium modifi- 
cations of the fragmentation functions are computed from leading twist and twist-4 contributions 
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(c) (d) 

Fig. 47. - Left: Gluon radiation from a single scattering (a) and from double scattering (b-d). Right: a 
sample diagram for quark-gluon rescattering processes with 3 possible cuts. Figures taken from [282]. 



to the leading order DIS cross-section, including diagrams with one elastic quark-nucleus scat- 
tering and one radiated gluon, see Fig. 47. Both the struck quark and the radiated gluon are 
allowed to fragment according to vacuum FF. The obtained modified FF, D, can be modeled to a 
good accuracy by shifting zh in the leading-twist fragmentation function 

(61) Dizh) -t-^dI 



l-Azh \l- Azh 

where Azh - AEq/v is the quark's fractional energy loss, and AEq 0.6{zg)v with the average 
fractional energy (zg) carried away by the radiated gluon computed diagrammatically [58, 282, 
283] 

(62) {zg} ^ a?(G-) C(e') RI ^ fgil - Zh) , 

where w^v is the nucleon mass and Ra the nuclear radius, and fg is a function of (1 -z/,) because of 
energy conservation [326]. The average gluon energy depends on one parameter, C{Q^), which 
represents the strength of parton-parton correlations in the nucleus. We can note the dependence 
on the square of the medium size, typical of the LPM effect in QCD. Inclusion of quark-quark 
double scatterings beside quark-gluon double scatterings leads, except for pions, to a different 
quenching of positively and negatively charged hadrons [327]: 

nK- , nK* . nP ^ jfP . p/l" ^ nil* 
Rm <Rm ' ^M^^M' < Rm 

These features agree qualitatively with HERMES data [22,23,44], and with the BDMPS model 
previously discussed. A generalisation of the higher-twist formalism to include heavy-quark 
energy loss has been discussed in [328, 329]. 

To apply the model to HERMES data, the parameter C is fitted to the overall suppression of 
unidentified charged hadron on a nitrogen target: C(2^) = 0.0060 GeV^, with OsiQ^) - 0.33 at 

= 3 GeV^, which corresponds to an energy loss dE/dz = -0.5 GeV/fm for a Am nucleus [58], 
or equivalently to a transport coefficient q - 0.12 GeV^/fm [319]. The multiplicity ratio for 
other targets and parton species can then be computed without further adjusting the parameter, 
see Fig. 48. See also Ref. [93]. One observes an overall agreement with data, but with a tendency 
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to overestimate the slope of Ruizh) for heavy nuclei. As discussed above in connection with the 
BDMPS approach, this can be caused by a too schematic treatment of the medium geometry. 
Another important remark, is that the computations of Refs. [58, 93] include only twist-4 dia- 
grams, i.e., up to one parton rescattering in the medium. This might not be sufficient for large 
targets such as Kr. A computation up to twist-6 (two rescatterings) has been carried out in [330], 
however without comparison to experimental data. (An all-twist resummation for radiative pro- 
cesses is unfortunately very hard to achieve, but the first preliminary steps have been taken in 
Refs [331, 332].) Higher-twist effects on the DGLAP evolution of the fragmentation functions 
have been studied in Ref. [333], and reduce the slope of Rm at large Zh- We should finally note 
that coherent multiple parton scatterings may lead to an additional suppression of quark produc- 
tion [334], which would also reduce the slope of Rm- 

63.3. Drell-Yan processes . - Since the lepton pair does not interact strongly with the nu- 
cleus, the Drell-Yan process offers a clean probe of the initial multiple scattering of the projectile 
parton (a quark to leading order in a,) in the target nucleus before the hard process, qq (^{~ , 
takes place. The effects of energy loss on the Drell-Yan cross section can be estimated from its 
production cross section at leading order. 



(63) , , = — — 2j 1«'7p(^i)^<m(jc2) + qiip{xx)qiiA{x2)\ , 



where the fractional momenta xi 2 are related to the observable Feynman Xf and dilepton mass 
M as in Eq. (18) (see Section 21). At large xi > 0.5 and for X2 = 0.1 - 0.3, where quo ~ qt/p, 
the ratio of the DY cross section on a nucleus A and on deuterium D can be approximated by 



(64) 



^^Ja q""'(^i) (l-xf)" 



dx\dx2l dx\dx2 qu/pixi) (1 — xi)'' 



where the large-x exponent is 77 ^ 3 from quark counting rules [318]. Assuming that the projectile 
quark experiences an initial-state energy loss per unit length, a = -dE/dz, we can set x^ = 
xi + a{L)A I E, where E is the quark's energy, for the sake of a qualitative discussion. Hence, one 
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Reference 


-dEldz (GeV/fm) 


Observable 


[283] 


~ 1.2 


Nuclear modification oie + A fragmentation functions 


[58] 


-0.5 


Nuclear modification ofe + A fragmentation functions 


[277] 


-0.4 


Pt broadening of p + A D-Y yield 


[335] 


0.20 ±0.15 


Nuclear dependence of 150 GeV n + A D-Y cross sections 


[100] 


< 0.44 


Nuclear dependence of 800 GeV p + A D-Y cross sections 


[336] 


1.12 ±0.15 ±0.21 


Nuclear dependence of 800 GeV p + A D-Y cross sections 


[337] 


0.95 ± 0.21 ± 0.21 


Nuclear dependence of 800 GeV p + A D-Y cross sections 


[338] 


-2.8 


Pt broadening of p -I- A jets 



Table IX. - List of published parton energy loss values in cold nuclear matter extracted from various e + A 
and h + A observables (adapted from Ref. [318]). 



obtains 



^ ^ (1X1(1X21 dx\dX2 iip(l-X^)' 

where Ep is the incident proton's energy in the nucleus rest frame. This equation shows that the 
energy loss increasingly suppresses the DY cross-section the larger the nucleus, and the larger 
xi, or equivalently the larger xp = xy- X2 when X2 ^ 1. A related observable is the dilepton pj 
broadening, which is related to the /j^-broadening of the incoming quark caused by the parton 
initial state multiple scatterings (see Eq. (52)). 

The wealth of experimental data, reviewed in Section 4'1, makes likely a precise estimate of 
the cold nuclear matter transport coefficient q'^°^'^, or more generally the amount of energy lost by 
fast quarks in heavy nuclei via expressions like Eq. (65). Note, however, that the nuclear mod- 
ifications of the parton distribution functions (shadowing, an ti- shadowing and the EMC effect 
depending on the typical values of x) [20, 109-1 1 1] play a role in the nuclear dependence of DY 
production, and would need to be better constrained. 

In Ref. [339], the E772 and E866/NuSea data have been analyzed as a function of the mo- 
mentum fraction xi as well as a function of the mass of the Drell- Yan pair, M^^_ . The DY process 
is calculated using a dipole approach in the rest frame of the target nuclei, which is seen as the 
bremsstrahlung of a massive photon from the fast going incident quark, q — > qj*. This models 
allows the authors to compute equally the nuclear modifications of parton density in nuclei when 
the momentum fraction of the target parton, X2, is small, X2 ^ 1 (shadowing). Performing a 
global fit of the E772 and E866/NuSea measurements, a huge energy loss has been extracted: 

(66) = 2.7 ± 0.4 ± 0.5 GeV/fm, 

dz 

that is much larger than the above BDMPS-based expectation -dE/dz 1 GeV/fm, but in fact 

this value includes "vacuum energy loss" [339] which should be removed to obtain the medium- 
induced energy loss. It was later shown, however, that the large uncertainty of the sea quark 
nPDF in the E772 and E866/NuSea data makes it difiicult a model-independent extraction of 
energy loss in cold nuclear matter from those measurements [335]. Rather, it was proposed to 
use the lower energy pion-beam data from the NA3 collaboration which is sensitive to valence 
quark densities (while anti-quarks from the sea are probed in p+A collisions) for which nuclear 
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Fig. 49. - Compilation of tlie different estimates for the magnitude of an incoming quark mean energy loss 
per unit length, (-d£/dz)i„, in a L = 5 fm nucleus. Taken from Ref. [319]. 



modifications are known to be small. From the statistical analysis of these data, a much smaller 
quark energy loss in nuclei was obtained [335] 

(67) -— = 0.20 + 0.15 GeV/fm, 

dz 

equivalent to § = 2{dE / dz) / L x 0.14 ± 0.11 GeV/fm [319], and in good agreement with the 
BDMPS estimate. 

Various estimates on cold nuclear energy loss extracted phenomenologically from a variety 
of observables, are summarised in Fig. 49 and Table IX. Most of them point out to a rather small 
energy loss, dE/dz < 0.5 GeV/fm. In order to clarify this issue, Garvey and Peng proposed [318] 
to measure the nuclear dependence of DY production in fixed-target p+A collisions at low beam 
energies, Eiah - 50, 120 GeV (respectively y/s^N - 9.8, 15 GeV), where the effects of energy 
loss prove the strongest. 

6'3.4. Hadron-nucleus and nucleus-nucleus collisions . 
The discussion of jet quenching in A + A is typically focused on parton radiation in the hot 
medium created in the collision, and neglects energy losses in the nuclear target. Similarly, in 
h + A collisions, where no hot and spatially extended medium is created, energy loss effects are 
totally disregarded. However, as discussed in Section 6 3, the cold QCD matter in the target 
nucleus can induce sizeable IS or FS parton energy loss in DY processes and e + A collisions. 
Here we discuss cold matter energy loss effects on hadron production in /; + A and A + A, and 
review the phenomenological studies performed in Refs. [43, 340]. 

The initial state parton suffers multiple scatterings and medium-induced gluon radiation. In 
a simple phenomenological model [340], the resulting energy loss may be accounted for by a 
shift of the incoming parton fractional momentum, xi — > xi(l - e), with e - a:A'^^ the fractional 
IS energy loss. The effect of such energy loss is felt in a kinematic region where the flux of 
incoming partons varies rapidly with x\, typically at large rapidity. In [340], the suppression 
factor RdAuipr) for charged hadrons was computed including coherent parton multiple scatterings 
(higher twist shadowing) and energy loss. A value e - 0.0175 was fitted to STAR data on RdAu 
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Fig. 50. - Effects of IS energy loss and shadowing (in the "high-twist" formalism) for hadron production in 
d + Aii collisions: RdAuipr) for forward pions measured in STAR at -sfsNN = 200 GeV (left) and pion NA35 
data at various rapidities at -yTsjviv = 19.4 GeV (right). Plots taken from Ref. [340]. 



for charged hadrons at /; = 4.1 and -^siviv - 200 GeV, see Fig. 50. Applying the same formalism 
to CERN NA35 data at yflfm - 19.4 GeV, one obtains a reasonable description of the rapidity 
dependence of RdAn, and sees that IS state energy loss becomes relevant only at forward rapidity 
y-ycm ^ 0. According to the rapidity shifts listed in Table IV, we may expect a similar conclusion 
to hold for y - y„„ > 2(5) at RHIC (LHC). 

As discussed in Ref. [43] and Section 2'2, a parton scattered at negative rapidity, y - y^m < 0, 
in a h + A collision travels in the same direction as the target nucleus: seen in the nucleus rest 
frame, it appears to move slowly and corresponds to a low value of v in the language of ^ + A 
collisions. Therefore, based on the observed suppression of hadron production in lepton-nucleus 
DIS [21,22,22,45,71] at low v (see the kinematic analogy between DIS and h + h collisions 
at LO in Section 2 2), we can expect non-negligible hadron suppression due to FS interactions 
in cold nuclear matter also in h + A and A + A collisions. Using Eq. (21), it is possible to plot 
the experimental nDlS data on from HERMES and EMC in terms of the h + h kinematic 
variables pj, y\ and z, which provides a rough estimate of final-state hadron attenuation inh+A 
collisions [43]. The obtained quenching is not small, and increases with decreasing rapidity 
y\ as expected from the kinematic analysis of Section 2'2. However, the pj range covered by 
HERMES and EMC is quite limited compared to the pj for which hadron production inh+A 
and A + A can be measured. Moreover the value of z ^ Zh is not experimentally accessible 
in hadronic collisions and not easily correlated to the measured variables [43]. Furthermore, 
the A-dependence of hadron quenching in nDlS is non trivial, and very different from a naive 
A" power law as often assumed, see Section. 8' 1.1. For these reasons, a theoretical estimate is 
needed for h+ A and A + A collisions. This can be obtained in the energy loss model described 
in Section 63.1 [43], and is displayed in Fig. 51, which quantifies final-state "cold" hadron 
quenching by the ratio 



(68) R'l/pT,y) - 



' A dpjdyidy2 



1 dcr 



pB->hX 



B dpjdyidy2 



The plots show a substantial final-state hadron quenching akeady for midrapidity hadrons at SPS 
and FNAL energy. At RHIC it is sizeable for y/, < -2, where it may play a role in understanding 
the evolution of the Cronin effect at backward rapidity, is still present at y/, = 0, where it is of 
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Fig. 51.- Energy loss model estimates of final-state hadron quenching in cold nuclear matter for midrapidity 
pions produced in p + A collisions at SPS and FNAL fixed-target energies, and several negative rapidities 
at RHIC. 



order 5% at pr > 10 GeV/c consistent with nuclear PDF modifications, and quickly disappears 
at forward rapidity. In p + Pb at the LHC, we may expect negligible final-state effects at yi, > -3 
because of the rapidity shift Ay in Table IV. The rapidity regions in rapidity (measured in the 
centre-of-mass frame) where we can expect IS and FS hadron quenching in li + A collisions are 
summarised in Table X. 

In A + A collisions the scattered partons can traverse two distinct QCD media: first the two 
colliding nuclei, then the hot medium formed after the collision. The two media are well sepa- 
rated in time since the time-scale for QGP formation, ^(1 fm) is (much) larger than the crossing 
time of the colliding nuclei: Tcros, = 2RA/y ~ 1.5, 0.15, 5 10"^ fm at SPS, RHIC and LHC re- 
spectively. The degree to which the parton interacts with the 2 cold nuclei depends on rapidity: at 
mid-rapidity to the same degree; at large rapidity more strongly with the comoving nucleus and 
more weakly with the opposite moving nucleus. A rough estimate of the energy loss in the cold 
nuclei can be obtained by multiplying the R'j-^ values in Fig. 51 at 77 and -rj. At RHIC midra- 
pidity, cold quenching is less than 10% at \y\ < 1, much smaller than the observed factor 4-5 
hadron quenching observed in central Au + Au collisions. At LHC, due to the even larger longi- 
tudinal boost, cold nuclear matter effects will be smaller than 10% in the \y\ < 4 rapidity range. 
At SPS energy, cold quenching at midrapidity is of order 50% and may in fact be comparable 
to hot quenching: they both need to be taken into account in any QCD tomographic analysis to 
characterise the properties of the produced matter. 

6'4. Jet quenching from AdS/CFT duality . - 

The energy loss computations discussed so far are based on specific models of the medium 
rooted in perturbative QCD. On the other hand, indications come from RHIC data that the 
medium produced in heavy-ion collisions is strongly coupled, so that non-perturbative effects can 
become important [341]. However, non perturbative methods like lattice QCD have an intrinsic 
difficulty in computing dynamical properties, which are connected to the Minkowski geometry 
of space-time, and not easily simulated in Euclidean space-time as required by most lattice QCD 
methods. 

A recent theoretical development, the Anti-de Sitter / Conformal Field Theory (AdS/CFT) 
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Table X. - Regions of rapidity where IS and FS effects in cold nuclear matter may play a role in quenching 
hadrons in p + A collisions at various energies. In A + A collisions the effects should be symmetrised in 
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duality [342-345], has given access to analytical calculations of the dynamical properties of a 
plasma in ^ = 4 supersymmetric SU(A'c) Yang-Mills (SYM) theory - rather than QCD - in the 
limit of strong coupling and large number of colours Nc, by relating them to a weakly-coupled 
string theory living in a 10 dimensional space. The AdS/CFT duality has been later extended to 
a large class of supersymmetric theories, and is also known under the name of gauge/string du- 
ality. Even though QCD and supersymmetric Yang-Mills theories are quite different in terms of 
running coupling and matter contents, recent experimental and theoretical considerations suggest 
that the many-body physics of the QGP near the phase transition, where experimental findings at 
RHIC suggest that the coupling is strong, is similar in the two theories [346]. 

More in detail, the computation of an observable in the quantum gauge theory at strong cou- 
pUng can be rephrased as a classical computation of a related observable in a higher-dimensional 
gravity theory at weak coupling, where it is perturbatively calculable: 

^ = 4 SYM <^ type IIB string on AdSs x S ^ 

Here, ^ - A means that the theory contains 4 supercharges, AdSs is a 5-dimensional space 
with constant and negative curvature, 5^ is a 5-dimensional sphere. For finite temperature gauge 
theories, relevant to the description of a QGP, the gravitational equivalent contains a special kind 
of black hole: 

J/ =A SYM at finite T <r-y type IIB string on (Schwarzchild-AdSs black hole) x S ^ 

The details of the computation of dynamical quantities such as the drag force coefiicient of a 
moving heavy-quark in the plasma, the quenching parameter q, and the jet transverse momentum 
broadening in the AdS/CFT duality framework are reviewed and summarised in [346-348]. Al- 
though stiU highly speculative, the connection of these theoretically computable strong-coupling 
quantities with the physics of real-life QGP is opening new ways of understanding the deconfined 
phase of QCD near the critical temperature at a fundamental level. 

7. - Hadron formation, propagation and interaction 

In hadron absorption models, hadronisation is typically assumed to happen in two stages 
as shown in Fig. 52. One effectively identifies the prehadron production time and the colour 
neutralisation time: (i) the struck quark neutralises its colour and forms a so-called "prehadron", 
which then (ii) collapses into the asymptotic hadron h (see Fig. 7 and Section 2'4). When the 
(pre)hadron is formed inside the nucleus, it can reinteract with the surrounding nucleons. The 
space-time evolution of the hadronising system can be computed using the non-perturbative Lund 
string model or in a pQCD-based approach. The various realisations of the two-stage model difl'er 
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Fig. 52. - Left: Sketch of the time evolution of the hadronising system in the two-step hadronisation model: 
a quark q is struck at point y which evolves into a prehadron and hadron at point / and /' respectively (in 
the Lund string model, /' = / + zhv/k). Right: Sketch of the effective cross section (70). 



by the assumptions concerning the prehadron production time tpreh, the hadron formation time f/j, 
the (pre)hadron interaction with QCD matter, and by the treatment or neglect of the subleading 
effects due to the quark propagation stage. Most models address absorption in cold nuclear 
matter, especially in ^ + A collisions. We will mostly specialise the discussion to this case, and 
briefly comment about applications toh+A and A+A collisions. 

71. Early string-based absorption models. - 

The first nuclear absorption models [349-351] assumed a one-step hadronisation process, in 
which the struck quark propagates in the nucleus and interacts with the surrounding nucleons 
with cross section cr^. After a time (f/,) oc oc v the hadron is formed and interacts with cross 
section cr/j. However, the EMC collaboration showed in [45] that for any choice of the parameters 
(Tq and (th) only a poor description of charged hadron data infi+A collisions could be achieved. 

The idea of an intermediate "prehadronic" stage stage between the quark and the hadron was 
proposed in Ref. [63] in the context of the Lund string model. Hadronisation was imagined to 
proceed first through the breaking of the colour string stretched between the struck quark and 
the nucleus, then through the evolution of the string pieces, whose end-point quarks take some 
time to come together and finally recombine to form the final hadron (Fig. 53), see Section 2 4 
for more details. Typically, it is assumed that the string breaking process does not depend on 
the nature of the target. The nuclear dependence comes from the interaction of the hadronising 
system with the surrounding nuclear medium. In principle, the hadronising system is allowed 
to interact at all stages with cross-sections o-g, o-preh, cth (Fig. 52, left). In the case of { + A 
collisions, the hadron survival probability Sa, i e., the probability that the produced (pre)hadron 
does not interact with the nucleus, can be approximated as [63] 



where pa is the nuclear density normalised to 1, and D(zh',y' - y) is the probability for a string 
breaking (i.e. prehadron formation) at a distance Ic = y' -y, called "constituent length", from the 



(69) 




X I dy' D(zh;y' - ym - I dya-e0,y')pA{b,y) 
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Fig. 53. - Left: Sketch of the independent Lund fragmentation model. Cy and F, are the constituent (pre- 
hadron) formation point and "yo-yo" (hadron) formation point, respectively. Right: Sketch of the string-flip 
model in cold QCD matter. The 5, are diquark configurations initially connected to the quark q/ in one of 
the nucleons. Colour exchange interactions rearrange the string end-points in several short stringlets and a 
leading one with reduced energy. Figures taken from [352]. 



photon-quark interaction point. In the Lund string model, the hadron is then formed at a distance 
ly = y' +Zhy/K-y, also called "yo-yo" length. The two-step dynamics is contained in the effective 
cross section 



depicted in Fig. 52 right. Note that the ansatz (69)-(70) neglects elastic scatterings of the system, 
which induce an energy loss and feed down to lower Zh, so that its vaUdity is confined to Zh ^ 0.4. 
This limitation can be overcome in Monte Carlo implementations, see Ref. [352] and Section 7'4. 

The early applications of the independent string fragmentation picture, [63, 353, 354] focused 
on the relevance of the constituent vs. yo-yo lengths (quark-to-prehadron vs. prehadron-to- 
hadron formation times) in describing hadron production data m p + A and ( + A collisions 
assuming CTpreh = o"/, or cr^ = CTpreh- The key point is that at large Zh, ly = h °^ ZhV, while 
oc (1 - Zh)v. EMC data on leptoproduction of hadrons confirmed that it is in fact the constituent 
length Ic which controls hadron attenuation on nuclear targets, with a negligible contribution 
from the struck quark interactions: data could be well fitted with cr^ < 0.75 mb and 0-^,,^/, - cTh ~ 
20 mb [45], see Figure 54. In other words, the hadronising system starts interacting inelastically 
with QCD matter already at the prehadronic stage, well before the final hadron is formed. 

The effect of the nuclear medium on the string breaking process has been considered in 
Refs. [352, 355]. The struck quark qo (see Fig. 53 right) can interact with other nucleons with a 
cross-section of the order of the hadronic one because of the non-zero transverse size of the string 
attached to it [355]. When it does so, there is a probability of colour exchange with the nucleon, 
so that the qo reconnects with the diquark q\, and the quark q\ with the diquark qo originally 
attached to q^, and so on. This process was called a "string-flip" [352], and leads to the creation 
of a number of different string configurations: 



Assuming a color-exchange cross-section cTce, the mean free path is Xce = (CccPa) . with pA ~ 



(70) (reff(y,y') = (Tqe{y'-y) + (Tpreh0(y- y') 0[y' 




yj + o-he{y-y' , 



7* A (qoqi) + (qiqi) + ... + (qn-iqn) + (qnqo) + (A - n - 1) 
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Fig. 54. - Comparison of EMC R'|^J data to thie one- and two-step models with varying quark-nucleus and 
prehadron-nucleus cross sections. Adapted from Ref. [45]. 



0.17 fm . The independent string fragmentation model is obtained in the limit Ace ^ Ra- In 
the opposite limit, A^e <K Ra, the struck quark qo exits the medium reconnected to q„ (leading 
string) close to the nucleus surface. Therefore all the hadron components are produced outside 
the nucleus and do not suffer nuclear absorption. The intermediate stringlets have typically a low 
mass M ^ yllKm^Ace ~ 2 GeV/c^, decay in a few low-energy hadrons, and do not contribute 
at large Zh- Suppression of large-z/, hadron production is then due to the energy lost by the 
leading string because of colour reconnections, and is independent of the details of the hadron or 
prehadron formation time. Defining the probability for a leading hadron of energy Eh = ZhV to 
fragment out of a string with average energy-loss e - AceK, and assuming it to be independent of 
the energy of the string, one obtains a medium modified fragmentation function 



For fragmentation functions decreasing with Zh, this leads to hadron attenuation: Da < D. 
Eq. (71) is a close analogue of the z;,-shifted medium modified fragmentation function, see e.g. 
Eq. (45), used in the gluon bremsstrahlung energy-loss models considered in Section 6. The 
hadron attenuation caused by colour exchange can be naively interpreted as quark absorption in 
the nucleus 



However, the effective quark cross section cTq is process-dependent. For lepto-production of light 
flavours it is rather small, in agreement with the phenomenological EMC fit [355]. 

Applications of the string-flip model to EMC data have been considered in Monte Carlo 
[352] and analytic [356] implementations in [352], a rather large cross-section cr^.^ - 30 mb 
was considered, assuming a string to have a typical hadronic transverse size, so that Ace » ^a- 
Therefore, the leading string is produced at the surface of the nucleus, and hadron attenuation 
is purely due to the energy loss induced by colour exchanges, irrespective of the prehadron 



(71) 




(72) 
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production time. EMC hadron attenuation could be well reproduced both as a function of v and 
of Zh- It leads to a slightly larger hadron suppression at v < 20 GeV and Zh ^ 0.5 compared to 
the independent fragmentation model using the constituent length, otherwise giving comparable 
results. In Ref. [356], the colour exchange cross section was identified with the constituent-quark 
cross section in the additive quark model: (Tee ~ crq = O-^ctjtn ~ 10 mb. The smaller cTce leads to 
an intermediate model between the independent string fragmentation and pure string-flip model. 
However, computations turn out to underestimate the suppression in the EMC kinematics. 

Finally, we should mention that extensions of string models to address double hadron attenu- 
ation have also been discussed [78, 357]. 

7'2. Modem string-based absorption models. - 

Modem string-based absorption models build upon the early models discussed in the pre- 
vious Section, exploring different possibilities for the effective cross section tr,^ and/or for the 
hadron survival probability (69). Their computations have also been compared to the more recent 
HERMES data on e" + A colUsions at = 12 - 27 GeV [21-23,44]. 

The AGMP model [61,64] is based on the Lund string model estimate of the formation 
times discussed in Section 2'4, and neglects the interactions of the struck quark (cr^ = mb) 
in agreement with fits to EMC data [45]. In the HERMES kinematics one typically finds a 
production time {tpreh) ~ 4 fm < and a hadron formation time (f/,) ^ 6 - 10 fm > Ra'. 
the hadron is typically formed at the periphery or outside the nucleus so that its interaction 
with the medium is negUgible. On the contrary, after its formation, the prehadron is allowed to 
interact with the surrounding nucleons with a cross section crpreh(y) = 0.80 cr/,(v) proportional 
to the experimental hadron-nucleon cross section cr/, (taking e.g. production data on a Kr 
target [23]). The probabiUty S^^iijz, v) that the (pre)hadron - produced from the fragmentation of 
a quark scattered at point {b, y) - does not interact, can be computed using transport differential 
equations [61]: 



where Al - k - Ipreh, and pa is the nuclear density. The hadron multiplicity is computed, at 
leading order in pQCD, as 



Here a-(f and o-(a are the lepton-quark and lepton-nucleus cross sections, ^/ is the /-quark 
distribution function, and its fragmentation function. The model is applicable for 0.4 < z < 
0.9, where it describes EMC [64] and HERMES [61,64] experimental data on a wide range 
of hadron flavours and targets, see Fig. 55. As it stands, the pure-absorption AGMP model 
does not predict a dependence of the attenuation ratio Rm on 0^, because neither the prehadron 
production time tpreh derived in the Lund model, nor its assumed cross section CTpreh depend on 
it. This assumption does not contradict the slightly rising RuiQ^) measured at HERMES [23], 
see Section 3'2. 
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Fig. 55. - Lund string model based absorption model [64, 321] (solid) and energy-loss model [43, 321] 
(dashed) compared to HERMES data on n* production [23]. The prehadron cross section is iTprrf, = O.Sct;, 
and the transport coefficient is q = 0.6 GeV-/fm. Experimental statistical and systematic errors have been 
added in quadrature. 



A different variation of the two-scale model has been proposed in Refs. [358, 359]. The 
authors argue that in 3-h1 dimensions, the yo-yo time is ill-defined because the end quarks of a 
string snippet never really meet together, and propose instead that the hadron is formed after a 
proper time for quark recombination, t^, following string breaking at a time f taken from the 
PYTHIA generator. Then, the hadron formation time is simply set to f/, - (tpreh + Tr)vstr7str, 
and boosted according to the string velocity v.,/, and Lorentz factor jstr- The system does not 
interact in the nucleus until hadron formation, at which time it starts interacting with the full 
hadron-nucleus cross section analogously to the AGMP approach. The model does a fair job in 
describing pion and kaon attenuation on and Kr [21,22] at Zh S 0.4, with a flavour-dependent 
recombination proper time Trin) - 0.8 fm and Tr(K^) - 0.4 fm. It is interesting to note that 
the fitted values are such that t, (/z) « l/m/,. However, this model overestimates the recently 
measured suppression on Xe targets [23], and no comparison to EMC data is provided. 

Extensions of the two-scale model based on Eqs. (69)-(70) have been considered in Refs. [78, 
357], and more recently compared to double hadron attenuation at HERMES [24] in Refs. [94, 
95,360]. 

7'3. The colour dipole model. - 

In Ref. [29] the formation of a leading hadron (zh S 0.5) is described in a pQCD-inspired 
approach based on Refs. [66-68], see also Section 2 4: a hard gluon radiated at the interaction 
point splits into a qq pair, and the q recombines with the struck q to form the leading prehadron, 
which evolves into the observed hadron. The time development of hadronisation is included in 
this picture by observing that the radiated gluon can be physically distinguished from the struck 
quark only after a coherence time tc = 2Ea(l - a)/k^ , where kj and a = Eg/v are the transverse 
momentum and fractional energy of the emitted gluon. If the leading hadron is produced at large 
Zh and contains the struck quark, none of the radiated gluons can be emitted with a > 1 - Zh 
by energy conservation. The time-dependence of the energy radiated into the emitted gluons is 
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computed as 
(75) 



AE(t,Zh,Q) 



1 dNg 

dzav @(t-tc) 
-zf, dki,dz 



where dNg / dk^dz - aAkj) I (37r) 1 / {zk^) is the Gunion-Bertsch spectrum of radiated gluons [69] . 
The upper limit is imposed by the fact that gluons with kj > Q should be considered part of the 
struclc quarlc [11]. Next, the gluon (with momentum kj) splits into a qq pair. In the large-A'c 
approximation, the antiquark and the struck quark form a colourless dipole, which is identified 
as a prehadron. The prehadron production time is identified with the coherence time of the 
gluon (rather than with the qq splitting time) and hadron formation is computed by the overlap 
of the qq dipole with the hadron light-cone wave function Assuming that the q and the q 
in the pair share the same amount of gluon energy and transverse momentum, one can compute 
the probability distribution W{t, Zh,Q^, v) that the prehadron is formed at a time t after the y*q 
interaction: 



, f da r^' dkl \ I „\ 

Jo a Ja2 kl [ \ 21 



(76) 





exp(~t/tc) 


V 


tc 



a 3 



Qitp[-Ng{zh,t,Q^,v)) , 



where, £^,(0 = v - AE(t), T*/, is the hadron light-cone hadron wave function, .^V^ is the number of 
gluons radiated within a time t, and is a normalisation factor. Numerical results are presented 
in Fig. IL Integrating W over t one obtains the vacuum fragmentation function for the leading 
hadron which compares favourably with global fit FF from Refs. [324, 325] in the range Zh = 
0.5 - 0.9 and = 2 - 10 GeV^. 

The qq dipole, which is assumed to be formed with a Gaussian transverse size around an 
average (R^) oc propagates through the nucleus and fluctuates in size. According to colour 

transparency, see Sect. 7'5, it interacts with the nucleus with a cross section o-qq = C(Eh)r^. 
All efl'ects of fluctuations and nuclear interactions are computed in a path-integral formalism for 
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dipole propagation in QCD matter [361,362]. Finally, the effective in-medium fragmentation 
function reads 

db" dypAily) dtW(t,Zh,Q\v)Tr(b,y + t) 

J-oo Jo 

where Tr(b,y + t,oo), called nuclear transparency in [29], is the probability for a prehadron 
not to be absorbed in the nucleus. The hadron attenuation ratio can then be approximated as 

Finally, medium-induced energy loss is included in the model as an increase in the parton 
energy loss proportional to the quark transverse momentum broadening Apj, according to the 
relation AEind = |ajA/?^L, derived in Baier et al. (BDMPS) [140]. (Note that this formula is 
valid for an asymptotically large medium, which is not the case in e -h A collisions; the authors 
of [29] argue that finite medium size corrections are small at HERMES energy, but their effect 
might in fact be non-negligible, see Fig. 46.) The induced energy loss modifies the production 
time distribution via 

AE ^ AE + A^ind 
(78) , , , 

Q^^Q^ + Apl, 

with Ap^ computed in the colour dipole formalism [141, 142]; Apj- = 2C(Eq)pAL. This leads to 
an additional, but subleading, hadron suppression compared to the effect of prehadron absorption, 
as shown by the two curves in Fig. 56 left. 

The parameters of the model are fitted to other processes than hadron production in / -i- A 
collisions and, in this sense, the approach can be considered parameter-free. The comparison 
to HERMES experimental data for ;r* production is shown in Fig. 56. The model can also 
describe the smaller suppression compared to ;r* (not shown here), as well as the EMC data 
at zi, > 0.5. Use of the colour dipole formalism allows also the computation of the Cronin effect, 
shown against data in Fig. 56 right. The dependence of Rm is discussed in Sect. 7 5 (see 
Fig. 61). An interesting consequence of this formalism is that the prehadron production time is 
inversely proportional to Q^: 



(79) (tpreh) oc (1 - Zh)^ . 



where v « Eg. Technically, this arises from the upper limit of integration in Eq. (76). Physically, 
it is interpreted saying that a quark which is struck by a photon of large virtuality radiates more 
intensely than for a lower virtuaUty: as a consequence, it will be able to travel only a shorter 
distance before hadronising at a given Zh- In p + A and A + A collisions, oc p^^ and for 
midrapidity hadrons Eq ^ prh, hence {tpreh) ^ ^Ipt- So, the prehadron production time actually 
shrinks with increasing pT, instead of increasing because of the large Lorentz boost. This would 
imply that prehadrons are formed well inside the hot medium produced in A -i- A collisions at any 
Pj. Experimental confirmation of Eq. (79) is thus very important. It can be accomplished in a 
direct way by measuring the hadron pr-broadening, or indirectly with a scaling analysis of Rm, 
see Sect. 81. 



7'4. The GiBUU transport Monte Carlo. - 
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Fig. 57. - Sketch of the time evolution of the (pre)hadronic cross section for arbitrarily chosen formation 
times. Solid lines are for constant crp„;,, dotted and dashed lines are Eq. (84) for n = 1, 2. A typical nuclear 
distance of 7 fm is indicated as a vertical dashed line. Figure taken from Ref. [363]. 



In Refs. [363, 364] a Monte Carlo event generator has been used to explore hadron attenuation 
in ^ + A collisions in the framework of a transport model based on the Boltzmann-Uehling- 
Uhlenbeck (BUU) equation. In this approach, the lepton-nucleus interaction is split into two 
parts: 

1 . The exchanged virtual photon produces a final state X which is determined using the Monte 
Carlo generators PYTHIA and FRITIOF based on the Lund string model. Nuclear effects 
like binding energies, Fermi motion, Pauli blocking and coherence length effects that lead 
to nuclear shadowing are also implemented. 

2. The state X is propagated through the nuclear target using a BUU transport model [364, 
365] which - through a probabilistic coupled-channel computation - accounts for particle 
creation, annihilation and elastic scattering in the final-state interactions. 

The latest implementation of the code is known as GiBUU [363]. The space-time evolution 
of hadronisation has been explored using two concepts of prehadron differing mainly in the 
treatment of the production time tpreh'- set to zero in transport models [366-368], and depending 
on the energy and momentum of the string fragments in the Lund model, see Sections 2 4 and 7' 1 . 

In Ref. [364], the "transport model" view is adopted. The string decay into colour neutral 
prehadrons is assumed to happen instantaneously, hence tpreji - fm. Hadrons are assumed to 
form after a fixed proper formation time t/ = 0.5 fm in the hadron rest frame, which is then 
boosted to the laboratory frame. 



(80) 



ZhV _ 

m 



Between the prehadron and hadron formation times, only beam and target remnants, i.e., hadrons 
containing the valence quark of the struck nucleon or of the resolved photon, are allowed to in- 
teract with the rest of the nucleus. Specifically, the prehadron cross section is assigned according 
to the constituent quark model: 



(81) 
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Fig. 58. - Multiplicity ratio in the GiBUU model as a function of pj. Computations are done with parton 
intrinsic momentum <4> = 0-44 GeV^ (left) and <4> = 1 GeV^ (right). HERMES data are from [22], Plot 
taken from [363]. 



where «org denotes the number of quarks or antiquarks in the hadron coming from the beam or 
target nucleon. This model can describe fairly well HERMES data for tt'^, n'^ and K"^ at large 
Zh <; 0.4. For smaller Zh some discrepancy for production on light nuclei emerges, related 
in part to the assignment of the pre-kaon cross section (81), in part to feed-down due to the 
decay of diffractively produced p mesons into K^K , or due to ttN interactions. The results 
for baryons are less good, especially for p production, which cannot explain the strong rise 
of experimental R'^ at zi, < 0.5. However, the biggest challenge for this space-time scenario 
is in connection to EMC data, whose suppression is vastly overestimated due to the unlikely 
assumption of an instantaneous conversion of string fragments into colour neutral prehadrons 
and of an instantaneous jump from the prehadron to the hadron cross-sections at time f/,. 

The solution proposed in Ref. [364] and fully worked out in Ref. [363] is to use the concept 
of prehadron coming from the Lund model, and time-dependent prehadron cross sections. The 
PYTHIA generator is used as main source of information on the string fragmentation process. 
In the Lund string model, we can actually associate three time scales to each hadron: the times 
tpi and tp2 at which the constituents of the hadron are created, and the time ?/, at which the 
constituents meet and form the hadron. The prehadron production time is then defined as 

(82) tpreh = min(fpi,fp2) ■ 

This definition is different from what is commonly used, see Sections 71 and 7'2. Most no- 
tably, the rank-1 hadron contains the struck quark, which preexists the 'y*N interaction: hence, 
' = 0. The determination of the prehadron and hadron formation time is caiTied out on an 
event-by-event basis by extracting the coiTesponding information from the PYTHIA generator as 
described in detail in [65]. An interesting result is that the prehadron and hadron formation times 
are quite insensitive to the hadron mass, contrary to the common assumption that they scale with 
a Lorentz boost factor oc 1/m/,. The prehadron cross section includes the main features of colour 
transparency, and is allowed to vary in time according to linear and quadratic evolution, or to a 
quantum diffusion picture, see Eqs. (84)-(85). The requirement of reproducing both HERMES 
and EMC data rules out the quadratic evolution, and slightly favours the quantum diffusion pic- 
ture with (Tpi-gi, oc at f = tpi-eh, instead of crp„/, = 0, see Figs. 59 and 60 left. The model 
performs well in describing the flavour dependence of Rm at HERMES, with the exception of a 
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Fig. 59. - Multiplicity ratio in the GiBUU model, compared to HERMES at £<. = 27 GeV [22] and EMC at 
Eg = 100 - 280 GeV [45] data. The upper and lower bounds of the shaded band corresponds to higher and 
lower EMC beam energy, respectively. The prehadron cross section evolution is constant (83), linear and 
quadratic in time (84), from left to right. Plot taken from [363]. 



slight underestimate of K~ attenuation over the whole Zh range, and the lack of a strong rise in the 
proton sector at Zh < 0.4. It can also describe well the Cronin effect on hadron pr distributions, 
without need for a broadening of the intrinsic parton transverse momentum (fcj-) as commonly 
assumed, see Fig. 58. 

An important lesson coming from these simulations is that the geometrical acceptance of the 
experiment does not fully cancel out in the double ratio which defines the hadron attenuation ratio 
R'Ij Eq. (28), especially at low Zh 5 0.4. This is shown in Fig. 60 right, where unidentified charged 
hadron attenuation on Kr at HERMES [22] is computed with a 4;7r detector (dashed line) and the 
HERMES detector geometric acceptance (solid line), which is important at Zh £ 0.4 [364]. In 
both computations the full kinematic cuts have been included; the dotted Une shows the effect of 
neglecting the £/, > 1 .4 GeV cut. 

7'5. evolution and time-dependence of the prehadron cross section. — 

In the context of string models, hadron attenuation is typically constant in Q^. However, a 
dependence of the prehadron cross section on is to be expected on general grounds due to 
colour transparency (see e.g. [369]), which is a direct consequence of the composite nature of 
hadrons being made of quarks so that hadrons fluctuate between quark configurations of different 
sizes. In [370-373] it was found that in scattering of colour neutral objects, configurations with 
small transverse size h have cross section oc \p-. The size of a colour neutral prehadron is likely to 
be determined by the hard scale of the process, so that we may expect, at least for the leading 
prehadron, cTpre/, oc XjQ^. 

After its production, the prehadron must evolve in size until hadron formation; correspond- 
ingly, its cross section will evolve in time from CTpreh to cr/,. Let us picture the prehadron as a 
qq colour dipole for simplicity. In Ref. [11], it was shown that in a classical picture the size of 

a pair of colour charges with average -^QB^) evolves linearly in time. Therefore, due to colour 

transparency, its cross section evolves quadratically: 0"^^;, oc t^ . If one includes the uncertainty 
principle in these considerations, the assumption of fixed transverse momentum is not valid, 
since as r — > we have kj — > oo. As a consequence, the rise of the cross section is softened 
and becomes linear in time: CTpreh o^ t. On the other hand, the struck quark is not a bare particle, 
but is still accompanied by those parts of the original nucleon field with Fourier components 
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Fig. 60. - Left: As in Fig. 59, but for the quantum diffusion scenario (85) compared to the hnear time- 
dependence scenario (dashed line). Plot taken from [363]. Right: Effect of the HERMES detector geometric 
acceptance on measurements of charged hadron attenuation. Solid line: full geometric and kinematic cuts. 



Dashed line: An detector with full kinematic cuts. 
Figure taken from Ref. [364]. 



Dotted line: An detector with > 1 .4 GeV cut neglected. 



satisfying > Q^. Hence, the time dependence of a-preh may be expected to lie between linear 
and quadratic. 

In the context of semi-inclusive hadron production in ^ 4- A collisions, this idea has been 
phenomenologically explored in Refs. [363, 364] using the GiBUU Monte Carlo generator (see 
Section 7 4 for more details) and in [360] in the context of the two-scale string model, see Sec- 
tion 7 2. The time evolution of the prehadron cross section is also implemented by means of 
the colour dipole model in the perturbative fragmentation model of [29], which also includes a 
2^-dependent production time tpreh (1 - Zh)vjQ^, see Section 7'3. 

In Ref. [363], four scenarios have been explored. In the first, the prehadron cross section is 
kept fixed, 

(83) a-preh - -CTh , 

where the proportionality factor has been chosen to reproduce HERMES data. (Note that differ- 
ent proportionality factors arise in different models due to different assumptions on the under- 
lying space-time development of hadronisation, e.g., <jpreh - 0.67cri, in [61], where prehadron 
interactions are treated differently.) In the next two scenarios, 

(84) (Tprehit) = « = 1,2 

corresponding to the quantum mechanical linear time dependence and the classical quadratic 
time dependence discussed in [1 1]. Note that at f = fp„/,, the prehadron is created with zero cross 
section, and that no dependence on is included. The last scenario implements the quantum 
diffusion picture of Ref. [374], which combines the linear rise in time with a non-zero and Q^- 
dependent value for the initial prehadron cross section: 

/ t - tpreh \ IGeV^ 

(85) O-prehit) = Cro -H (1 - CTq) CTh CTq = riead , 

\ f/i - tpreh I Q 
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Fig. 61. - Left and bottom right: dependence of hadron attenuation in the colour dipole model of 
Ref. [29] for Ne and Kr, respectively. Solid lines include medium-induced gluon bremsstrahlung, dashed 
lines do not. Points are preliminary HERMES data [23]. The two plots are taken from [29]. Top right: 
dependence of hadron attenuation in the absorption model of Ref. [360]. Solid lines computed in the 
improved two-scale model, with time evolution of the prehadron cross section. Dashed lines computed in 
the two-scale model without time evolution of the prehadron cross section. Plot adapted from [360]. 



where riead is the number of the struck nucleon valence quarks contained in the produced hadron 
divided by its number of valence quarks. E.g, a leading pion has riead = 1/2 and a subleading one 
has riead - 0, see Section 7'4. In all scenarios, the prehadronic cross section is before t = tpreh 
and reaches the hadronic cr/, at f = f/,, and includes the basic features of colour transparency. 
These cross sections are sketched in Figure 57. While HERMES data alone do not give enough 
leverage to distinguish between these time evolution scenarios [360, 364], the combined analysis 
of EMC and HERMES data singles out the linear evolution scenario [363], see Fig. 59. However, 
it cannot distinguish between an initial o-prei, = or o-prei, 1 / g^, despite a large difference in 
the average {Q^} in the two cases, even though the latter is slightly preferred, see Fig. 60 left. 
This result is in line with the slow variation of the attenuation ratio Rm with observed at 
HERMES [23]. 

An analysis of the dependence of 7?^ has been performed in [29, 360]. Based on colour 
transparency, one expects Rm = a + b/Q^ + d'i l/ Q*) with negative b. However, if time evolution 
of the prehadron cross section to reach the full hadronic one is sufficiently fast, its depen- 
dence can become negligible. This is shown in Fig. 61, top right, as difference between the solid 
and dashed line, computed with a linear time dependence (Eq. (85)) and no time dependence 
respectively [360]. In [29], where colour transparency and the time evolution of the prehadron 
cross section are naturally implemented in the light-cone dipole formalism, two further elements 
contribute to the dependence of R''^. The first is the dependence of the production time 
tpreh 1 / Q^, which increases the in-medium prehadron path length and consequently its at- 
tenuation with increasing Q^, working in the opposite direction as colour transparency. The net 
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effect is an almost Q^-independent attenuation ratio for light nuclei, which gains a small slope on 
heavy targets or for low y cuts, see Fig. 61. The second is medium-induced gluonbremsstrahlung, 
which induces an additional attenuation at small Q^, but whose effect disappears at large Q^, see 
Fig. 61. A comparison of model computation with HERMES data on distributions in hght 
and medium-heavy nuclei {Ne, Kr, Xe) from Ref. [23], which show a pretty slow variation of 
Rm with Q^, and the forthcoming JLab multi-differential measurements of z and 2^ distributions 
will shed light on this problem. 

In Ref. [61, 64], a different kind of dependence has been explored, namely, the possibiUty 
of a modification of fragmentation functions induced by a partial deconfinement in the nuclear 
wave function [375-377]. It results in a rescaling of the FF: 



(86) D (z„, Q^) ^ D [zh,UQ^)Q^) UQ^) = 




where is the strong coupling constant computed at the scale at which DGLAP evolution 
starts, Ao is the nucleon radius in the vacuum and Aa is the radius in the medium. The rescaling 
factor depends on the amount of nucleon overlap in the nuclear wave functions, and depends 
on because of DGLAP evolution [376, 377]. 

The phenomenological consequence is a multiplicity ratio which shghtly decreases with 
contrary to what is experimentally observed: partial deconfinement does not seem to affect frag- 
mentation functions (even though it might modify the parton distribution functions, and lead to a 
possible explanation of the EMC effect [109]). This might have to do with the fact that hadrons 
are typically formed outside the medium in the HERMES kinematics, so that deconfinement is 
not affecting them. At JLab, or for smaller values of v at HERMES, hadrons will be in part 
formed inside the medium and partial deconfinement may affect their fragmentation. 

In summary, the 2^ dependence of hadron attenuation is very sensitive to the underlying dy- 
namics of the hadronisation process. However, since different physics inputs in theory models 
can have similar effects on the computation of R^^, this observable must be used in conjunc- 
tion with others like the hadron transverse momentum broadening, and with the requirement to 
describe attenuation in both low-energy experiments (CLAS,HERMES) and high-energy exper- 
iments (EMC), to sort out the physics underlying hadron attenuation. 

7'6. Prehadron absorption in hot QCD matter. 

Ill Ref. [378, 379], prehadron absorption in hot nuclear matter as a mechanism for jet quench- 
ing in A + A colhsions, has been considered in the context of a transport model similar to the 
GiBUU model described in Sect. 7'4. The medium formation and evolution is described in the 
Hadron-String Dynamics (HSD) approach [367,380,381]. The prehadron formation and evo- 
lution is the same as in the GiBUU model, with the additional assumption that hadrons are not 
allowed to be formed if the medium energy density at their point in space and time is above the 
critical energy density for a QGP formation, - 1 GeV/fm"'. This requirement mimics the effect 
of a QGP on high-energy hadrons, even though the produced medium in the HSD simulations 
is always a hadron gas. Note, however, that the energy-density cut is not assumed to affect the 
prehadronic states. In addition, the model phenomenologically includes the Cronin enhancement 
at intermediate pxh by a nuclear enhancement of the intrinsic parton transverse momentum. 

Two different scenarios for the space-time evolution of prehadrons have been used in the 
simulations. In the first scenario, the leading prehadron cross section is constant and the hadron 
formation proper time is kept fixed at t = 0.5 - 0.8 fm. The resulting hadron suppression is 
mostly due to the interaction of the leading prehadron with the medium, and is sufficient to 
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Fig. 62. - Suppression factor R^^ as a function of pj = pn,, for midrapidity charged hadrons in 0-5% 
central Au + An collisions at RHIC top energy, -\fsNN = 200 GeV. The solid curve is the hadron quenching 
computed when the leading prehadron cross section rises linearly in time according to Eq. (84) with tprei, = 
and n = \. Experimental data are from Refs. [149, 156]. The figure is taken from Ref. [378]. 



explain most of the hadron quenching at RHIC. This is interpreted by the authors as evidence 
that prehadronic final-state interactions can result in the observed jet quenching results. However, 
when the leading prehadron cross section grows linearly with time, as requested by a consistent 
description of HERMES and EMC data for nuclear DIS [363], the resulting hadron quenching is 
clearly insufficient to explain data for RHIC central Au + Au collisions, see Fig. 62. Therefore, 
prehadronic final state interactions seem ruled out as main mechanism for hadron quenching at 
RHIC, while they may in fact be responsible for hadron quenching me + A collisions. 

The reason may be that, contrary to the assumptions in the HSD model, prehadrons cannot 
be formed in the QGP phase (e > e^) because of the screening of the colour force needed to 
bind them. If this is the case, the struck quark remains bare and interacts with the medium un- 
til it escapes from it or until the medium enters the hadron gas phase. Only at this point the 
hadronisation process can start but at this stage the medium is too dilute and the prehadronic in- 
teractions are too weak to substantially contribute. On the other hand, in the QGP phase, partonic 
energy loss can be important because of the large colour charge density. If this explanation is 
correct, it would be interesting to confront the HSD prehadron absorption model with Rphphipr) 
data at SPS energies (note, however, that we currently lack high-p^ p + p spectra at around 
yfsNN = 17 GeV [159]) where the QGP is formed in a more restricted range of centralities, if at 
all, and is shorter lived than at RHIC. A good description of peripheral collisions data may val- 
idate the model, so that a substantial overestimate of central collision data might be interpreted 
as evidence for the formation of a QGP. 

8. - Challenges and outlook 

In this last Section, we discuss observables that can help to better understand parton propaga- 
tion and fragmentation in QCD matter and, in particular, the relative role of perturbative parton 
energy loss (Section 6) and prehadron absorption (Section 7) mechanisms. We conclude with a 
review of proposed future facilities and experiments. 
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Fig. 63. - Results of the production-time r scaling analysis of HERMES data on n* production on Kr 
targets [22]. Left: z- and v- distributions as a function of t with corresponding as a function of the 
exponent A. Right: Scaling exponent and x~ per degree of freedom for "1ow-t" and "high-r" data sets 
compared to the full data set for pion production on Ne targets corresponding to "medium-r". Black disks: 
from preliminary HERMES data [383]. Red circles: from the AGMP absorption model [61, 326]. 
Blue triangles: n* from the energy loss model of Refs. [320, 326]. 



8'1. Partonic versus (pre)hadronic energy loss in nuclear DIS. - 

Despite very different assumptions on the parton lifetime, the models presented in Sec- 
tion 6 (resp. 7) based on partonic (resp. prehadronic) degrees of freedom, usually reproduce 
quite well the experimental results on the multiplicity ratio R'^^ in e + A collisions. It is therefore 
necessary to study observables which are more directly sensitive to the space-time evolution of 
the hadronisation process and can discriminate between partonic and pre-hadronic contributions. 
Candidates suggested in the literature include the atomic mass A dependence of the nuclear at- 
tenuation, its production time scaling, and the hadron transverse momentum broadening. 

8' 1 . 1 . Mass number dependence. - The mass number dependence of the nuclear attenuation 
ratio 1 - R'I^ had been suggested as an observable sensitive to the different mechanisms involved 
in the hadronisation process. To leading order in A^^^, one expects 1 - oc A^^^ in energy 
loss models because the average energy loss AEq oc (L^^ oc A^^^ for partons with asymptotic 
energies, due to the LPM interference in QCD [140]. On the other hand, in absorption models 
the survival probability is proportional to the amount of traversed matter, so that 1 -R'Ij oc {La) oc 
A^l^ [349-351]. However, it has been recently shown [61,321,363,382] that this simplistic 
argument does not really hold for absorption models, where additional powers of A'^^ need to 
be introduced to account for non-zero values of {tpreh), thus predicting R'^j oc A^^^ + G (A) as in 
energy loss models. Moreover, hadron attenuation is not likely to be describable by a simple A" 
power-law for nuclei with A > 80 [326]. 

81.2. Production time scaling in nDIS. - In Ref. [384] it is conjectured that Rm should depend 
on a scaling variable r defined as 
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rather than on Zh and v separately, i.e. 



(88) 



Rm = RM[TiZh,v)]. 



The scaling exponents A and /i can be obtained by a best fit analysis of data or theoretical compu- 
tations, whereas the proportionality constant C cannot be determined from the fit. Both absorp- 
tion and energy loss models can be distinguished by the value of the scaling exponent: a positive 
A value is characteristic of absorption models, which assume a short (tpreh) with a functional 
form similar to Eq. (87); while a zero or negative /I < is predicted within energy loss models 
due to energy conservation (a parton cannot radiate more energy than (1 - z/i)v). 

The fits performed in Ref. [384] and similarly in Ref. [23] (with yu = 1 taken from analytical 
Lund model computations [63, 64]) show that t is indeed a good scaling variable to analyze the 
HERMES data, see Fig. 63 (left), and furthermore it lives up to the promise of separating energy 
loss and absorption models. The charged pion data exhibit globally a clear A = 0.5± 0.15, which, 
being positive and different from zero, indicates that pre-hadronic effects significantly contribute 
to and allows to interpret t » (tpreh)- Moreover, it has been shown by applying different 
cuts on Zh and v that the A exponent derived from experimental data decreases with increasing t, 
thus suggesting an increasing contribution of partonic energy loss, see Fig. 63 (right). Finally, the 
sUght pion-charge and nuclear-mass number dependence of the A exponent observed in Ref. [384] 
may indicate additional dynamics beside production time and prehadron absorption. Fits to other 
hadron flavors, such as kaons, protons and antiprotons, are less conclusive due to the limited 
statistics. 

In conclusion, the r-scaling analysis of Rm seems a promising tool for studying energy-loss 
and formation-time contributions to hadron attenuation. However, new higher precision data sets 
in a wider range of y, mass number and hadron flavors are necessary to determine a possible 
scaling with /x^ I and to fully exploit this observation. 

8'1.3. pt -broadening and prehadron formation time. - The scaling analysis for Rm just de- 
scribed gives only indirect evidence of the pre-hadron formation time effects, but does not allow 
to measure the absolute scale of {tpreh)- A more sensitive observable to the pre-hadron formation 
time is the hadron transverse momentum broadening A(p^) in nDIS compared to a proton or 
deuteron target [29,70]. The pr-broadening is commonly believed to arise essentially during 
parton propagation because the prehadron is the quark momentum broadening A;?^ is propor- 
tional to the quark path-length in the nucleus [140], and the prehadron is supposed to have a 
negligible elastic cross section [29]. Under these assumptions, and remembering that hadron and 
parton transverse momenta are related by pt = zkj [385] with z~ Zh^A LO, we obtain 



with q the partonic transport coefficient, and the prehadron formation time obtained from the 
discussed scaling analysis: 



with C setting the overall scale of the production time to be determined from broadening data. 

The seemingly linear in A'^^ dependence at HERMES (Figure 21) leads to the conclusion 
that the prehadron is formed close to the surface or outside the heaviest nucleus, at least at the 
average kinematics (zh) = 0.41 and (v) = 14.0 GeV. Using an average in-medium path length 
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[GeV^] 


<v> 
[GeV] 


{z) 


2mjv<v> 


{tp) 
[fm] 


{I^pIi) vs a 
Ne (2.3 fm) 
Kr (3.7 fm) 
Xe (4.3 fm) 


2.4 
2.4 
2.4 


13.7 
13.9 
14.0 


0.42 
0.41 
0.41 


0.09 
0.09 
0.09 


4.2 
4.2 
4.3 


(APrft) vs z 


2.4 
2.4 
2.3 
2.2 


14.6 
13.3 
12.6 
10.8 


0.30 
0.53 
0.74 
0.92 


0.09 
0.10 
0.10 
0.11 


4.5 
3.7 
2.3 
0.7 




2.1 
2.5 
2.6 
2.4 


8.1 
12.0 
15.0 
18.6 


0.48 
0.42 

0.40 
0.36 


0.14 
0.11 

0.10 
0.07 


2.4 
3.7 
4.6 
5.8 


<Ap2,) vs 2^ 


1.4 
2.4 
4.5 


14.0 
14.1 
14.5 


0.41 
0.41 
0.39 


0.06 
0.10 
0.16 


4.2 
4.2 
4.3 



Table XI. - Average HERMES kinematics for the pi -broadening results [386, 387]. In parentheses (next to 
the target nucleus symbol) is the average in-medium path length of the hadronising system (L^i) » (3/4)i?^, 
withRA = (1.12/»2)A''''. The production time is computed according to Eqs. (89)-(90) with C = O.S fin/GeV. 
The average {xb) is very well approximated by {Q^}/{mt^{v}). 



{Lxe) ~ i3/4)Rxe = 4.3 fm in Eq. (90) one finds 

(91) C > 0.8 fm/GeV . 

The resulting prehadron production time is plotted in Fig. 64 left, and computed for each exper- 
imental bin in Table XL However, the size of the /7j--broadening requires q s; 0.03 GeV^/fm, 
which is much smaller than the q = 0.60 GeV^/fm obtained by fitting HERMES data on Rm with 
a pure energy loss model. Part of this discrepancy is due to the fact that the energy loss model 
used in the fit of a? is known to systematically overestimate the transport coefficient compared 
to other implementations [166]. However, this can hardly explain the factor ~20 discrepancy 
found here. One possible explanation is that (tpreh) is in fact smaller than just estimated. 

Turning to the Zh dependence (Figure 20), the observed decrease as z;, ^ 1 is qualitatively 
compatible with the production times just derived (because (tpre) as z;, — > 1), but also 
with long formation times and a pure energy loss scenario (because energy conservation requires 
AEq < (1 - z;,)v). However, in either case it is very hard to reproduce the detailed shape in 
Zh- Indeed, due to the zl factor in Eq. (89) in both cases (Ap^) would have a peak at Zh ^ 0.7, 
shifting to larger Zh values as A decreases. (Note that this behaviour of the peak with A is rather 
independent of the detailed shape of (tpreh)-) This is at variance with HERMES data, which show 
a fixed peak around Zh ~ 0.5, whose size increases with A. One scenario not considered so far in 
Uterature involves short formation times with large transport coefficients, to explain the shape in 
Zh, and a non-negligible prehadron elastic cross section to explain the A dependence 

On the other hand, the basically flat in v dependence of the broadening would be compatible 
with the production time estimate of Eq. (91). However, one should keep in mind that (x) de- 



PARTON PROPAGATION AND FRAGMENTATION IN QCD MATTER 



91 




0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 







- 


- ^ 




i 




ji* in Xe (HERMES preliminary) 






n in Xe (HERMES preliminary) 






71* in Xe calc. 


□ - 




jt" in Xe calc. 


1 



[GeV^] 



Fig. 64. - Left: Pre-pion formation time per unit v estimated in Ref. [384]. Right: hadron -broadening at 
HERMES as a function of compared to the computations of Ref. [390]. 



creases with v so that NLO contributions increase; since they tend to decrease the average energy 
of the fragmenting quarks they might in fact contribute to tilt down the slope. 

Finally, the linear increase of A{pj) with observed at HERMES (see e.g. Fig. 64) remains 
a challenge to current theoretical models. This dependence is in sharp contrast with the inverse 
power dependence from the colour dipole model of Ref. [29], thus ruling out the predicted f^rc/i 
EprehlQ^ behaviour due to the stronger gluon radiation. A flat predicted within most string- 
based models is also at variance with these data. Since a growing in prehadron formation 
time seems unlikely, dynamical effects at the parton level are needed in addition to prehadron 
absorption to explain the Q^-behaviour. Three mechanisms have been suggested in Ref. [388]: 

1 . Medium-enhanced DGLAP evolution. The longer medium-enhanced DGLAP evolution 
[303-305, 389] at larger would imply a stronger gluon radiation, hence a larger pT- 
broadening than at low Q^. 

2. Next-to-leading order processes. At NLO, gluon fusion process, y* -\- g ^ q + q, imply 
Eq < V, hence a smaller production time and p/- -broadening than at LO. The increasing 
importance of LO vs. NLO processes as xb = Q^/(2mi^v) increases could lead to the 
observed /57-broadening dependence on Q^. 

3. Prehadrons with very short production time, and inelastic cross section o'prehitp) 1/2^ 
slowly evolving in time. Under these conditions, they will be the less absorbed the larger 
Q^, suffer a smaller surface bias and contribute more to A{pj). 

Reference [388] also argues that studying the pr-broadening simultaneously binned in x and 
allows one to factor out the trivial v = /(Im^XB) kinematic coiTelations and to distinguish the 
proposed scenarios. 

The effect of medium modification of the DGLAP evolution has been explored quantitatively 
in Ref. [390], which proposes 

(92) MpIXQ^) = MpIXQ^) + zlvq 

The first term is the p^-broadening from quark multiple scatterings, computed at - 2.5 GeV^ 
using {tp} = 1.19z° ^'(l - z,,)i "V'^ from [61,64] and a transport coefficient ^ = 0.032 GeV^/fm 
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calculated in the dipole model for quark-nucleon scattering. The second term models the effects 
of the medium modified DGLAP evolution from Ref. [305]. Numerical results are presented in 
Fig. 64 right. Note that Eq. (92) predicts a plateau in Hp\) at 2^ > 4 GeV^, which should be 
experimentally verified. 

8' 1.4. A path forward. - As we have seen, HERMES data suggest relatively short production 
times, such that the prehadron is formed in average close to the nucleus surface, or even much 
more inside as the Zh dependence of the pr-broadening seems to require. The details of the v de- 
pendence of the /7r-broadening and its linear increase with demand for parton-level dynamics 
typically not included in theoretical models. 

To clarify this complex interplay between parton dynamics and hadronisation, one needs 
detailed theoretical computations, combining energy loss and prehadron absorption up to NLO, 
and accounting for surface bias effects mostly neglected in these simplified considerations. A 
production time scaling analysis of high precision Rm data coupled to a detailed study of the pj 
broadening dependence on all kinematic variables is likely to reveal the size of the production 
time and yield information on parton versus (pre)hadron dynamics. Experimentally, progress 
will be achieved by multi-dimensional binning, improved statistics for different hadron flavours, 
and using an extended set of targets up to Pb. 

The CLAS experiment is ideally suited for such a programme, but has a reduced range in 
V < 5 GeV, which enhances the role of prehadron absorption over parton energy loss, and in 
xb ^ 0.1. Its preUminary data follow the same quaUtative behaviour of HERMES data, but 
the measured pr-broadening is larger than what would be naively extrapolated from HERMES 
results based on the different average (v). This might be due to a significant prehadron elastic 
cross section, cr^'^^^ "c o-,^', additionally contributing to the hadron /jj-broadening at CLAS, but 
negligible at HERMES [390J . The origin of the CLAS/HERMES discrepancy needs clarification. 

Finally, the Electron-Ion Collider (EIC) [46,47], able to span from very small to large Xb, 
may clarify the role of LO and NLO processes. Furthermore, it will provide a large range in v, 
which wiU allow for the study of purely partonic in-medium propagation and the calibration of 
energy-loss models. 

8'2. ^r" vs. T] attenuation. 

\ lie measurement of jf' and rj suppression, which have a similar valence quark content but 
different masses and different hadronic cross sections, can provide an interesting cross-check 
of the hadronisation picture and time-scales for (pre)hadron formation. In Au + Au collisions 
at RHIC, the PHENIX experiment [147] has shown a similar suppression Raa{Pt) ratio for the 
two species up to pr ~ 10 GeV/c, with the n^/rj ratio independent of the collision system or 
centrality of the Au + Au collisions within uncertainties, see Fig. 30 (right). Such a result is 
naturally explained if the suppression takes place at the partonic level before the quenched parton 
fragments into final hadrons [146, 147], i.e., hadronisation occurs on timescales larger than the 
typical medium size. Yet, hadronisation in a deconfined QGP could be delayed compared to 
cold QCD matter at least until the medium temperature falls below the critical one. In fact 
results from e +A colhsions at HERMES, where the quark energy v « priRHic but the virtuality 

- 2.5 GeV^ ^2 _ ^(jqO GeV^), seem to favour much shorter hadronisation timescales 
for hard partons (Sect. S'l). 

It would be thus important to measure and compare n and tj attenuation in nDIS experiments, 
where clean measurements are possible and no QGP is formed. These measurements are planned 
for the CLAS experiment and will also be very important in a hadronisation physics programme 
at the proposed EIC. 
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8'3. Baryon formation . - 

Protons, and in general 3-quark baryon species, seem to follow a different nuclear modifica- 
tion pattern than mesons in both cold and hot nuclear matter. In nDIS, the difference observed 
in proton vs. mesons multipUcity ratios resembles the anomalous baryon enhancement reported 
at intermediate pr ~ 1.5-8 GeV/c in proton-nucleus and heavy-ion collisions. However, an- 
tiprotons follow the mesons attenuation pattern, unlikely what happens in p + A collisions where 
both baryons and antibaryons are anomalously enhanced. Contrasting the baryon nuclear modi- 
fications in nDIS and p+A colUsions will clearly lead to a more profound understanding of the 
hadronisation mechanism. 

8'4. High-pT hadrons at fixed-target energies. 

'1 lie detailed interpretation of the nucleus-nucleus and proton-nucleus data and the relative 

roles of high-pr hadron suppression and Cronin enhancement at around ^|sNN - 20 GeV is 
obscured by the current absence of a direct measurement of the corresponding p + p reference 
spectra at these energies. Use of various p + p interpolations [159, 391] compared with the A+A 
data measured at top SPS [160, 161] or low RHIC [176] energies, seems to indicate the onset 
of parton energy loss in this energy range. Precise tests of the mechanisms of pt broadening 
m p + A and of parton and/or (pre)hadron energy loss m A + A require the measurement of the 
baseUne vacuum production of various hadron species in a perturbative regime {pT <: 2 GeV/c) 
mp + p collisions at ^Jsm = 20 GeV. Such large statistics measurements can hopefully be made 
in coming low energy scans at RHIC (and at higher luminosity in the planned RHIC-II upgrade). 

8'5. The heavy flavour puzzle . - 

The solution of the intriguing heavy-flavour results at RHIC, namely the larger than expected 
suppression of heavy flavour production in central A+A colUsions, is likely to provide important 
insights in our theoretical understanding of parton propagation and hadronisation, and on the 
nature of the Quark-Gluon Plasma. 

On the theory side it has already prompted a revision of established ideas and stimulated 
many new ones, see Section 5'2. However, no consensus has been reached yet on the physical 
mechanism behind the large heavy-flavour suppression. Interestingly, heavy quark energy loss 
and in-medium propagation are processes that can be addressed in AdS/CFT-based descriptions 
of the QGP, and they are very promising to investigate the relevance of string-theory methods for 
discussing QCD near the deconfinement phase transition [392]. 

On the experimental side, new clues will be offered in the near future by experimental up- 
grades at RHIC, aimed at directly identifying D and B mesons and at measuring their individual 
suppression factors [194]. At the LHC, heavy-flavour production will be more abundant, and its 
nuclear modification one of the main areas of experimental interest because of their potential role 
in identifying and studying the QGP properties [393,394,418]. At the Electron-Ion Collider, a 
clean measurement of D and B suppression will also be able to discriminate the proposed mech- 
anisms, and remove the uncertainties in the interpretation of the data due to the modelUng of the 
medium created in A -i- A colUsions, and to possible in-medium regeneration processes. 

8'6. Dihadron and photon-hadron correlations. - 

High-pr two- and three-hadron correlations have proven a rich field to characterise medium 
modifications of parton fragmentation in A -i- A colUsions, see Section 5'3 . One can more directly 
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Fig. 65. - Left: A leading-order channel for y-TiP production. Right: y~n^ imbalance distributions using 
p^"' = 10 GeV/c and p^"' = 50 GeV. The solid lines show the (rescaled) quark fragmentation functions used 
in the calculation. 



access medium-modified FF in heavy-ion collisions by tagging large-/?^ particle production with 
prompt photons, as suggested in [281] and later investigated in detail at RHIC and LHC en- 
ergy [228, 395]. As shown in Fig. 65 left, at large pj and leading order in a,, a photon is 
produced directly in the hard subprocess, back-to-back to e.g. a quark which loses energy in the 
medium before fragmenting into a pion. Because of momentum conservation, the experimen- 
tal y-7r" momentum imbalance variable, z^^ = -p^^ ■ p^^ /IPt,P' reduces at LO to the theoretical 
hadron fractional momentum z which can therefore be experimentally estimated. In Fig. 65 right, 
the -y-TT^ distributions is computed in QCD at LO in p + p and Pb + Pb collisions at the LHC. The 
observable tracks the quark FF into pions until the onset of the photon-fragmentation channel at 
large in a wide range of z^^- This simple connection between z^^ and z, is however complicated 
by the effects of higher-order corrections or initial- and final-state soft gluon radiation. The most 
important background channel, to be reduced with appropriate kinematics and isolation cuts, is 
photon production by the collinear fragmentation of a hard parton. 

In contrast to heavy-ions, only dihadron correlations in Zh have been measured in e -i- A colli- 
sions, see Section 3'5. They are the analog of the "same-side" (intrajet) dihadron correlation in 
A -I- A integrated over the azimuthal angle, but are not affected by the complication of hot QCD 
matter effects. In the HERMES kinematics, (x) a; 0.1 where the gluon density is comparable to 
the quark density. Hence the processes y+q — > q + g and y + g^q + q may compete with the LO 
y + q^q scattering. For NLO processes, the smaller parton energy implies a larger nuclear sup- 
pression for either hadron, which in turn biases the hard scattering towards the nuclear surface. 
Moreover, the subleading hadron may be produced with larger z than the leading hadron, which 
complicates a naive interpretation of the observable. For these reasons, the relative importance 
of NLO compared to LO hadron production may be amplified in dihadron correlation compared 
to single-inclusive observables. An observable which is directly sensitive to ff{as) partonic cross 
sections, and could cross-check the dihadron Zh correlations, is the dihadron azimuthal correla- 
tion relative to the virtual photon direction. 

Photon-hadron correlations will be one of the key measurements at CLAS, and very interest- 
ing for a hadronisation programme at the EIC. Moreover, the photon-hadron transverse momen- 
tum imbalance with respect to the virtual photon axis will allow one to study the analog of the 
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Fig. 66. - Centre-of-mass energy vs. luminosity of past, present, approved and proposed lepton-nucleus 
facilities and experiments. 



same-side (intrajet) hadron correlations with the added advantage of measuring the transverse 
momentum of the jet. 

8'7. Future perspectives . - 

The challenges just discussed can be addressed and likely solved at future facilities planned 
(or proposed) to become operational within the next decade and beyond. The LHC heavy-ion 
programme, the RHIC-II upgrade, the 12-GeV upgrade at Jefferson Lab, and the Electron Ion 
Collider are examples of planned new facilities with capabilities that will significantly increase 
the existing experimental reach in terms of energy and luminosity, see Fig. 66. 

8'7.1. The 12-GeV Jefferson lab upgrade. - The 12-GeV accelerator upgrade of JLab will 
significantly extend existing experimental capabilities by doubling the accelerator electron-beam 
energy and by increasing the luminosity by at least an order of magnitude for large-acceptance 
fixed-target measurements. The CLAS12 spectrometer will operate with a continuous luminosity 
of 1-2x10^^ cm"^.s"' for nuclear target experiments with an 11 -GeV electron beam. The planned 
particle identification scheme, which will directly identify pions, low-energy charged kaons, 
neutral kaons, protons, and electrons, can be upgraded to include direct detection of higher- 
energy charged kaons as well. With the increased luminosity comes access to rarer hadrons 
whose attenuation in the nuclear medium has never before been explored. 

Table XII lists hadrons accessible with CLAS12 that are stable over nuclear distance scales. 
Nuclear attenuation measurements can be performed for all of these hadrons, and transverse mo- 
mentum broadening will be directly accessible for a number of the listed hadrons. Estimations 
for the CLAS12 geometric acceptances for these particles are plotted in Fig. 67. As can be seen, 
a systematic study of the mass dependence for two-quark systems is feasible for eight meson 
species with a range of masses spanning nearly 1 GeV/c^. Further, new insights into the poorly- 
understood hadronisation mechanisms for baryons can be gained by studying the series of eight 
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Table Xll. - Final-state hadnms potentially accessible for formation length and transverse momentum 
broadening studies in CLAS12. The rate estimates were obtained from the LEPTO event generator for an 
11-GeV incident electron beam. (The criteria for selection of these particles was that ct should be larger 
than the nuclear dimensions, and their decay channels should be measurable by CLAS12.) 



baryon species indicated in this table. While experimentally challenging, the baryon measure- 
ments have the potential to open up a completely new realm of studies of hadron formation for 
three-quark systems. 

87.2. Other fixed target experiments: DIS and DY. - The ongoing fixed-target COMPASS 
experiment [396] with its up to 190 GeV muon beam, can cover a much larger v and range 
than HERMES or JLab. In particular, due to the higher reach in v, one can select partons hadro- 
nising outside the medium, and cleanly test parton energy loss models by studying the hadron 
/77--broadening. Therefore, a yU + A experimental programme at COMPASS would be very de- 
sirable to improve our knowledge of the space-time evolution of hadronisation and anticipate 
some of the results expected at the future electron-ion colliders discussed in Sect. 87.4. At the 
same time the availability of a high-energy proton or pion beam would allow measuring nuclear 
modifications of DY processes and studying the interplay of initial state and final state parton 
energy loss within the same detector. 

Dedicated high-statistics medium energy DY experiments have been proposed in [318] to 
precisely measure initial-state parton energy loss in nuclear targets. The authors considered using 
the proton beams from the 120 GeV Fermilab Main Ring Injector and the future 50 GeV Japan 
Hadron Facility (JPARC). The kinematic acceptance of a dimuon spectrometer proposed for the 
Fermilab Experiment E906 [397] and a similar spectrometer considered [398, 399] at JPARC. 
The nuclear partons are probed in the 0.1 < x < 0.5 region where the anti- shadowing and the 
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Fig. 67. - Acceptances for mesons and baryons in the CLAS12 experiment at the proposed 12-GeV JLab 
facihty. 



EMC effect on nuclear PDFs interferes with energy loss suppression of DY events. However, in 
this jc-region nuclear PDFs are well determined in global fits [400-402]. 

The interplay of hadron attenuation and /97-broadening in nDIS and DY nuclear modifications 
in /;> + A collisions will be vital to constrain energy-loss models and characterise the stopping 
power of cold nuclear matter. 

8'7.3. Other fixed target experiments: high-px hadron spectra. - The study of high-/?7- hadron 
spectra in p + p to A + A collisions in the SPS energy range is fundamental for precisely detect 
and study QGP formation in heavy-ion collisions. Such a programme is one of the aims of 
the NA61 experiments at CERN [403], which upgrades the NA49 spectrometer, and also by 
the NA60 [404] and SELEX [405] experiments, originally focused on quarkonium and charmed 
hadron production and nuclear modifications. 

The peculiarity of NA60 and SELEX is that they mount multiple nuclear targets on the same 
beam: from Beryllium to Uranium in NA60; Copper and Carbon in SELEX. As emphasised 
in [162] this target setup allows a high-precision measurement of the ratio of differential cross- 
sections for inclusive hadron production on two different nuclei, R'^^ipB ~ J dp^dy l ^ d^dy ^^^^"^ 
B is the lightest available nucleus. Indeed, the beam luminosity as well as part of the systematic 
errors cancel in the ratio. The A dependence of the ratio is analogous to the centrality dependence, 
but eliminates the large experimental uncertainties due to the determination of the centrality and 
to the normalisation of the Cronin ratio. Moreover, without need of centrality cuts the statistics 
may be sufficient to probe the high-/9r region and test the mechanism underlying the Cronin 
effect. 

8'7.4. Electron-ion colliders. - An electron-ion collider would significantly complement and 
add new dimensions to the experimental studies at fixed target facilities, by extending the range 
of accessible lepton, E^, and ion, Ea, energies. Proposed designs for an electron-ion col- 
lider currently have a fairly wide range of luminosities and centre-of-mass energies, see Ta- 
ble XIII. The Electron-Ion Collider (EIC) is a proposed US-based facility, which would provide 
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Table XIII. - Beam energies E, center of mass energies ^fs, and luminosity for e + A collisions at the 
ENC (at FAIR), the EIC (BNL and JLab designs), and the LHeC ( ring-ring mode) proposed facilities. 



Ee - 3 - 20 GeV and Ea - 15-100 AGeV [46,47]. There are currently two complementary 
conceptual proposals [47,406]: (i) eRHIC would add an electron beam to the existing RHIC 
ion accelerator at BISIL - the option of staging it, and first realise a lower energy version with 
Eg -2-4 GeV, based on an energy recovery linear accelerator for the electron beam is being ac- 
tively discussed [407,408]; (ii) EIC @ JLab would add an ion accelerator to the CEBAF upgraded 
12-GeV electron ring at Jefferson Lab [409,410]. The eRHIC concept emphasises the energy 
range; the EIC@JLab concept will reach up to 100 times more luminosity, but with lower max- 
imum energy. In both cases polarised proton and light-ion beams will be available. The Large 
Hadron-electron Collider (LHeC) is a proposed upgrade of the LHC at CERN [411,412], and 
will reach much higher energies than the EIC still with a good luminosity, but only with unpo- 
larised hadron beams. Finally, the idea of building a low-energy electron-nucleon (ENC) collider 
at FAIR has recently been advanced [413]: the goal is to have a centre-of-mass energy between 
HERMES and COMR\SS, but up to a factor 100 higher luminosity and polarised hadron beams. 
Here, we will focus on the EIC capabilities, but most of the discussed measurements will be in 
principle possible also at the LHeC and, with the exception of jet measurements, at the ENC. 
The EIC will provide virtual photons with energies in the v = 10 - 1600 GeV range, large 
up to 1000 GeV^, and low x > 10"^. Its high luminosity will allow access to rare signals, 
multi-differential measurements, and dihadron and y-hadron correlations. It will be able to cross- 
check HERMES and CLAS measurements, while offering many more channels and an extended 
kinematic range to study hadronisation inside the nucleus at low v, as well as testing basic QCD 
processes such as non-Abelian parton energy loss and the space-time evolution of the DGLAP 
shower at high v. EIC will be unique, compared to HERA (HERMES) and JLab (CLAS), in the 
foUowing areas: 

• At the large v accessible at the EIC, hadrons will clearly be formed outside the nuclear 
medium, so that effects due to the propagation of the struck quark can be experimentally 
isolated. One wUl have new access to pr-broadening studies, which can cleanly probe 
the parton radiative energy loss as predicted by pQCD at asymptotic energy. It will also 
be possible to study in detail the interplay of radiative and collisional parton energy loss, 
medium modifications of the DGLAP evolution, and test factorisation for the fragmenta- 
tion functions. 

• Heavier mesons like rj and will be more abundantly produced than at HERMES/JLab. 
At both medium-low and large v, comparing their quenching to that of their lighter n and 
K counterparts will provide important clues about parton propagation and hadronisation, 
see Sect. 8"2. 

• The EIC excellent low-x coverage provides increased production of heavy flavour and 
quarkonia. In particular, it will be possible to study the heavy quarks energy loss through 
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the B and D meson suppression. Theoretical mechanisms proposed to explain charmonium 
suppression in p + A and A +A collisions can be put to the test in a clean experimental 
envirormient by studying J/i/f, iff' mdx spectra. 

• The large coverage will access a truly perturbative QCD regime, whose prediction can 
be more confidently tested against the data, in particular colour transparency eff'ects and the 
2^-dependence of the observables discussed in this review, like the hadron pj-broadening. 

• Baryon production through parton fragmentation will also be accessible, because of the 
collider mode and the accessible final-state invariant masses. This will allow studying 
baryon transport in cold QCD matter and the baryon anomalies observed in fixed target 
e +A collisions at HERMES and in heavy-ion collisions at RHIC. The ability to identify 
a good variety of baryons, including the strange and charmed sector, wiU be a key to this 
programme. 

• For the first time, jet physics in DIS with a nuclear target will be experimentally accessible. 
In particular, medium modifications of the jet shape, and the comparison of light-quark to 
heavy-quark and gluon initiated jets will shed light on the mechanisms underlying parton 
energy loss. These studies can also be extended to dijet or y-jet correlations. 

In summary, the collider kinematics and associated detectors with excellent calorimetry, particle 
identification, vertex detection, and rapidity coverage will allow for a comprehensive programme 
to better understand how parton energy loss and fragmentation occur in cold QCD matter. 

By contrast, the ENC will cover a smaller the y = 10-50 GeV range, which might or 
might not be sufiicient to ensure hadronisation outside the nuclear medium. The high luminosity 
and rather symmetric kinematics would however ensure the possibility of a rich hadronisation 
programme complementing and extending the HERMES results. 

8'7.5. RHIC II . - In order to undertake extensive high-statistics studies of processes with low 
cross sections in A -i- A collisions, an upgrade of the RHIC luminosity will be necessary as well as 
a comprehensive new detectors [414]. Key measurements will include two-particle correlations 
up to a pt range where fragmentation can be safely tested and extended to include heavy iden- 
tified hadrons. In addition, a comprehensive set of measurements of photon-jet correlations, and 
heavy-quark tagged jets can provide valuable information on parton energy loss and hadronisa- 
tion mechanisms. The capabilities of the proposed new detector include: (i) excellent charged 
particle momentum resolution up to pr - 40 GeV/c in the central rapidity region, (ii) complete 
hadronic and electromagnetic calorimetry over a large phase space, (iii) particle identification up 
to large pr (20 - 30 GeV/c) including hadron and lepton separation in the central and forward 
regions, and (iv) high rate detectors, data acquisition, and trigger capabilities. 

87.6. LHC . - The possibility of colliding heavy nuclei with high luminosity at the Large 
Hadron Collider (LHC) ofl'ers a unique opportunity to investigate in detail the behaviour of 
strongly interacting matter under extreme conditions of density and temperature. The factor 
30 increase in -\/snn relative to RHIC corresponds to a huge increase in kinematic and statistical 
reach for hard probes, and new channels wiU become available [172] such as detailed jet studies 
(jet shapes, medium-modified fragmentation functions) [137]. 

ALICE: - The ALICE experiment was designed specifically to tackle as many measurements 
as possible in the high multiplicity environment of heavy-ion collisions. ALICE aims firstly at 
accumulating sufficient integrated luminosity in Pb + Pb collisions at ^[sm =5.5 TeV, as well 
as to carry out studies oip+p and p + A collisions in order to establish the benchmark processes 
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under the same experimental conditions. ALICE can measure heavy-flavour production down to 
very low pt because its unique low transverse-momentum cutofl" for particle detection. This can 
be achieved by using inclusive large impact-parameter lepton detection, and by reconstructing 
exclusive charm-meson decays at relatively low pj. As a result, the measurement of the total 
heavy flavour cross section will require a smaller extrapolation and will thus show improved pre- 
cision. As discussed in 8'6 a very attractive methods to study hadron jets is to tag high energy 
hadrons or jets with prompt photons emitted opposite to the hadrons or jet directions. The com- 
bined use of the PHOton Spectrometer (PHOS) - with excellent energy resolution although in a 
limited acceptance - and the ElectroMagnetic Calorimeter (EMCal) - about seven times larger in 
acceptance but with slightly worse resolution and coarser granularity - and the central tracking 
system for charged particles, will allow one to recover a large fraction of the jet energy, thus 
reducing sensitivity to the specific pattern of fragmentation, and to perform y-hadron and y-jet 
correlations measurements. 

ATLAS: - ATLAS is the largest particle detector ever constructed, and its design, like CMS, 
is aimed at addressing a broad variety of physics processes [415]. The detector design has many 
features that make it ideal for studies of heavy-ion collisions [416]. A central solenoid magnet 
and large outer toroidal field magnets provide momentum analysis for the tracking systems. The 
large acceptance electromagnetic and hadronic calorimeters have longitudinal segmentation and 
fine transverse segmentation, covering the range in pseudorapidity I77I < 3.2 (electromagnetic) 
and \r]\ < 4.9 (hadronic), providing excellent photon and jet physics coverage. The high preci- 
sion siUcon tracking system covers the range \ri\ < 2.5 and the external muon spectrometers cover 
\ri\ < 2.7, for large-acceptance dimuon studies of heavy quark systems. A fast, three-level trig- 
ger system and a high-rate data acquisition system are designed for triggering on rare, high-pr 
particles and jets. 

This instrumentation provides excellent tools for study of several observables relevant to 
heavy-ion coUisions. Study of hard interactions and jet quenching will be a strong focus, and 
are closely connected to the topics in this review. The acoplanarity of di-jet pairs should provide 
sensitivity to the angular diffusion ofhigh-pj partons in the medium. Jet properties that will be 
measured include transverse momentum spectra, hadron fragmentation functions, and jet shapes; 
di-jet angle and energy correlations may be able to separately quantify coUisional and radiative 
energy loss of hard-scattered partons. In addition, measurements of quarkonia in bottom and 
charm systems, as well as of D and B mesons, will open up detailed studies of heavy quark jet 
quenching at high Ej. 

CMS: - CMS is a general purpose experiment at the LHC designed to explore the physics at 
the TeV energy scale [417] including comprehensive measurements in Pb + Pb collisions [418]. 
The CMS detector features a 4 T solenoid surrounding central silicon pixel and microstrip track- 
ing detectors and electromagnetic (ECAL, I77I < 3) and hadronic (HCAL I77I < 5) calorimeters, 
and muon detectors (|?7| < 2.4) embedded in the flux return iron yoke of the magnet. CMS is 
the largest acceptance detector at the LHC with unique forward detection capabilities with the 
HE and CASTOR calorimeters (5.1 < I77I < 6.6). CMS is, by design, very well adapted to de- 
tect and reconstruct high-pj and high-mass particles. The experiment can significantly extend 
the Pt reach with respect to RHIC, thanks to the large hard cross sections at 5.5 TeV, the large 
acceptance of its tracking system ([//I < 2.5), and its high-pT- triggering capabilities. The silicon 
trackers have excellent reconstruction performances starting above pr > 0.3 GeV/c with a low 
fake track rates in central Pb + Pb. Within \r]\ <1, the pj resolution is better than 2% up to the 
highest Pt values reachable. The leading hadron suppression can thus be measured with low 
uncertainties aU the way up to 300 GeV/c, allowing us to clearly discriminate between various 
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model predictions (Fig. 32, left). 

Full jet reconstruction in Pb + Pb collisions can be performed in CMS using the ECAL and 
HCAL calorimeters. Jets start to be distinguishable above the background at Ej ~30 GeV and 
can be fully reconstructed above 75 GeV (efficiency and purity close to 100%) with a good 
energy resolution (better than 15%). The expected jet Ej spectrum after one month of Pb + Pb 
running (0.5 nb"') - taking into account High Level Trigger bandwidths and quenching effects as 
implemented in HYDJET [3 14] - reaches up to Ej ~ 0.5 TeV in central Pb + Pb, and will give 
access to detailed differential studies of jet quenching phenomena (jet shapes, energy-particle 
flows within the jet, ...). The possibility of CMS not only to fully reconstruct jets, but to tag them 
with prompt y and/or Z bosons [419] and to carry out detailed studies in the c, b quark sector 
will be very valuable to clarify the response of strongly interacting matter to fast heavy-quarks, 
and will provide accurate information on the transport properties of QCD matter. 



9. - Conclusions 

The physics of parton propagation, interaction and fragmentation in both cold and hot strongly 
interacting matter has been reviewed. The most recent theoretical and experimental results 
on hadron production in deep inelastic lepton-nucleus scattering (nDIS), high-energy hadron- 
nucleus (h+A) and heavy-ion (A + A) collisions as well as Drell-Yan processes in hadron-nucleus 
interactions have been discussed. The main motivation of such studies is to provide new insights 
on the mechanisms of parton fragmentation and on the space-time evolution of quark and gluon 
hadronisation. Those are complex dynamical processes characterised at some stage by energy 
scales for which the QCD coupling becomes large and thus cannot be completely addressed with 
standard perturbative QCD tools. 

The operational approach consists in comparing the modifications in hadron production ob- 
served in a "cold QCD matter" environment (as found in nDIS oi p + A collisions) as well as in 
"hot QCD matter" (created in A + A colUsions) with respect to the measurements in the "QCD 
vacuum" (in DIS off protons or in proton-proton collisions). The quantitative comparisons are 
carried out using various phenomenological ratios (hadron multiphcity ratio in nuclear DIS, 
or nuclear modification factor RpAM in nuclear collisions) as a function of the relevant variables: 
hadron momentum, fraction of the photon energy, virtuality, hadron flavour, size of the medium, 
etc. In the case of cold nuclear matter, the detailed study of the observed deviations with respect 
to vacuum production allows one to constrain various models describing the space-time evolu- 
tion of parton fragmentation and, in the case of nucleus-nucleus collisions, also to characterise 
the (thermo)dynamical properties of the hot and dense QCD matter produced in the reaction. 

A detailed comparison of the kinematics and phase-spaces in nDIS experiments at CERN, 
HERA and JLab and in proton-nucleus (FNAL, SPS, RHIC) and heavy-ions (SPS, RHIC) col- 
lisions, indicates that the range of relevant hadron energies measured is comparable, £/, » 2- 
20 GeV, as well as typical values of fractional momentum of the parent parton energy carried 
away by the leading hadron, z ~ 0.4-0.9. The same is not true, however, for the parton virtual- 
ities Q involved which are at least a factor of five larger for the hadrons produced at RHIC top 
energies (Q < 20 GeV) than those typical at HERMES/JLab (g < 3 GeV). Therefore a direct 
comparison of hadron suppression data from nDIS to h + A and A + A colUsions is not possible, 
and theoretical models are needed to connect both colliding systems. 

The main theoretical approaches to account for parton propagation and fragmentation in QCD 
matter can be roughly divided in two groups: (i) models rooted in perturbative-QCD studies of 
non-Abelian parton energy loss via gluon Bremsstrahlung or in-medium modifications of the 
DGLAP equations for gluon radiation, followed by hadronisation of the final state hadrons in 
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the vacuum, and (ii) approaches mostly based on modifications of the non-perturbative Lund 
string fragmentation model that account for the interaction in the nuclear medium of a colorless 
pre-hadronic stage. 

The range of application of these two classes of models depends on the value of the parton 
lifetime, i.e., by the time tpreh needed to produce a prehadron. Simple perturbative estimates 
indicate that the typical parton lifetime at HERMES or RHIC strongly depends on the hadron 
mass. For pions, one obtains lifetimes tpreh ~ 20- 30 fm, larger than the nuclear radius, al- 
though modified versions of the Lund and dipole fragmentation models indicate smaller life- 
times, tpreh ^ 6 fm. For heavier particles such as K, rj, p, one finds tpreh ~ 6 - 9 fm, and even 
smaller values for heavy-flavor mesons, tprehiB, D) < 1 fm. Therefore it is in general necessary 
to incorporate prehadron-medium interactions in the theoretical models, and to look for experi- 
mental observables able to distinguish them from purely partonic interactions and measure tpreh- 

A simple way to test the different models and their underlying assumptions is to study hadron 
production in nDIS and Drell-Yan processes inh+A collisions, where the traversed medium is 
static, if not uniform, and its properties are well known. This cold nuclear matter information 
serves then as baseline for A + A collisions, where a number of new mechanisms modifying 
high-pT hadron production can be used to "tomographicaUy" study the thermodynamical and 
transport properties of the hot and dense matter created. 

Lepton-nucleus DIS 

The recent HERMES and older EMC data on hadron attenuation in nDIS, encoded in the 
v; A) ratio, support a picture where the space-time evolution of the fragmentation process is 
modified by the surrounding matter, with a sizeable inelastic cross section between the colorless 
pre-hadronic state and the nuclear medium. However, when the detailed geometry and finite size 
of the target nucleus are taken into account, nDIS data can also be described in terms of parton 
energy loss alone. 

Additional observables sensitive to the underlying dynamics of the hadronisation process 
have been investigated. In particular, the observed slow variation of Rm with disfavours 
mechanisms Uke partial deconfinement and puts additional constraints on other theoretical mod- 
els. The rising and seemingly linear behavior of the hadron pr-broadening (Apr) vs. A at HER- 
MES supports the conclusion that the prehadron is formed close to or just outside the nucleus 
surface (at (v) =s 14 GeV and ^ 2.5 GeV^) and allows to estimate the prehadron production 
time. The details of the v dependence of the pr-broadening, and especially its linear increase 
with are better interpreted in terms of medium-modified DGLAP evolution, and demand for 
NLO computations so far neglected. However, the shape of (Apr) vs. Zh cannot be qualitatively 
understood by current prehadron absorption or energy loss computations, and is challenging the- 
ory models. 

Hadwn-nucleus collisions 

The observation of a Cronin enhancement {RpA > 1) of the hadron yields in the range 
Pt =2-S GeV/c is a characteristic feature of high-energy proton-nucleus collisions at fixed tar- 
get energies (FNAL, SPS, HERA-B). At RHIC top energies, intermediate-p^ proton, but barely 
meson, production is enhanced in d+Au collisions and the typical iJ^A maximum around 2 GeV/c 
is progressively suppressed as the rapidity increases. The usual interpretation of the Cronin ef- 
fect is based on multiple scatterings of the parton prior to its fragmentation. Both initial- and 
final-state interactions must play a role in order to explain the various experimental observations. 
On the one hand, the disappearance of the enhancement at forward rapidities is easily explained 
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by non-linear QCD evolution of the initial-state parton distributions (though may also be due 
to effects occurring at the edge of phase space). On the other hand, DY data suggest that par- 
ton rescatterings contribute only a small amount to the hadron enhancement and, in addition, 
final-state coalescence of the scattered parton with other soft partons in the nuclear environment 
naturally explains the larger enhancement observed for baryons compared to mesons. The Cronin 
efi^ect helps thus to study the hadronisation mechanism in hadronic interactions at intermediate 
Pt, which may be dominated by parton recombination rather than by independent parton frag- 
mentation. Detailed control of the role of all these mechanisms is important to help identify the 
onset of high-/7r hadron suppression in A + A colUsions in the = 10 - 100 GeV range. 

Nucleus-nucleus collisions 

The two most notable experimental results in Au + Au collisions at RHIC are (i) the observed 
factor of ~5 suppression of high-pj hadrons, and (ii) the strongly distorted azimuthal distribu- 
tions of secondary hadrons emitted in the away-side hemisphere of a high-;??- trigger hadron. 
Their properties such as magnitude and light flavour "universality", and their dependence on 
Pt, centrality, path-length, and ^/Jnn are in quantitative agreement with the predictions of mod- 
els based on non-Abelian gluon radiation ofi" hard scattered partons traversing the dense QCD 
medium. The confrontation of these models to the data permits to derive the initial gluon density 
dN« Idy =s 1400 and transport coeflicient q = GeV^/fm) of the produced medium at RfflC. 
Yet, other observations such as the same suppression factor for Ught- and heavy-quark mesons, 
and the less suppressed baryon compared to meson production indicate that the standard fac- 
torisation assumption of vacuum hadronisation after in-medium radiation, implicit in all parton 
energy loss formaUsms, may well not hold. At lower collision energy, with shorter propagation 
times and a less dense final state medium, hadron suppression in the cold nuclear targets may 
compete with that originating from the hot medium. 

Outlook 

We have reviewed a set of observables and new or improved measurements that can answer 
several of the open questions mentioned here and discussed in detail in this review. They in- 
clude multi-dimensional pj- -broadening measurements, hadron-hadron and hadron-photon cor- 
relations, heavy flavor and jet shape modifications in nDIS, improved DY measurements and 
large-pr hadron spectra mh + A colUsions, heavy quark tagging in A -i- A colUsions. These ob- 
servables can be addressed at new facilities planned for the near or medium-term future, such 
as RHIC-II, JLabl2, and the Electron-Ion Collider (which in particular will open the study of 
purely partonic in-medium processes in DIS). They can also be addressed by new experiments 
at existing facilities, such NA61 at SPS, or by a creative use of data taken at closed experiments 
such as NA60 or SELEX, among others. The combination of new results, theoretical develop- 
ments, and an open mind to combine information from traditionally different fields such as DIS 
and hadronic collisions will doubtlessly lead to a more profound knowledge of parton propaga- 
tion and fragmentation, quark and gluon hadronisation, and the properties of strongly interacting 
QCD matter. 
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